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FOREWORD 


NASA TM-82478 and NASA TM-82501 (Volumes 1 and 2) entitled, "Space and 
Planetary Environment Criteria Guidelines for Use in Space Vehicle Development, 
1982 Revision," provide information relative to the natural environment for 
altitudes greater than 90 km above the surface. (NASA TM-82473 entitled, 
"Terrestrial Environment (Climatic) Criteria Guidelines for Use in Aerospace 
Vehicle Development, 1982 Rivision," provides natural environmental informa- 
tion for altitudes below 90 km.) NASA TM-82478 and NASA TM-82501 are a 
compendium of up-to-date information on the physical characteristics of bodies 
in the solar system and a guide to recent scientific literature. They are 
intended for NASA engineers, for scientists who design spacecraft and space 
experiments, and for teachers and others interested in the space program. 

Figures and tables reprinted from copyrighted material are included with 
permission of the copyright holder. 

These volumes are a follow-on to NASA TM-78119, 1977. However, NASA 
TM-78119 is still considered to be of use and of interest to scientists and 
engineers in some areas, but any conflict in the data will be resolved by 
using the values in the new documents. Furthermore, there is no intent to 
automatically change any references to previous documents in contract Scopes 
of Work by the issuance and acquisition of either TM-82473, TM-82478, or 
TM-82501. 

These guidelines are recommended as a tool for use in the development 
of space vehicles and associated equipment, and the information presented is 
based on data and models currently considered to be accurate. However, an 
environment specialist should be consulted in critical design interface appli 
cations for the most current information and scientific interpretation. 

Sections on the Sun, Terrestrial Space, the Moon, Mercury, Venus, and 
Mars constitute Volume I. Volume II contains chapters on Jupiter, Saturn, 
Uranus, Neptune, Pluto, Comets, Asteroids, and Interplanetary Dust. 

Various programs of NASA's Office for Space Transportation Systems, 

Office for Aeronautics and Space Technology, Office for Space Science and 
Applications, and Office for Space Transportation Operations provided resources 
required for the preparation of this document. 
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SUMMARY 


This document provides a consolidated presentation of space and planetary 
natural environment data for use as design criteria guidelines in space vehicle develop- 
ment programs. 

Specifically, information is provided in the disciplinary areas of atmospheric and 
ionospheric properties, radiation, geomagnetic field, astrodynamic constants, and 
meteoroids for the Earth's atmosphere above 90 hm, interplanetary space, and the 
atmospheres and surfaces (when available) of the Moon and the planets (other than 
Earth) of this solar system. The current Marshall Space Flight Center (MSFC) upper 
atmosphere model is described in detail. 


INTRODUCTION 


A knowledge of environment parameters is necessary for the establishment of 
design requirements for space vehicles and associated equipment. Such data are 
required to define the design condition for fabrication, storage, transportation, test, 
preflight , in-flight, and on-orbit design conditions and should be considered for both 
the whole system and the components which make up the system. The purpose of 
this document is to provide guideline data on space and planetary environmental con- 
ditions relative to various disciplinary areas which are applicable to the design of 
space vehicles and associated equipment for NASA. 

Good engineering judgment must be exercised in the application of the environ- 
ment data to space vehicle design analysis. Consideration must be given to the over- 
all vehicle mission and performance requirements. Knowledge still is lacking on the 
relationships between some of the environment variates which are required as inputs 
to the design of space vehicles. Also, interrelationships between space vehicle 
parameters and enAuronment variables cannot always be clearly defined. Therefore, 
a close working relationship and team philosophy should exist between the design/ 
operational engineer and the respective organization's environment scientists. Although 
a space vehicle design should accommodate all expected operational environment condi- 
tions, it is neither economically nor technically feasible to design space vehicles to 
withstand all extremes. For this reason, consideration should be given to protection 
of space vehicles or critical subsystem for some extremes. This document does not 
specify how the designer should use the data in regard to a specific space vehicle 
design. Such specifications may be established only through analysis and study of a 
particular design problem. 

Assessment of the natural environment dn early stages of a space vehicle develop- 
ment program will be advantageous in developing a space vehicle with a minimum opera- 
tional sensitivity to the environment. For those areas of the environment that need 
to be monitored prior to and during tests and operations, this early planning will 
permit development of the required measuring and communication systems for accurate 
and timely monitoring of the environment. 

The environment criteria data presented in this document were formulated based 
on discussions and requests from engineers involved in space vehicle development and 
operations; therefore, they represent responses to actual engineering problems and 
are not just a general compilation of environmental data. This report is used exten- 
sively by various government and private space vehicle development organizations in 
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design and operational studies. Inquiries may be directed through appropriate organ- 
izational channels for subsequent communications to the Atmospheric Sciences Division, 
Marshall Space Flight Center, Alabama 35812. 

Based on known and projected user requirements, it was decided that it is 
more advantageous to present the natural environment parameters grouped according 
to spatial regions rather than scientific disciplinary areas. These spatial areas are 
interplanetary space and those surrounding the individual planets and the Earth's 
Moon. The outer limits depend upon the component of the environment being dis- 
cussed. For clarity in disciassion, limits have been arbitrarily taken to be 10 radii 
above the surface of the planet or Moon. In the discussion of the meteoroid environ- 
ment, however, the limits extend to a point where the gravitational attraction of the 
planet or Moon becomes negligible. 

Some of the space and planetary environmental data available today are specu- 
lative and will remain so until additional satellites , planetary probes , and manned 
space flights provide more information. Therefore, care must be used in the inter- 
pretation and use of available space environment data for a specific design decision. 
Although the data in this document provide valuable guides for preliminary design 
studies and analyses, their use for final program decisions depends upon the specific 
design problem Involved. This report was prepared primarily for use in NASA space 
vehicle development and advanced study projects. The data contained in this docu- 
ment are reviewed on a continuing basis, and revision or amendments will be pub- 
lished as necessary. The numbers given in brackets refer to references given at the 
end of each section. 

This document does not include the subject of environmental test procedures. 
Reference should be made to Department of Defense MIL-STD-8ipc Environmental Test 
Methods (1975)* available from the National Technical Information Service, Springfield, 
Virginia, 22161. This document covers procedures for: Low Pressure (Altitude), 

High and Low Temperature, Temperature Shock, Temperature Altitude, Solar Radiation, 
and Space Simulations (Unmanned Test). An excellent comparison of the various 
international environmental testing standards may be found in the Journal of Environ- 
mental Sciences, Vol. XXIV, Number 2, March/April 1981. 


*Revision being considered by DOD. 



TECHNICAL MEMORANDUM 82478 


SPACE AND PLANETARY ENVIRONMENTCRI HR! A GUIDELINES FOR 
USE IN SPACE VEHICLE DEVELOPMENT, 1982 REVISION 

(VOLUME I ) 

SECTION 1. THE SUN 


1.1 INTRODUCTION 

The Sun is the star about which all bodies in our solar system revolve. It 
provides the light and heat upon which we depend for life, and its gravity field 
determines the movement of the planets and other bodies in their orbits. It is a main 
sequence star of type G-2. In any other place in the universe, it would scarcely be 
worthy of special study, but its central role in our solar system makes it a unique 
object for research. 

Short and long wavelengths of solar electromagnetic radiation along with ener- 
getic particles emitted by the Sun are absorbed by the Earth's atmospherev Moreover, 
the extremely bright central solar disk interferes with observation of the corona except 
at times of solar eclipse or, less effectively, when the coronagraph is used to block 
the central image artificially. Much new information has been learned in the past 
several decades through observations from platforms above the atmosphere: balloons, 

high-flying aircraft, rockets, satellites, and spacecraft. Besides observations of 
electromagnetic radiation at all wavelengths, such platforms have permitted direct 
measurements of solar particle emissions that had previously been inferred from geo- 
magnetic and other effects observable at the surface. 

Like all main sequence stars, the Sun's energy output derives from the nuclear 
fusion conversion of hydrogen to heliiim and heavier nuclei in its core. Photons 
diffuse outward until, after several thousand years, they reack the convective 
envelope (Fig. 1-1). Another byproduct of the fusion reaction is a flux of neutrinos 
which freely escape the Sun. (Reference 1-1 discusses attempts to detect the neu- 
trino flux and interpretations of the puzzling results. ) 


1.2 BRIEF QUALITATIVE DESCRIPTION 

At the. center of the Sun is the core, where energy is produced by nuclear 
reactions. Surrounding the core are concentric shells: the radiative zone;, the con- 

vective envelope at about 0.8 Rg ; the photosphere, the layer from which visible radia- 
tion emerges; the chromosphere ; and the corona, the Sun's outer atmosphere. Energy 
is transported outward through the convective envelope by convective motions which 
are organized into cells. The Sun's "lower atmosphere" or photosphere is the region 
from which energy is radiated directly into space. This radiation approximates the 
Planck distribution at 5800°K; most of the energy falls within the visible part of the 
spectrum (Fig. 1-2). 


i-1 



1-2 



Figvire 1-1. A cross section of the Sun indicates the inner energy - generating 
core (0. 2 in radius) , the radiative envelope which extends out to 

0.8 R , and the outer convective shell [1-12] . 
o 
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Figlire 1-2. Normally incident solar radiation at sea level on very clear days, solar spectral irradiance 
outside the Earth's atmosphere at 1 AU, and blackbody spectral irradiance curve at T = 5762°K 

(normalized to 1 AU) [1-2]. 
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The chromosphere, extending for a few thousand kilometers above the photo- 

4 

sphere, has a maximum temperature of around 10 °K. The corona, which extends 

0 

several solar radii above the chromosphere, has temperatures of some 2 x 10 °K. 

These regions emit electromagnetic radiation in the UV, EUV, and X-ray spectral 
regions. This radiation, though of relatively small total energy, is the major factor 
in forming ionospheres and in planetary upper atmosphere photochemistry. 

A consequence of the heating of the Sun's outer atmosphere is its expansion into 
the surrounding interplanetary space. This continuous outward flow of plasma is 
called the solar wind. It consists of protons, electrons, and alpha particles, as well 
as small fluxes of heavier ions. Typical velocities are 300 to 400 km /sec, but its 
velocity may be as high as 1000 kmVsec. 

Although the total energy output of the Sun is remarkably steady, the surface 
of the Sun exhibits many types of irregularities, including sunspots, faculae, plages 
(bright areas), filaments, prominences, and flares. All are believed to be ultimately 
the result of interactions between ionized gases in the solar atmosphere and the Sun's 
magnetic fields. Most solar activity follows the sunspot cycle. The appearances of 
individual sunspots on the solar disk cannot be predicted, but they statistically follow 
a fairly regular periodic pattern (Fig. 1-3). Sunspots are discussed in greater detail 
in the next section. 



Figure 1-3. The Sunspot Number Chronology as derived from direct measurements 
and historical documents. (Figure courtesy of Dr. J. A. Eddy, 

High Altitude Observatory.) 


1.3 PHYSICAL PROPERTIES 

Table 1-1 lists the nhvsical properties of the SilH. 
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TABLE 1-1. PHYSICAL PROPERTIES^ 


F” Quantity 

Value 

Comments 

1 Distance from the Earth 

1.496 X 10^ km 

1 AU 

Mass 

1.991 X 10^^ kg 


Radius 

6. 960 X 10^ km 

^Siarth 

Angular Diameter at Earth 

1919 arc sec 


Luminosity 

S3 

3.86 X 10 erg/sec 


Equivalent Blackbody Temperature 

5780°K 


Gravity at Visible Surface 

2. 74 X 10^ cm7sec^ 

^ ^Earth 

Rotation Period^ 

Approximately 27 d 


General Magnetic Field 

Approximately 1 gauss 
at surface (local fields 
can be much greater) 


Magnetic Period 

Overall reversal 
approximately every 
11 years (22-year 
period) 


Tilt of Solar Equator to Ecliptic 

7 deg 



a. The information presented is based primarily on Reference 1~3. 

b. Note there is differential rotation between different latitude zones: a 

26. 9- day period at the equator and a 29.5-day period at 45 deg 
latitude (paragraph 1.4.3), 


1.4 SOLAR EMANATIONS - DESCRIPTIVE 
1,4.1 The Nature of the Sun’s Output 

The end result of the energy-producing fusion reaction in the solar core is the 
emission of energy and mass from the Sun into the solar system and beyond into inter- 
stellar space. This solar output occurs in the form of electromagnetic radiation, 
plasma, and energetic particles. Table 1-2 [1-4] provides a summary. 

The bulk of the energy export occurs in the form of electromagnetic radiation. 
Although neutrino flux is potentially the second largest mode of energy loss from the 
Sun, the theoretically predicted neutrino flux has not been confirmed experimentally. 
The discrepancy has not been resolved, since, as Bahcall and Davis [1-1] emphasize, 
experimental errors are sufficiently large to allow measured neutrino fluxes from zero 
to one-third of the predicted values. Neutrinos interact very weakly with all 
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TABLE 1-2. TOTAL ENERGY FLUX OF SUN 


-2 -1 

(erg cm sec t at 1 AU) 


Radiant Energy 

solar constant S = 1. 37 x lo” 

Neutrinos 

4 

predicted flux = 5 x 10 

Solar Wind 

1 

Magnetic Field 

0.01 


materials, and they can be ignored in spacecraft design. In terms of sensible solar 
outputs, solar radiation at visible wavelengths and the solar wind plasma are the most 
Important energetically. 

Values for the fluxes of the solar outputs will be stated at one astronomical 
unit (AU) , because the most reliable measurements have come from spacecraft near the 
Earth. Electromagnetic radiative fluxes are easily scaled for other positions in the 
solar system by using the inverse square law, but scaling plasma and particle fluxes 
is difficult because of interactions that occur in interplanetary space. These mapping 
problems will be discussed later 


1.4.2 The Solar Cycle 


The solar magnetic cycle is the basic underlying cause of all the observed fluc- 
tuations over long (years) and short (days to even minutes and seconds) time scales. 
The solar cycle Was first seen in the systematic rise and fall of the number of visible 
sunspots. The apparent cyclic nature of solar activity with an 11-year period can be 
seen in the sunspot number record (Fig. 1-3) extending bach to the beginning of the 
18th century. 

Figures l-4a and l-4b show more detailed variations in sunspot number for 
solar cycles 18, 19, 20, and 21 (now only about halfway completed). One-year run- 
ning averages of the monthly mean of the Ztirich sunspot number : 


- K (10 g d- f) 


are given, where g is the number of spot groups observed on the solar dish at one 
time, f is the number of individual spots, and K is a normalization constant which 
varies from observing site to observing site. 


The model of the Earth's upper atmosphere described in Appendix A requires 
inputs of the level of solar activity expected to prevail at the time the thermodynamic 
properties of the upper atmosphere need to be hnown. These inputs are the 162-day 
means of the 10.7-cm solar radio noise flux, F^^q and the 3 -hourly values of the 


geomagnetic, Epfa^) , index. 


Since the measured F 


10.7 


data base is relatively short 


(1947 to date), the data base was extended bach to 1749 using the following equation 
to convert recorded smoothed (R2) data to smoothed solar flux (Fj^q y) data: 
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^ = 49.4 + 0.97 + 17.6 exp (-0.035 R^) 

The recorded data for F^^ ^ and R^ from 1947 to the present were used to 

derive the equation which has a linear correlation coefficient of 0. 98. This equation 
was used to construct smoothed solar flux F^^ ^ data for the period 1749 to 1946. 

After 1946, the measured values of F^^ ^ are used to compute the F^^^ 

Figure l-4c shows the 97.7, 50,0, and 2.3 percentile values of this F^^^ ^ data 

base arranged in a monthly sequence beginning at solar cycle minimum and extending 
through a mean solar cycle period, 11,08 yr. The values on the figure are the 
extreme and median values of the data set for each month since extreme value statistics 
has shown that they are most representative of the quoted percentile values for small 
data sets, 

Hale's pioneering experiments at Mt. Wilson to measure solar magnetic fields in 
sunspots at the beginning of this century established the regular reversal of field 
polarities in sunspot pairs every 11 years. Later, in the 1950s, Babcock working at 
Mt. Wilson further established the fundamental character of the magnetic cycle through 
his discovery that the Sun's general dipole field as measured at the solar poles also 
reversed near the maximum of solar activity every 11 years. Therefore, the most 
fundamental process appears to be the 22-year periodicity in magnetic fields produced 
by dynamo action deep in the solar interior, but this in turn results in an 11-year 
periodicity in the degree of solar activity and intensities of those solar outputs orig- 
inating in active regions surrounding simspots. 



Figure l-4c. Solar cycle variations of F^^ superposition of cycles 1-21. 
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1.4.3 Variation in the Sun's Output 


Both radiation and particle flixxes vary over times much shorter than 11 years. 
The Sun rotates with a period of approximately 27 days; consequently, any emission 
from, for example, a bright active region will be visible to a stationary observer for 
13. 5 days and then be inAdsible while the emitting region traverses the far side 
of the Sun, The pure square wave rotational modiilation is smoothed by the system- 
atic darkening of the active region near the two edges of the solar disk (i.e., limb 
darkening effects), and also by the slow evolution of active regions. Active regions 
where sunspots form have lifetimes of up to ten solar rotations, and during this time 
the emission shows the effects of both slow morphological change in the local magnetic 
field structure and the much faster fluctuations (seconds to minutes) caused by the 
occurrence of solar flares. We typically see more than one active region on the disk 
at a time; therefore, the observed waveform for the solar emission seen from a single 
point in the solar system may appear to be quite complicated and even aperiodic at 
times. The observed fluctuations are further complicated by the differential rotation 
of the Sun; i.e, , equatorial regions of the Sun rotate faster than regions near the 
poles. Sunspot observations give the rotational period as 26. 9 days for the equator 
and 29.5 days for 45 deg solar latitude. Nevertheless, autocorrelation analyses of 
such signals over spans of several rotations show the 27-day recurrence effects very 
clearly. Flares occur randomly in time and have a wide range in emission strength, 
but the large flares produce very large excursions in the Sun's radiation at wave- 
lengths below 100 nm. 

In considering variation of solar electromagnetic radiation caused by solar 
activity, we usually identify the active region effects with increases in emission 
strength, but quite a different situation exists for sunspots and the solar wind. 

Dark sunspots are obviously locales where visible radiation is reduced substantially; 
as discussed in the next section, the large sunspots appear to decrease the solar 
constant slightly as they cross the visible solar hemisphere. Similarly, times of 
increased solar wind speed are associated with the so-called "coronal holes," which 
appear dark in X-ray images of the Sun and which do not contain active regions. 
These explanations of variation in solar output levels became possible only through 
measurements made from space platforms over the last two decades. The most fre- 
quent cause of variation is the rotation of different types of solar regions across the 
visible disk, already described, together with the occurrence of solar flares and 
coronal transients. 


1.5 SOLAR ELECTROMAGNETIC RADIATION 

In presenting the data for the Sun's electromagnetic radiation, it is necessary 
to discuss the spectral distribution of this energy and the total amount of energy 
under the curve. The latter integral quantity measured at 1 AU is the Solar Constant 
-2 -1 

(units of watts in sec ). Measurements of the spectral distribution yield values for 
the Sun's spectral irradiance (units of energy flux in a spectral band of specified 
width). 

The spectral energy distribution for the Sun resembles a Planck curve with an 
effective temperature of 5800°K, which means that the bulk of the solar energy lies 
between 150 hm and 10 ym with a maximum near 450 nm. This portion of the 
solar spectrum is emitted primarily by the solar photosphere, the deepest level 
observable. There are, however, striking and variable solar emissions outside this 
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wavelength band that originate in the chromosphere and corona as well as in transient 
events such as flares. When we study the solar constant then, we study primarily 
the photosphere, which is dominated by convective motions in an ionized medium with 

T = 5800°K and N = 10^^ electrons cm 


The state of the art in the measurement of the Sun's electfomagnetic spectrtun 
is given in considerable detail in White's [1-6] compendium on the solar output and 
its variation. The following discussion draws heavily on the individual review chapters 
[1-6], but new material from more recent measurements will be emphasized. 


1.5.1 Measurements of the Solar Constant 


Adopted value: 1371 + 5 W m 

Some of the most intriguing recent results in solar physics have come from 
spacecraft measurements of the solar constant with cavity radiometers capable of pre- 

• 4 

cisions of one part in 10 . This t3rpe of radiometer was developed at the Jet Propul- 
sion Laboratory during the 1970s [1-7]. These devices are simple and rugged enough 
for flight outside the Earth's atmosphere onboard rockets and orbiting spacecraft. 

To appreciate the impact of the new measurements, consider that Frolich [1-8] 

_ 2 

derived a consensus value of 1373 ± 20 W m from all solar constant measurements 
made between 1965 and 1975. Willson and Hickey [1-9] immediately challenged this 
value on i«the «g^unds that it did not take into account the much higher precision 
of the very recent measurements. Their weighted mean from all credible aircraft, 

_2 

balloon, and space experiments conducted is 1371 3: 1 W m . New measurements con- 
tinue to come from the Nimbus 7 spacecraft [1-10] and the Solar Maximum Mission 
(SMM) satellite [1-11], The Nimbus 7 experiment began on November 16, 1978, and 

- 2 

Hickey et al. [1-12] report a mean solar irradiance value of 1371 ± 5 W m from the 
Earth Radiation Budget (ERB) radiometer. Willson's experiment on SMM began opera- 
tion in February 1980, and WiRson et al. [1-11] report a total irradiance value of 
_2 

1368 ± 6 W m . Comparison of these four values for the solar constant shows the 
state of the art in the absolute determination of the total irradiance. Using 
these data it is difficult to make a case for a secular change in this solar output, but 
the prevailing trend in the measurements has been downward since 1966 [1-13]. 


1,5.2 Short-Term Fluctuations in the Solar Constant 


The Nimbus 7 and SMM experiments mentioned in the preAdous paragraph give 
the first credible empirical evidence for variation in the total solar irradiance. Both 
experiments [1-10,1-11] show maximum downward excursions of 0.2 percent in the 
solar flux level, and the larger amplitude decreases are well correlated with the 
passage of large sunspots across the solar disk. These excursions have durations of 
approximately 10 days. Although these variations are quite small, they are important 
in solar physics because we are now forced to accept the idea of a net radiative loss 
in a sunspot; i.e. , the missing sunspot flux is not made up by a slight increase in 
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the brightness of the quiet photosphere surrounding the spot. Since the sunspot 
geometries and positions are known, this effect can be modeled by using mean sunspot 
brightnesses. Preliminary attempts at simulation account for only about half of the 
variance in the solar constant data. Nevertheless, the first '-order effect appears to 
be the disk passage of dark sunspots, which is the simplest explanation possible. 


1, 5. 3 The Solar Spectral Irradiance 


Because the solar constant is so heavily weighted by the photospheric emission, 
its solar cycle variation is probably quite small (less than 0.5 percent), but the situa- 
tion is quite different for the long and short wavelength! ends of tlie Sun's spectral 
energy distribution. These extremes of the solar spectrum and strong absorption 
lines show effects of the solar activity cycle very clearly. This apparent different 
sensitivity to solar activity occurs because the emission in these bands originates in 
higher layers of the Sun's atmosphere, i.e., the chromosphere and corona, or in 
very energetic transient phenomena such as solar flares. 

Absolute measurement of the solar spectrum is a difficult photometric task; con- 
sequently, our present pictiare of the spectral energy distribution is not completely 
satisfactory. Chapter IV in White's [1-6] compendium summarizes the situation as of 
1976, but since that time new experiments have been started to measure the absolute 
flux levels from the time of solar minimum in 1976 to the maximum of solar cycle 21 in 
1979. Heath [1-14] summarizes many of the new results which are just now being 
published. Table 1-3 defines several bands in the spectrtim and lists typical flux 
levels and their variability. Figure 1-5 shows the solar spectral energy curve as 
constructed from measurements available in 1976. 


TABLE 1-3. FLUX LEVELS IN THE SOLAR SPECTRUM 


Spectral Region 

Wavelength Band 

Flux Level^ 
Across Band 

Variability in Cycle 

Radio 

A > 1 mm 

10'^ - 10'^^ 

lOOX 

Far Infrared 

1 mm > \ >10 ym 

10-2 

V 

Infrared 

10 ym > A >0. 75 ym 

10 - 10^ 

9 

Visible 

0. 75 ym > A > 0. 3 ym 

10® 

<1 percent 

Ultraviolet 

0. 3 ym > A > 0. 12 ym 

10^ - 10® 

1 percent to 2X 

Extreme UV 

0. 12 ym > A > 0. 01 ym 

10^ 

lOX 

Soft X-ray 

0. 01 ym > A > lA 

10^ - 10' 

lOOX 

Hard X-ray 
& Gamma ray 

1 A > A 

10'^ - 10'® 

10 to lOOX 


- 2 - 1 

a. In units of erg cm sec y 
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Figure 1-5. The Solar spectriam as derived from measurements 

available in 1976 {1-6]. 

The energetically dominant spectral region from 1200 A to 10 ym (the ultraviolet, 
visible, and infrared bands in Table 1-3) contains 99.5 percent of the total radiant 
energy from the Sun and shows variations comparable to current measurement errors, 
except at short wavelengths near 1000 A where the variability increases rapidly to 2X. 
This increased variability occurs because the blackbody radiation cwve for the photo- 
sphere decreases rapidly to small flux values in the UV, and, as a result, the spectrum 
switches from a continuum-dominated spectrum to a pure emission line spectrum near 
1400 1 . This emission line character continues all the way to the gamma ray region. 

The measured variability during a solar rotation and over the cycle increases because 
these emission lines originate in the corona, chromosphere, and flares rather than in 
the much brighter photosphere. 

In the long wavelength region beyond 10 ym into the radio spectrum, we have 
a weak continuous spectrum again, showing increased variability as we go further 
from the visible spectrum. Once again, this increase arises from the increasing 
dominance of coronal emission as we go into the radio region. The low density (N 

Q _ O 0 

^ 10 cm ) and a high temperature (T^ ^ 10 K) solar corona displays large changes 

in its structure over a solar cycle as well as many short-lived transient events asso- 
ciated with individual active regions on the Sun. 
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Even though our picture of the solar spectrum is still incomplete in terms of 
absolute flux levels measured over one solar cycle, the situation is improving. Since 
many of the newer experiments began near solar minimum in 1976, by 1987 we should 
have our first continuous record of the variation of the spectriim over one solar 
activity cycle. 


1.6 SOLAR PLASMA EMISSION 

The Skylab epoch from 1972 to 1974 led to a significant improvement in our 
understanding of the solar corona and solar wind because for the first time the direct 
link between the solar atmosphere and the interplanetary medium was found in the 
so-called coronal holes. Hundhausen [1-15] summarizes the phenomenological descrip- 
tion of the corona and its effect on the interplanetary medium. This early discovery 
epoch came in the declining phase of solar cycle 29; consequently, confirmation of the 
phenomenological model circa 1978 awaits measurements still to be made throughout the 
current cycle. Even though this description is not fully confirmed by observations 
over a full solar cycle yet, it succeeds in answering many long-standing questions 
about the origin of interplanetary disturbances and magnetic sector structure. 


1.6.1 Properties of the Mean Solar Wind 


The properties of the mean solar wind are summarized in Table 1-4. 


TABLE 1-4. PROPERTIES OF THE MEAN SOLAR 
WIND AT 1 AU 11-16] 

Radial Velocity = 468 i. 116 km. sec’^ 

Particle Density — 8. 6 3: 6* 6 cm’ ^ 

- fi 

Mean Field = 6. 2 ± 2. 9 x lo gauss 

Proton Temperature - 1.25 ± 0.9 x lo^ deg K 

The uncertainties are 1-sigma limits. For comparison, note the properties of 
the four high-speed streams in Table 1-5. 

TABLE 1-5. PROPERTIES OF HIGH-SPEED 
STREAMS (1973-1974) [1-15] 


Radial Velocity - 

600 - 700 km sec ^ 

Proton Flux — 

8 “2-1 
4 x 10 proton cm sec 

Energy Flux = 

« -2 -1 
2 erg cm sec 

Field - 

^5 

4 X 10 gauss 
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1.6.2 The Solar Wind and the Interplanetary Magnetic Field 

The mean solar wind is a flow of solar plasma past the Earth at a speed of 
approximately 450 km /sec. This flow was predicted by Parker [1-17] as the natural 
consequence of the outward pressure gradient in the high- temperature solar corona. 

It is not everywhere a simple radial flow because the Sun rotates and has magnetic 
fields extending into the corona. The Sun's weak dipole field of 1 to 10 gauss field 
strength emerges near its poles and extends out into space where it acts as a cage 
for the solar wind flow, but the solar wind pressure is sufficient to drag the field 
lines out radially at the solar equator and create an equatorial neutral sheet separating 
the regions of N and S magnetic polarity. Figure 1-6 shows this distortion of a 
dipole field as emanated from a simple axially symmetric model. These interplanetary 
fields of solar origin are stretched along spiral tracks because the solar wind carries 
the field lines outward while the Sun continues to rotate with its 27-day period. The 
ability of the solar wind to produce neutral sheets and flows along spiral tracks that 
extend outwards into the solar system is the reason for much of the structure observed 
in the solar wind plasma at 1 AU. 

2 


I 


0 


I 2 3 4 

Figure 1-6. Distortion of a simple dipole field by a steady radial solar wind. 
Solid lines show field configxiration resulting from flow, and dashed lines 
show the closed dipole field without flow 11-15]. 


0 
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The strong magnetic fields of solar active regions complicate this structure 
through large-scale topological changes in the solar field structure as it extends into 
the interplanetary mediiom. In these closed field regions, strong fields (>500 gauss) 
prevent flow of active-region plasma into the solar wind and thereby create zones 
in the solar wind where the plasma velocity is smaller. More important for inter- 
planetary conditions, however, are the unipolar regions lying between the active 
regions and at the solar poles. The magnetic field lines from these quiet regions are 
open; i.e., they extend far out into space and provide an easy path for the solar 
wind flow (Figs. 1-6, 1-7, and 1-8). Such regions on the Sun are identified as the 
so-called coronal holes that appear to be the sources of high-speed solar wind streams 
with plasma velocities of up to 800 km7§ec or more. 

The properties of the interplanetary magnetic field are principally determined by 
separation of the N and S polarity regions, the equatorial neutral slieet, and tlie 
localized high-speed flow of plasma out of coronal holes. Measurements of the inter- 
planetary field in the ecliptic show the abrupt polarity changes that define the inter- 
planetary sector boundaries. These sectors of one dominant magnetic polarity arise 
because of the deviation of the neutral sheet from the plane of the ecliptic. In the 
idealized case of a perfect magnetic dipole, the neutral sheet woTild be planar and 
tilted relative to the ecliptic plane because of the 7 deg inclination of the Sun's rota- 
tion axis (Fig. 1-7). As the Sun rotates, this idealized neutral sheet would appear 
to wobble, and the field polarities observed at the Earth would cTiange sign each time 
the intersection of the neutral sheet and the plane of the ecliptic sweeps across the 
observer. For the simple model shown in Figure 1-7, the observer would see two 
sector boundaries separated by 13.5 days in the interplanetary field. However, 
tightly closed magnetic fields from active regions warp the effective neutral sheet and 
produce additional intersections of the sheet with the plane of the ecliptic, which in 
turn leads to more frequent magnetic sector boundary crossings at the Earth (Fig. 

1-8). An instrument flying well out of the ecliptic at heliographic latitudes greater 
than about 30 deg would probably never see a sector boundary l it would always see 
one magnetic polarity unless the sense of the global solar field changes. The mapping 
of the solar wind flow and interplanetary fields to high solar latitudes remains as an 
unfinished task in the exploration of the solar system. 

Interplanetary magnetic field lines and, hence, sector boundaries, follow an 
Archimedes spiral pattern because the field lines are rooted in the solar surface and 
are carried outward by the radially flowing solar wind as the Sun rotates. To an 
observer located above the ecliptic plane, the pattern _would resemble a slowly spinning 
pinwheel. This is indicated in Figure 1-9, which shows sector boundaries inferred 
from magnetic polarity data taken by spacecraft at 1 AU. 

To gain some insight into the spiral structure, consider a stream with a 500 
km /sec radial flow. A parcel of plasma in the stream takes 18 days to travel to 
Jupiter's orbit, and its trajectory will be a spiral of 0.7 of one turn. At Saturn's 
distance, the trajectory is wound 1.25 turns and increases to 5 turns at the orbit of 
of Pluto. Consequently, the spiral pattern is not tightly wound inside the orbit of 
Saturn. 


1.6.3 High-Speed Streams 

In addition to the large-scale effects of the neutral sheet, solar wind structure 
is strongly influenced by high-speed solar plasma streaming out of coronal holes. The 
solar surface field topology is such that a region with an open magnetic field can form 
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Figure 1-7. A schematic drawing of a simple solar wind and field 
configuration defining the Sun's equatorial neutral sheet and its 
tilt relative to the plane of the ecliptic [1-15]. 
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Figure 1-8. Schematic drawing of a realistic structure for the solar corona 
and neutral sheet, typical of the declining phase of the solar cycle. 
Note that the neutral sheet as stown crosses the ecliptic 
more than twice [1-15]. 
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Figure 1-9. Observed polarity of the interplanetary magnetic field at 1 AU 
and the spiral-shaped magnetic sector boundaries inferred from 

field observations [1-181. 





and allow hot coronal plasma to stream directly into interplanetary space. The stream- 
ing lowers the corona density to produce the so-called "hole" observed in coronal 
X-ray emission maps. These streams also have an Archimedes spiral shape but are 
less tightly wound than the ambient solar wind because of greater flow speed. 

The overtaking of slower moving interplanetary material by a high-speed stream 
eventually produces a shock front along the leading edge of the stream (Fig. 1-10). 

At those times during the solar cycle when coronal holes and high-speed wind streams 
are present, the shocks and spiral tracks of the long-lived streams play an important 
role in particle acceleration and modiilation of the flux of energetic particles in the 
outer regions of the heliosphere. 



Figure 1-10. Schematic drawing showing how shock fronts 

are produced 11-19], 


1-20 



Hundhausen [1-15] summarized the properties of four high-speed streams in 
Table 1-5. Although the proton density in a high-speed stream is lower than for the 
mean solar wind, the velocity is higher and the particle flux remains about the same 
8 _ 2 - 1 

(a few X 10 protons cm sec ) . The magnetic field strength in streams is much 
higher than in the mean solar wind. 

Figure 1-11 shows the observed variation of flow speed, density, and tempera- 
ture with time as a high-speed stream passed by the Vela 3 spacecraft in 1966. 

Notice the shock front and rarefaction region trailing behind it. 



a. iQ^ I i ^ ^ 1 

APR. 2! 22 23 24 25 26 27 

TIME OF OBSERVATION 


Figure 1-11. Seven-day time profiles of velocity, density, and temperature for a 
well-defined high-speed solar wind stream. Observations (dots) made by VELA- 3 
in April 1966, (The curves are from a theoretical model.) [1-19] 
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Hundhausen [1-15] suggests that the high-speed streams are at the high-speed 

limit, while the flow near a neutral sheet is the low-speed limit (300 km sec ^). Since 
the solar poles have permanent coronal holes for most of the solar cycle , they may be 
the sites of long-lived, high-speed streams perpendicular to the ecliptic; the existence 
of such solar streams is only conjecture since all observations to date have been con- 
fined to ±30 deg of the ecliptic plane. 


1 1.6.4 Coronal Transients 

Solar material, at times, is suddenly ejected into space with velocities on the 
order of 1000 to 1500 kmVsec, producing shock fronts as the fast material snow plows 
through the ambient wind, just as occurs for high-speed streams. In many cases 
I such transients take place in the solar corona without the occurrence of a flare in 

I the lower solar atmosphere. Many appear to be associated with the eruption of large, 

I normally stable solar prominences. So, added to the picture of a wobbling neutral 
I sheet and spiraling streams in the heliosphere are parcels of coronal plasma traveling 

I at speeds of about 2X the solar wind which occur randomly in time. Models of the 

I heliosphere with such parcels distributed along spiral tracks have been studied in the 

I cosmic-ray modulation problem, but their role in the acceleration and diffusion of 
I energetic particles is not yet understood. These transients are called coronal trans- 

I ients; they appear to be a true coronal phenomenon, possibly arising from large-scale 

I topological rearrangement of the magnetic fields that extend into the outer solar 
I atmosphere. 


1.6.5 Spatial Variation of Solar Wind Properties 

Our conception of the solar wind structure is strongly biased by the predomi- 
nances of measurements in the plane of the ecliptic and near the Earth. However, 
estimates of the solar wind density at heliographic latitude up to 26 deg have been 
obtained from data taken on the Viking spacecraft during the last solar minimum in 
1976 (Muhleman and Anderson [1-20]). They find the density to decrease according 

-2 7 - 2 

to a linear combination of two power laws with r ‘ and r radial dependences, but 

the r " term disappears at latitudes above 20 deg, leaving only the r dependence. 
This behavior is expected from the simple models shown in Figures 1-6 and l-<7 where 
the plasma is confined in the equatorial region by the dipole field. At highier latitudes 

the flow straightens the field lines to allow the density to follow an r dependence 
expected from number density conservation. 

The mean solar wind speed increases with solar latitude, and the current esti- 
mates of the latitudinal gradient range from 2 km sec ^ deg ^ to 15 km sec ^ deg 
Stable plasma streams are believed to flow from the solar poles with velocities on the 

order of 700 km sec 


1.6.6 Variation of the Solar Wind During the Solar Cycle 

In situ measurements of the properties of the solar, wind show variability at all 
time scales. Changes in averages over time periods of the order of the 11-year solar 
activity cycle have been attributed, at least in paT*+ . tr> the relatively greater number 
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of high-speed streams that appear during the declining phase of the cycle. The mag- 
nitude of the interplanetary magnetic field remains relatively constant, though its 
components are extremely variable. The helium abundance appears to vary in phase 
with the solar cycle. 

The spatial and temporal structure of the solar wind also varies during the solar 
cycle. The pattern of change is determined by the cyclic variation of the surface 
magnetic fields on the Sun. Hundhausen’s picture of the variation begins with the 
reversal in the solar magnetic field polarity that occurs approximately 2 years after 
the maximum of solar activity. At this time the Sun’s polar fields stren^hen, and 
the two polar coronal holes reappear as stable structures lasting for the next 6 years 
or so. As the number of active regions declines, the solar magnetic field pattern at 
low latitudes evolves into larger and larger structures of the same magnetic polarity. 
This in turn leads to the appearance of small, short-lived coronal holes at latitudes 
below 45 deg. The solar wind now begins to show short-lived streams with moderate 
plasma flow speeds. This magnetic evolution to larger scale structures continues 
through the declining phase of the cycle until the low -latitude coronal holes connect 
with the polar holes and appear as low-latitude extensions of the magnetic polar caps. 
During this time the solar wind is dominated by high-speed solar wind streams lasting 
for several solar rotations. Empirically, these effects appear as the 27-day recurrence 
patterns commonly seen in geomagnetic activity and solar wind speed and magnetic 
sector structure. 

Figure 1-12 shows the daily averaged polarity of the interplanetary magnetic 
field at 1 AU as inferred from ground-based observations of the geomagnetic field. 

The data are arranged into rows of 27 days so that recurrent sectors associated with 
long-lived features on the Sun appear as columns of predominantly one polarity. 

Notice the bimodal nature of the display for the 2 or 3 years before the minimum 
of solar cycle 20 (early 1976). This has been clearly identified with two large and 
long-lived coronal holes which extended to the equator from the north and south polar 
holes at roughly opposite sides of the Sun. These holes and their associated high- 
speed wind streams are shown in Figure 1-13. Notice that the magnetic field polari- 
ties of the streams (indicated by + and - symbols and taken from the data of Fig. 

1-12) agree with the polarities of the holes as inferred from measurements of the 
underlying photospheric fields (indicated similarly on the coronal maps). 

As the active regions associated with the next solar cycle appear at helio graphic 
latitudes of about 45 deg, the polar cap fields begin to weaken, and the scale of low- 
latitude coranal holes decreases. During the rising part of the 11-year cycle, high- 
speed streams occur less frequently. The sixrface field structure on the Sun becomes 
more chaotic, and the frequency of transients such as solar flares increases. This 
evolution toward a more chaotic condition continues until solar activity peaks, and 
the process begins to repeat itself. 


1.7 HIGH-ENERGY PARTICLES IN THE HELIOSPHERE 

At the beginning of the 20th century, experiments on the conductivity of air 
led to the discovery of a flux of high-energy particles at the Earth. These early 
studies established that the most energetic particles arrive at the Earth from all 
directions and, for this reason, they were called galactic cosmic rays. During the 
1940s, however, lower energy particles observed at high altitudes were definitely 
identified as solar in origin. They are frequently associated with large solar flares. 
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CARRINGTON LONGITUDE 


Figure 1-13. Top: Solar wind speed observed by spacecraft at 1 AU. + and - 

indicate magnetic polarity [1-15]. Bottom: Coronal brightness contour 

map showing the north and south polar coronal holes in 1974. 
(Carrington longitude refers to a coordinate system fixed on the Sun which 
rotates rigidly with a period of 27.14 days as seen from the Earth.) 
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The mean energy of the galactic component is on the order of 10 BeV, with energies 

9 

as high as 10 BeV. Solar particle energies rarely exceed 1 BeV, and they have a 
mean energy on the order of 100 MeV. 

All of the elements from hydrogen up to iron in the periodic table have been 
identified in cosmic rays, but their abundances depart substantially from the standard 
cosmic abundance; cosmic rays are abnormally rich in carbon and the heavier elements 
such as iron. This abnormal composition is attributed to the ultimate sources of the 
particles in the galaxy. Rgure 1-14 shows the .galactic cosmic-ray energy spectra 
for various species. Cosmic rays remain objects of research in cosmology, solar 
physics, high-energy physics, magnetohydrodynamics, and stellar structure. 



Figure 1-14. Integral energy intensities for galactic 
cosmic ray nuclei [1-21]. 

The most fascinating mystery about galactic cosmic rays is their origin. Local 
source conditions must be extreme to accelerate particles to BeV energies, and 
supernovae are prime candidates. 


1.7.1 High-Energy Particles from the Sun 


There are two classes of high-energy particles of solar origin, both closely 
associated with very large solar flares and eruptions often occurring in the declining 
phase of the 11-year solar cycle. The frequency of their occurrence tends to follow 
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the degree of solar activity. The two components are distinguished by their arrival 
time at the Earth and, hence, their energies. The highest energy particles arrive at 
the Earth approximately 20 min after the occurrence of a large flare; they appear to 
be accelerated particles directly emitted from the flare volume. The second component 
is observed approximately 24 hr later in association with a magnetized plasma cloud 
passing by the Earth. These late arrivals appear to be trapped within the plasma 
cloud and carried to the Earth at speeds of about 1000 km7sec. They enter the 
Earth’s atmosphere isotropically as the cloud envelopes the Earth, whereas the higher 
energy component is strongly directional, coming predominantly from the western 
hemisphere of the Sun. This east to west asymmetry for location of source regions 
on the Sun shows that the particles have spiral trajectories similar to the pattern of 
interplanetary magnetic field frozen into the radial solar wind flow, and this immediately 

establishes the ability of relatively weak magnetic fields (3 x 10 gauss) to influence 
the particle diffusion over long distances in the heliosphere. 

Because of their close association with large flares, the energetic particle bursts 
are relatively infrequent. For example, only 34 events occurred with particle energies 
greater than 30 MeV in solar cycle 19 (1954 to 1965), while for cycle 20 (1965 to 1976) . 
there were fewer than 20 such events. For this reason, such high-energy particle 
events are unimportant in total energy relative to the much slower but steady solar 

/ 2 ** 1 ” 1 

wind. Typical proton fluxes at 1 AU are 10 to 100 protons cm ster sec with 

-2 -1 -1 

energies above 10 MeV, but this falls to below 10 protons cm ster sec for par- 
ticles above 30 MeV, illustrating the steepness of the energy spectra of particle 
bursts associated with flares 11-21]. The total number of protons associated with 

7 8 "2 

single events is 10 - 10 protons cm on the average, but the particularly intense 

10 -2 

bursts deposit 10 protons cm at energies greater than 10 MeV. 

The higher energy particles are easier to detect, and most of our experience 
to date is for energies above 10 MeV. Improved detectors for particles with energies 
less than 10 MeV show the existence of a fluctuating background. The Sun is 
apparently a continuous emitter of these lower energy particles. Tkey appear to come 
from active regions, but the continuous fluctuation in flux level prevents association 
with particular events such as sub-flares. 

Figure 1-15 shows energy spectra for protons and alpha particles for seven 
major events between 1967 and 1969. 

Lanzerotti [1-21] discusses the solar cycle cosmic ray variation; his figure for 
variability of the solar component is shown as Figure 1-16. The solar component 
follows the activity cycle, with the largest number of events occurring near solar 
maximum when the flare production rate is the highest. For reasons not yet under- 
stood, the most intense particle events tend to occur approximately 2 years after the 
maximum . 


1.7.2 Influence of Solar Activity on Energetic Particle Fluxes 


Even though galactic cosmic rays themselves are not related to the Sun, their 
flux at the Earth is modulated by solar plasma effects. These very high-energy par- 
ticles arrive at the Earth uniformly from all directions. Records made at the time of 
great solar flares show a decrease in the observed galactic flux level beginning 24 to 
48 hr after the flare. These decreases, known as Forbush decreases, are on the order 
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1-15. Event -integrated proton and a particle spectra for 
seven major events between 1967 and 1969 Tl-21]. 













Figure 1-16. Variability of the solar component [1-4]. 


of 10 percent and may last for several days before the galactic intensity recovers its 
preflare level. This effect coincides with the passage of a cloud of solar plasma 
which disturbs the local conditions for diffusion and' propagation of high energy 
particles . 

Figure 1-17 shows galactic cosmic rays plotted together with the Zurich snnspot 
number [1-22]. The galactic component is strongest at solar minimum and weakest at 
solar maximum. This modulation is almost 100 percent for energies below 500 MeV, 
according to Lanzerotti [1-21], but drops to a few percent at energies of 100 BeV. 

It is believed to arise from the interaction of the high-energy particles with high- 
speed plasma clouds and shocks produced by transients on the Sun. The modulation 
effect is the subject of intense study and controversy at the present time. 
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SECTION 2. lERRESTRIAL SPACE (THE EARW) 


2,1 INTRODUCTION 

The region of space discussed in this chapter extends from approximately 80- 
km altitude above Earth's surface, the lower boundary of the ionosphere, to the 
boundary of the Earth's magnetosphere, where the unperturbed interplanetary medium 
is reached. Thus it extends to the bow shock approximately 100,000-km radial dis- 
tance (16 Earth radii, R^) in the sunward direction, and down the magnetospheric 

tail in the anti-sunward direction, to a larger distance whose actual value depends on 
the effect of interest. This distance may extend beyond the lunar orbit at 60 R^ 

(approximately 385,000 km). "Terrestrial Space" includes the orbits of lower artificial 
Earth satellites to beyond the geostationary orbit (approximately 7 Rg). 

A spacecraft operating in this volume is affected by the following: 


Effects 


Acceleration, torque 
Heating and drag, torque 
Heating and power, drag, torque 

Torque, surface charges, induced electrical potential 
Electric potential 

Charging, change of refractive index for em waves 
Ionization, radiation damage 
Mechanical damage 

It is not possible in a single chapter to explain all the interactions which 
influence the Earth's atmosphere, ionosphere, and magnetosphere. Indeed, the I 
detailed quantitative understanding of Terrestrial Space in terms of fxmdamental 
physics is still on the frontier of research. However, much useful information has 
been gained from ground-based and space measurements and partial theoretical under- 
standing has been achieved. 

The following paragraphs briefly describe the physical relationships between 
different parts of Terrestrial Space. Succeeding paragraphs describe the most 
important factors quantitatively using empirical or parametric models. For an excel- 
lent general reference review the AFCRL Handbook of Geophysics [2-1]. 


Environmental factors 
Earth's gravity 
Neutral atmosphere 
Sunlight and albedo light 
Magnetic field 
Electric field 
Thermal plasma 
Fast charged particles 
Meteoroids 


2.2 GENERAL DISCUSSION 

The gravitational field results from the mass of the solid Earth and reflects 
the distribution of that mass. It traps the neutral atmosphere. 
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Above approximately 90-km altitude the ion-electron pairs, formed in sunlight 
by photoionization and, to a lesser extent, colHsional ionization by cosmic rays and 
am-oral particles , recombine so slowly that in quasi-equilibrium an appreciable plasma 
density forms. While the plasma density remains less than the neutral density below 

3 

an altitude of 10 km, electrical forces dominate the behavior of the plasma above 
approximately 90 km where the electron cyclotron and plasma frequencies (1 to 3 MHz) 
exceed the collision frequency of electrons with neutrals. The ionospheric plasma is 
most dense at 200 to 400 km and supplies plasma to higher altitudes by bulk flow. 

At the 100 to 200-km level the ionosphere has some electrical conductivity perpendicular 
to the magnetic field, B, because of collisions with neutrals. (At higher and lower 
altitudes this conductivity is very small.) Parallel conductivity is larger than the 
perpendicular at all altitudes. 

The magnetic field has two sources: currents inside the Earth and currents 

in the magnetosphere. At the surface approximately 99 percent of the field is due to 
internal currents. Without the solar wind flowing past the Earth, this intrinsic 
geomagnetic field would decrease with distance with the dominant dipole term giving a 

3 

I7r dependence. It would have approximate cylindrical symmetry about the dipole 
axis, inclined approximately 11 deg to the Earth’s rotation axis. In fact, the solar 
wind plasma cannot penetrate the geomagnetic field freely. It compresses the field to 
about 10 Rg distance at the subsolar point and flows around the volume in which the 

field is thus confined. The magnetic field lines from the polar regions of Earth, 
which would have reached an equatorial distance greater than 10 R^, are dragged 

along by the solar wind to form a tail. Figure 2-1 shows the resiilting shape of the 
magnetosphere. Another way to describe this phenomenon is to note that the solar 
wind plasma in flowing around the magnetosphere produces electric currents which, 
when added to the Earth's internal currents, produce the field of Figure 2-1. Figure 
2-2 shows the bow shock, approximately 4 R^ sunward of the magnetopause. The 

supersonic solar wind, flaring back around it, is slowed down and partially 
thermalized. 

Although most impinging solar wind flows around the magnetosphere, some 
flows in at the cusps, as magnetic flux is peeled off the sunward side and laid back 
along the tail. More plasma, energy, and momentum diffuse in along the flanks of 
the tail. (Details of these mechanisms are not established.) Thus there is general 
anti-sunward flow of plasma and magnetic flux over the polar caps and along the tail 
in the magnetopause boundary region. The plasma and magnetic flux diffuse inward, 
coupled with migration of the field lines tied to the rotating Earth. Moving within 
the plasma sheet, they convect sunward (toward Earth) to approximately 4 R^, where 

positive and negative particles then move in opposite directions around the Earth and 
out through the sunward magnetopause where they join the westward flowing ring 
current . 

Inside the plasmapause', at approximately 4 in the equatorial plane on the 

antisunward side, the magnetic field is dipole-like and co-rotates with the Earth. The 
particle density in the plasmasphere , which is filled with ions and electrons drawn up 
from the ionosphere, is relatively greater than in the convection region outside. 

At any point where the plasma convection velocity is y and the magnetic field 
B there is an electric field E = y xB normal to B. Because the conductivity is very 
high along B and very low perpendicvdar to it (with certain local exceptions) , the 
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Figure 2 - 1 . Schematic view of the magnetosphere (courtesy of W. Heikkila) [2-2] 




Figure 2-2. Plasma convection through the magnetosphere. (Question marks 
designate the assumed connection of geomagnetic and interplanetary field 
lines. Hollow arrows represent flow directions. The view is of a 
north-south meridional plane, with the Sun to the left; the dawn 
side of the Earth is in back, dusk side is above the page 12-3].) 

electrical potential difference between any two magnetic field lines is constant so that 
electric fields generated by convective flow in one region are mapped along B to other 
regions, including the ionosphere. The ionospheric conductivity across B is~relatively 
high so that currents flow there and dissipate energy by Joiile heating. The iono- 
spheric currents are supplied from the generating region of convection by field-aligned 
(Birkeland) currents. 

Another Important manifestation of the extraction of energy from the solar wind 
is the acceleration of ions and electrons to high, even relativistic, energies. Some 
portion of these charged particles precipitate into the atmosphere to form the aurorae, 
which are centered at approximately 65 deg magnetic latitude around each polar cap. 
The rest escape the magnetosphere, and some are trapped for a time as the Van Allen 
radiation . 

The processes described above fluctuate substantially. Fluctuations are pro- 
duced by variations in solar wind speed and density and changes in the imbedded 
interplanetary magnetic field (IMF). Enhanced solar wind'flow, compressing the 
magnetosphere, may lead to a world-wide magnetic storm, which gives rise to changes 
in the magnetic field at Earth’s surface and in the size of the various magnetospheric 
regions. Another source of fluctuations is sudden release of magnetic energy stored 
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in the tail producing enhanced polar ionospheric currents and auroral particle pre- 
cipitation, This latter set of phenomena is called the polar magnetic substorm. The 
role of interplanetary medium fluctuations in triggering this particular mode df release 
of the stored energy is the subject of much controversy. That external effects, such 
as a southward turning IMF, can trigger the release is clear. However, it is likely 
that the energy buildup and rapid release is an intrinsic property of the magneto- 
spheric system and woiald proceed even if external conditions were unchanging. 

To stimmarize, the solar wind and IMF streaming past the magnetic field of the 

13 

the Earth produce the structure shown in Figures 2-1 and 2-2. About 1.4 x 10 W 
from the solar wind impinges on the front of the magnetosphere of which 2 to 3 per- 
cent is absorbed [2-4]. Of this 1-3 x 10^^ W reaches the ionosphere, 85 percent 
dissipates and Joule heating and 15 percent as aurorae. Incidentally, the radiant 

1 6 

energy input of sunlight onto the top of the atmosphere is 1.75 x lO"*^ W. 

In the succeeding pages we give ^rief descriptions of each environmental factor 
and some numerical data. Reference 2-5 gives a general discussion of magnetos pheric 
modeling. 


2.3 THE SOLID EARTH AND ITS GRAVITY 

The Earth as a planet in the solar system is described by the constants in 
Table 2-1. 

TABLE 2-1. EARTH CONSTANTS (Epoch 1960.0) [2-6, 2-7, 2-8, 2-9] 


Distance from sun (average) 

1.496 X 10^ km (1.0 AU by definition) 

Eccentricity of orbit 

0.0167295 

Orbital period (sidereal) 

365.25636 days 

Radius (equatorial) 

6378.140 km 

Mass 

5.977 X 10^"^ g 

Flattening (dynamic) 

1:298.256 

Average density 

5.517 g/cm^ 

Rotation rate 0 

0.72921 X 10~‘* see~^ 

Gravitational constant, GM 

3. 986012 X 10^ km^7sec^ 

Inclination of equator 

23.45 deg 

Period of rotation (sidereal) 

23. 934 hr 


2.3.1 The Earth's Gravitational Field 

The gravitational acceleration g can be written as the gradient of a scalar 
potential U, 
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g = vu . 

It is convenient to write U in a spherical coordinate system as 

00 n 

1 - Z O') + E E 

n=2 n=2 m=l 

(Cnijjcos mA + S^^^sin mA) 

where 

r = radius from center of the Earth, km 
(^' = geographic latitude 
A = geographic longitude 
= Legendre function 
P = associated Legendre function 
GM = gravitational constant of Earth 
a^ = equatorial radius of the Earth, km 

Values for the zonal harmonics (J^) are as follows: 

Jg = (1082.7 + 0.1) X 10^® 

Jg = (-2.56 ±0.1) X 10“® 

= (-1.58 ± 0.2) X 10^® 

Jg ^ (-0.15 ± 0.2) X 10’® 

Jg - (0.59 ± 0.2) X 10’® 

Jy = (-0.44 ± 0.2) X 10’® 



(sin (j)*) 
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nm’ nm 

follows : 
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0.004 X 10"® 
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These coefficients have been determined from surface gravity measurements and 
satellite orbit observations [2-7,2--8]. The zonal harmonics (J2S J33 J4 •••) have a 

greater effect on the orbit of a satellite than the tesseral liarmonics. The tesseical 
harmonics cause oscillatory disturbances that change sign rapidly, but the effect of 
zonal harmonics is cumulative. The even coefficients, J2, ...» can be determined 

from the regression of an artificial Earth satellite orbit's node and the rotation of its 
perigee. Reliable gravitational coefficients have been obtained from low- altitude 
artificial satellites. 

A significant geodetic result thus obtained is the reliable determination of Jg 

and, hence, the flattening f. Presently accepted values of 17f range from 298.2 to 
298. 3. This flattening of the Earth causes the largest deviation of the gravitational 
field of the Earth from that of a homogeneous sphere. 

Niimerous higher-order spherical harmonic expansions have been derived from 
analyses of artificial satellite motions, but because they yield divergent results, only 
the first few terms are tabulated here. They make the largest contribution, with 
additional terms giving diminishing returns. 


2.3.2 An Approximate Geopotential Field 

The gravitational field can be well approximated using the first three zonal 
harmonic coefficients (J2, Jg, and J^) and the main tesseral 'harmonic coefficients 

(C22 and 822)* The following expression for the gravitational potential function of 

the Earth is recommended for space vehicle design studies. 
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1 - (a /r)^ (3 sin^ - 1) - (a /r)^ (5 sin^ cf)' - 3 sin cj)') 

- (a^/r)^ (35 sin"^ cj)' - 30 sin^ (j)’ + 3) 

+ 3(C22 cos 2X + S 22 sin 2A) (a^/r)^ cos^ (j)') 

2.4 NEUTRAL ATMOSPHERE AND WINDS 
2.4.1 Neutral Atmosphere Structure 



Atmospheric conditions affect Earth satellites at 90 to 1000 hm altitude by pro- 
ducing drag and, at the lower altitudes, frictional heating. Mass density and winds 
are the important factors and only neutral density need be considered. Near the 
lower limit of the range , where atmospheric density is greatest , an unpowered space- 
craft will de-orbit in a very short time; near the upper limit, the drag effect on 
orbital lifetime is almost negligible. 

Most of the density values for the atmosphere between 90 and 2500 km have 
been derived from analyses of changes in the periods of orbiting satellites. Figures 
2-3 and 2-4 show mass and number density of the 1976 Standard Atmosphere; Figure 
2-5 shows temperature [2-10]. It should be noted that as tke collision rate decreases 
with height a point is reached where atmospheric constituents are not completely mixed; 
each tends to assume its. own scale height, kT7mg, with m the molecular weight. 

When the mean free path for collisions exceeds the scale heigltt, individual molecules 
must be regarded as in orbit aix)und the Earth, though subject to collisions [2-11], 


2. 4. 1.1 Variations 

Density and other parameters vary with solar activity as the absorbed solar 
energy changes. Note that solar radiant flux at extreme liltra- violet (EUV) and other 
short wavelengths (absorbed above 100 km) varies much more than solar flux in the 
visible. 

Geomagnetic activity (indirectly related to solar activity) also affects tlte neutral 
atmosphere through the Jovile heating gaused by ionospheric currents and by precipi- 
tation of particles whose energy is absorbed in the 90- to 200-km range. 


2.4.2 Neutral Atmosphere: Chemical Composition 

In the Earth's homosphere, which extends from the surface to an altitude of 
90 to 100 km, the atmospheric gases mix thoroughly. However, above 90 km, and 
particularly near 105 km, EUV solar radiation causes molecular oxygen to dissociate. 
The atomic oxygen thus produced, together with diffusive transport processes, results 
in a changing atmosphere. Accordingly, the atmospheric chemical composition above 
the homopause at 90 to 100 km is a function of the variable amounts of EUV radiation 
received from the Sun. 


2-8 




Figure 2-3. Atmospheric density versus geometric height for the 1976 U.S 
Standard Atmosphere and for four earlier U.S. Standard and model 
atmospheres [2-10]. (For references, see original article.) 
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Figure 2-4. Atmospheric density [2-10]. 
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Kinetic Temperature (K) 


Figure 2-5. Kinetic temperature versus altitude [2-10]. 



The temperature lapse rate is influenced by solar radiation. In the lower 
thermosphere (100 to 300 km), solar radiation in the EUV band (40 to 2000 A) is 
absorbed, causing the temperature to increase steadily with altitude. Above 300 km, 
where little or no solar radiation is absorbed, tlxe temperature becomes isothermal as 
shown in Figiu*e 2-5. The isothermal temperature, which is designated as the "exo- 
spheric temperature," varies diurnally, seasonally, and with solar and geomagnetic 
activity, from about 650 to 2100°K. 


2.4.4 Neutral Atmosphere: Density Variations 

Variation in density above the homopause is closely related to the amount of EUV 
received from the Sun. Since EUV cannot be measured at the Eartli's surface, early 
investigators made use of the approximate correlation between solar activity and 
radio radiation at 10.7 cm, which is received at the surface. Tke mean daily solar 
flux at 10.7 cm measured by the National Research Council, Ottawa, Canada, has 
been adopted as an indicator of the amount of EUV radiation reaching the top of the 
Earth's atmosphere. 
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The principal periodic variability found in the solar flux at 10.7 cm occurs with 
a cycle of 27 days, corresponding to the 27-day solar rotation. The same periodicity 
is reflected in atmospheric density. Semiannual, diurnal, and solar cycle variations 
have also been identified in the upper atmospheric temperature. 

Upper atmosphere density variations also can be correlated . with fluctuations in 
the magnetic field at the Earth's surface, since geomagnetic variations are strongly 
correlated with precipitating electrons and protons in the polar regions, and pre- 
cipitating particles heat the lower thermosphere in the polar regions with resulting 
changes in the global thermospheric temperature and density. This correlation between 
changes in the geomagnetic field and atmospheric density variations is usefxil in pre- 
dicting density variations. Figures 2-6 and 2-7 show typical density profiles under 
various conditions. 



Figure 2-6. Typical daytime maximum and nighttime minimum atmospheric 
density profiles for high and low solar activity 
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Figure 2-7. Density profiles for seven values of exospheric temperature [2-14] 


Atmospheric density also varies between the winter and summer hemispheres in 
the 60 to 80 deg latitude range. Some investigators attribute the observed higher 
upper atmosphere density in the winter hemisphere to an increased concentration of 
heliiim . 

For engineering applications an upper atmosphere computer model is needed to 
predict the neutral gas properties of the atmosphere between 90- and 2500-km altitude 
for orbital dynamics, lifetime, and control and guidance analyses. Only nominal 
atmospheric density data is needed below approximately 90 km. The Global Reference 
Atmosphere Model is used with the MSFC Lifetime Prediction Computer Routine at MSFC 
for calciilating orbital lifetime. This model is a combination of the 4-D Model Atmos- 
phere, the Groves Model, and the Jacchia 1970 Atmosphere Model [2-12], Solar 
activity, either observed or predicted, is an input. Reference 2-13 provides for one 
of the many techniques for predicting future levels of solar activity for inputs to the 
atmospheric model. 


2.4.5 Neutral Atmosphere: 2,500 to 65,000-km Altitude 

Gas pressure decreases exponentially with increasing altitude above 2,500 km 

_ 1 c 2 1 

until it reaches the interplanetary value of 10 N/cm (10 mb) near 20,000 km. 

It then remains relatively constant with increasing altitude. 

Gas density decreases exponentially above 2,500 km to the interplanetary value 

_ oo o 

of approximately 10 gjcm near 20,000 km and then remains relatively constant 
with altitude. 
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Kinetic temperature increases above 2,500 km to the interplanetary value of 

5 

about 2 X 10 °K at 20,000 km and then remains relatively constant. 

Additional information on the structure and variability of the very high thermo- 
sphere is given in References 2-15 and 2-16. 


2. 4. 6 Winds in the Neutral Atmosphere 


Information about upper atmosphere winds is derived from satellite drag data 
and ground-based optical and incoherent-scatter-radar data. These are fit into a 
theory based on the following considerations: the upper atmosphere is heated by 

absorbed solar ultraviolet energy and, in the auroral zone, by particle fluxes plus 
the Joule heating of the auroral electrojet. Some heating also results from waves in 
the lower atmosphere [2-17]. Heated air rises and flows towards cooler regions where 
it descends, just as in the troposphere. Above 300- km viscosity tends to make the 
motion uniform, i.e. , slablike. Collisions with ions tend to reduce the velocity 
difference between ions and neutrals; the effect is strongest in the sunlight ionosphere 
where the ion density is highest. The ions are tied to the geomagnetic field or drift 
where there are electric fields, as in the polar regions. 

Based on these ideas, Geisler [2-18] and Kohl and King [2-19] have computed 
wind fields from Jacchia's [2-20] static diffusion model, J64, of the thermosphere, and 
Geisler [2-21] has computed winds using Jacchia and Slowey's [2-22] modifications to 
J64. In these studies, the temperature, density, and mean molecTilar weight fields 
are used to derive pressure gradients. These are combined with ion drag forces to 
generate wind fields. Figures 2-8 and 2-9 depict the distribution of winds as a func- 
tion of time and latitude, respectively, assuming only solar heating. The predicted 
meridional winds become large at night, approaching 225 m/sec above 300 km at 
0200 LT , and then decrease to about 100 m/sec above 300 km at 1400 LT because of 
the increase in ion drag. 

Analyses of satellite drag data show that, at least in some cases, winds in equi- 
nox blow towards the equator, not away. Calculations by Dickinson et al. show that 
during geomagnetically active times , auroral heating at bigh latitudes dominates the 
solar heating and forces circulation from pole to equator [2-23], In geomagnetically 
quiet times the flow is away from the sunlit regions. Figures 2-10 and 2-11 show 
calculated zonally averaged flow. Wind speeds at 300-km altitude and 40-deg latitude 
average 50, 20, and 5 m/sec equatorward at summer solstice, equinox, and winter 
solstice . 


2.4, 6.1 Super-Rotation 

King-Hele deduced from an analysis of satellite orbits that the atmosphere is 
rotating faster than the Earth above 140 km to approximately 450 km [2-24,2-25]. 

The average eastward winds are inferred to increase from about 40 m7sec at 200 km 
to 160 m/sec at 350 km at 30 deg latitude. Above 450 km, the' wind direction appears 
to reverse. Challinor believes that this apparent mean rotation of the upper atmos- 
phere is caused by an eastward daytime motion which is primarily confindd to equa- 
torial latitudes and notes that the satellite observations are weighted toward low 
latitude and daytime conditions [2-26]. A nxunber of other explanations have been 
offered for the inferred super- rotation of the upper atmosphere, and it may even be 
that the phenomenon is spurious. 
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Figure 2- 8. 


Contoixrs of meridional wind speed at 45 deg N latitude (m/sec) 
positive values directed toward the north [2-22]. 
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Figure 2-9. Wind vectors at the 300-km level. (Largest arrow 
represents wind speed at 225 msec ^ [2-22].) 
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Figure 2-10. Schematic diagram of the zonal mean meridional circulation 
in the thermosphere during equinox for various levels of atororal 
activity. ContoTjrs illustrate mass flow. 
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Figure 2-11. 


Same as Figure 


2-10, but at solstice 12-23]. 
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2. 4. 6. 2 Anomalous Strong Winds Associated with Geomagnetic Storms 

Increases in winds as well as atmospheric density in the upper atmosphere 
during geomagnetic storms have been observed many times. On one occasion near the 

auroral oval, a horizontal wind speed on the order of 1.5 km sec ^ was observed to 
occur at approximately 150- km altitude simultaneously with a large increase in the 
ambient density during a very large geomagnetic storm. Theoretical calculations 
suggest the possibility of associated vertical wind speeds on the order of 50 to 75 

m sec Large changes in the composition of the neutral upper atmosphere have also 
been observed during geomagnetic storms. 


2.5 SUNLIGHT AND EARTHS HINE 

A satellite in Earth's orbit is illvuninated by sunlight and earthshine. Solar 
radiation is treated elsewhere (Section 1.0). Earthshine consists of reflected sunlight 
and thermal radiation from the Earth and its atmosphere. 

Thermal radiation can be expressed by blackbody times the emissivity of the 
radiating surface. At wavelengths for which the atmosphere is transparent, radiation 
at orbital altitude coming from the Earth's surface is that of a blackbody at approxi- 
mately 288°K. Where the atmosphere is opaque, radiation is received from the colder 
upper atmosphere as shown in Figure 2-12 [2-27,2-28]. 



WAVELENGTH (pm) 


NOTE: The 288°K blackbody curve approximates the radiation from the Earth's 

surface , and the 218°K blackbody curve approximates the radiation from 
the atmosphere in those spectral regions where the atmosphere is opaque. 

Figure 2-12. A typical spectral emissive power curve for the thermal 
radiation leaving the Earth [2-27,2-28]. 
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The spectral distribution of sunlight is approximately blackbody at 5 7 60 °K, 'and 
so peaks in the visible, 0.4 to 0.7 ym wavelength. The energy flux at the Earth 
(the solar constant outside the atmosphere) is 

1.97 cal/(cm^ min) = 1371 ± 5 W m ^ = ttF . 


The ratio of total energy reflected to the total incident, called global albedo, equals 
0.39. If we assume that the Earth is a uniform, gray, diffusive reflector, having 
reflectance independent of incident angle, then the spectral radiant flux reflected 
from an element of Earth's sunlit surface is 

A 2 

= 0.39 cos 8 ergs7(cm sec ster wavelength) if cos 0 > 0 , 

= 0 if cx)S 6 < 0 , 


where 

0 = solar zenith angle at element of surface 
2 

F^ = erg7(cm wavelength) from Sun . 

Obviously, tt F^ cos 0 is the incident energy per unit surface area and ttL^ is the 
energy reflected per unit area into a hemisphere. 

Actually depends in a complicated way on the angles of incidence and reflec- 
tion and upon the nature of the surface. We may define the local albedo as the ratio 
of the energy reflected at all wavelengths in all directions from a "unit area to the 
incident energy at all wavelengths. Some values are 

0. 10 to 0. 80 clouds 
0. 05 to 0, 45 land 
0.03 to 0.20 water. 


Table 2-2 shows some typical variations of local albedo and thermal emission encoun- 
tered in Earth orbit over various time periods. 


2. 6 METEOROIDS 

The solid objects encompassed by the term "meteoroids" range in size from 

microns to kilometers and in mass range from < 10~^^ g.to > 10*^® g. Those less 
than 1 gram are often called "micrometeoroids." If objects of more than approximately 
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TABLE 2-2. VARIATION OF THE TIME AVERAGE EFFECTIVE 
VALUES ABOUT THEIR GLOBAL VALUES (2-2 7) 



Earth Emitted Thermal 



Radiation Variation 

9 

Time Increment 

Albedo Variation 

W/m^ 

At < 0.3 hr 

„ +0.30 

+»•*« -0.16 

237 

-97 

0.3 < At < 3 hr 

0.30 ± 0.10 

237 

-48 

At > 3 hr 

0.30 ± 0.05 

237 ± 21 


- fi 

10 g mass reach Earth's atmosphere they are heated to incandescence, producing 
the visible effect called a "meteor.” If the initial mass and composition permits some 
of the original meteoroid to reach Earth's surface unvaporized, the object is called a 
"meteorite." 

Meteoroids are thought to derive primarily from comets and asteroids with peri- 
helia near or inside Earth's orbit. The original objects were supposedly broken down 
into a distribution of smaller bodies by collisions. Meteoroids recently formed still 
tend to be concentrated near the orbital path of their parent body. These "stream 
meteoroids" produce the well known meteor showers which occur at certain dates and 
from particular directions (Table 2-3). 

Meteoroids may be classified by composition: stony, iron, and, perhaps, icy. 

From their composition the type of parent body can be inferred. 

Meteoroids are attracted by the Earth's gravity field so that the flux from 
allowed directions in near-Earth orbit is increased by approximately 1.7 over the 
interplanetary value. The Earth also shields certain arrival directions. 

The total mass infall to Earth is estimated to be approximately 10^^ g/year. 
Figure 2-13 shows the distribution of number with mass, where N (> m) is the number 
flux with mass > m [2-29]. The flux is low and, therefore, difficult to measure. 
Evidence includes: spherules on the sea floor and the polar icecaps, impacts detected 

with special sensors on satellites , meteor trails in the atmosphere observed visually 
and by radar, lunar crater counts, and zodiacal light. 

The fluxes of Figure 2-13 are probably uncertain by a factor of 10. The units 

2 13 . 

may be converted to particles/m sec by division by 3.155 x 10 = antilog 13.499. 

2 

(To convert to interplanetary intensity, particles/m sec ster, multiply by 2 to correct 
for Earth shielding, divide by 2 to correct for gravitational focusing, and divide by 

TT.) The data are of the form N (> m) = const /m'^, with a slowly changing. There is 
some evidence that the flux in Earth-lunar space is greater than the general level 
along the Earth's orbit by a factor between 1 and 2. The interplanetary flux is 
higher in the asteroid belt than at 1 AU. 
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TABLE 2-3. MAJOR METEOROID STREAMS [2-14] 


Name 

Period of 
Activity 

Date of 
Activity 

F 

max 

Maximum 

Geocentric 
Velocity 
(km /sec) 

Quadrantids 

Jginuary 2 to 4 

January 3 

8.0 

42 

Lyrids 

April 19 to 22 

April 21 

0. 85 

48 

r|-Aquarids 

May 1 to 8 

May 4 to 6 

2.2 

64 

0-Cetids 

May 14 to 23 

May 14 to 23 

2.0 

37 

Arietids 

May 29 to June 
19 

June 6 

4.5 

38 

5-Perseids 

June 1 to 16 

June 6 

3.0 

29 

3-Taurids 

June 24 to 
Jiily 5 

June 28 

2.0 

31 

5 -Aquarids 

July 26 to 
August 5 

July 8 

1.5 

40 

Perseids 

July 15 to 
August 18 

August 10 to 14 

5.0 

60 

Orionids 

October 15 to 25 

October 20 to 23 

1.2 

66 

Arietids , southern 

October through 
November 

November 5 

1.1 

28 

Taurids, northern 

October 26 to 
November 22 

November 10 

0.4 

29 

Taurids, night 

November 


1.0 

37 

Taurids, southern 

October 26 to 
November 22 

November 5 

0.9 

28 

Leonids , southern 

November 15 to 
20 

November 16 to 17 

0.9 

72 

Bielids 

November 12 to 
16 

November 14 

0.4 

16 

Geminids 

November 25 to 
December 17 

December 12 to 13 

4.0 

35 

Ursids 

December 20 to 
24 

December 22 

2.5 

37 


F_ is the ratio of average maximiim cumulative stream to average sporadic flux 
niHx 

for a mass of 1 g and a velocity of 20 kmisec. 
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Figure 2-13. Terrestrial mass-influx rates of meteoroids. N is the flux 
of particles with mass greater than m 12-29], 

Figure 2-14 shows a recent compilation of data for near-Earth space derived by 
various means, over a more restricted mass range than Figure 2-13. (The fluxes 

-12 

shown in Fig-ures 2-13 and 2-14 are 1-year averages.) The flux for m < 10 g is 
rather uncertain. There have been other recent estimates of micrometeoroid flux a 
factor of 10 higher than those shown in Figure 2-13 [2-30]. This appears to repre- 
sent a real uncertainty. 


2.7 PLASMA AND ELECTROMAGNETIC FIELDS 
2.7.1 Magnetic Field 

The Earth's internal field is measured on the surface and by low-altitude satel- 
lites. It changes measurably over a period of a few years and its strength and rate 
of change have been remeasured several times in tbe past 15 years. Precision, as 
well as temporal and spatial resolution, have improved as field mapping finds increas- 
ing application in geophysical prospecting and other areas. 

For the present purpose the internal field from the surface to approximately 
2000-km altitude may be represented by the gradient of a scalar magnetic potential 
(because the field has zero curl in this region of space) . 
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Figure 2-14. Cumulative particle fluxes from various data soxirces [2-31]. 
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r = radial distance 


0 = CO -latitude 

(() = longitude east . 


The strength of the multipoles is given by the Schmidt coefficients, ^ ^n*^‘ 

Detailed discussions, values of the coefficients up to n = 10, their time derivatives, 
and computer programs for nivnerical evaluation of g are given in References 2-32, 
2-33, 2-34, and 2-35. The residuals between the field so calculated and the measured 

surface field are approximately 180 y (1 y = 10 gauss) and 0.15 deg of inclination 
and declination [2-36]. Valties of the first few coefficients in gammas according to the 
1965 IGRF fit are 


n 

m 

&n 

h 

n 

1 

0 

-30339 

— 

1 

1 

-2123 

5758 

2 

0 

-1664 

— 

2 

1 

2995 

-2086 

2 

2 

1567 

130 


The dipole terms dominate even at r -s: R . 

e 

Note that these coefficients are written for use with tlie Schmidt functions, 

which have a different normalization than the usual Legendre pol 3 momi£^s P . The 
relationship is ™ 


P ™ = 

n 

[- 

(n - 
(n + 

ni)!l^ p 
m)! J nm 

"m = 2 

for 

m 

^ 0 

= 1 

if 

m = 

= 0 . 


m 

The mean square value of integrated on a sphere is [2n 1] 

The components of the field are: 

B = i- ^ 

north" r 9 9 

B ^ L_ 

east r sin 6 9 (f) 
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Also the major dipole term is = cos 0. 

To visualize the field, it is convenient to choose a new spherical coordinate 
system as follows: 

The geographic system is rotated, so that the new polar axis coincides with 
the "dipole axis" making g^' = h^’ = 0 in the new system. 

The system is translated to minimize the n — 2 term, leaving n = 1 
unchanged. (Higher order terms are altered, however.) 

This is called the eccentric dipole coordinate system. The rotation carries the pole 
11.5 deg from the rotation axis in the meridian of 69 deg W longitude. The transla- 
tion carries the origin 342 km in the direction of 6.5 deg N latitude, 161.8 deg E 

25 3 

longitude. The dipole strength is approximately 8.1 x 10 gauss cm . 

In describing the magnetosphere it is convenient to label the field lines with 
their equatorial crossing distance in units of Earth radii; this is called the "L value," 
and it may be defined to take account of the real internal field including the higher 
order multipoles. It is also useful to define the invariant latitude. A, at which a 
field line reaches Earth's surface. Using its real L value and the dipole relationship 


1 

cos A ■= 

/L 


Beyond an altithde some fraction of an it is necessary to include the effects 

of the external currents; magnetopause surface currents, the distributed ring cur- 
rents, and tail currents. Attempts to calculate these in a self-consistent manner 
including plasma convection are of great research interest at present [2-3,2-21]. 
Various mathematical formalisms have been devised for representing the external field 
mialytlcally. (Like the scalar potential for the internal field they can represent the 
field but do not specify the source currents uniquely.) 

A useful procedure has been given by Olson for calculating the quiet-time field 

inside the magnetopause on the dayside and inside the tail out to approximately 

200 R [2-37]. Updated models allowing for dipole tilt have since been derived [2-5]. 

® 1 
A computer program to calculate B can also be provided by Olson. Figure 2-15 shows 

the strength of this external field. 

An abbreviated form of the representation is 


B = 0.8535Z - 0.00507 xz + (2.1143z - 5.246 xz) exp (-0.06 r^) 


1. Dr. W. P. Olson, McDonnell-Douglas Aeronautics Co., Code A- 3- 208-13-3, 5301 
Bolsa Avenue, Huntington Beach, CA 92647. 
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Figure 2-15. The AB along the Sun-Earth line [2-37]. 

By = - [0.01382 + 67.7682 exp (-0.06 r^)]yz 

B = 8.196 + 1.363X - 0.01524z^ 0.0088y^ + 0.0609x^ 

z 

+ [-8.803 + 3.578X - 1.7927z^ - 6.0491y^ - 6. 8595x^] exp (-0.06r^) 


To this must be added the internal field for which a dipole antiparallel to the z axis 
is adequate. The coordinates are 


X toward the Sun 
z toward the north ecliptic pole 
y toward dusk meridian. 


The origin is at Earth’s center and the units are R^. 
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2.7.2 Electric Field 


In the presence of a magnetic field the electric field that is observed depends 
upon the frame of reference in which the field is measured. If E^^ and are observed 

in frame 1 and E 25 §2 in a second frame moving with velocity y << c with respect to 

the first, then in cgs units 


«2i= 




~2il ■ Bill 


where II and i are measiired with respect to B. In more useful units the first equa- 
tion becomes: 


y [m/sec] 

~2i “ §11 [volts/m] + — - X B Igauss] . 


Since above a few hundred kilometers altitude the plasma is in most cases free to 
move, making E = 0 in its own rest frame, in any other frame with velocity v 


V 

Ef = — 2 X B 
10 ^ ~ 


E 


II 


0 . 


Thus a spacecraft moving with y develops a potential across it and may draw current 
from the plasma depending upon its resistance and the ability of the plasma to supply 
current (paragraph 2. 9). Energy dissipated by the current is at the expense of the 
kinetic energy of the satellite. A satellite in low-Earth orbit woiild develop (8 x 

10^/10^) X 0.5 — 0.4 V7m, for example. 

In the lower ionosphere the high collision rate partly ties the plasma to the 
neutrals so that in the neutral frame of reference an electric field may exist. Then 
an electric current J flows; given by 


Ej^ X B 

“p 5i" ®H * '’ii!ii ■ 
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where e and are the Pedersen and Hall conductivities and a., the intrinsic or 
p H II 

parallel conductivity [2-38]. The perpendicular current occurs because the different 

ratios of cyclotron to collision frequency for ions and electrons makes the particles 

drift differently. Figure 2-16 shows typical values of conductivities. 

Since Ey = 0 (except perhaps in certain time- varying cases and above auroral 

arcs) the perpendicular field maps along B so that the electric potential between two 
field lines is constant along their length, regardless of their separation. Thus, the 
electric field known in one plane, may be mapped elsewhere. 

In an Earth- centered but non-rotating coordinate system it is convenient to 
regard the field as having three parts: 

a) Dawn to dusk polar cap field. A few tens of millivolts per meter gives 
approximately 50 kV potential across the polar cap from dawn to dusk, driving the 
magnetosphere. 

b) The plasma sheet convection field which maps into a dusk-dawn field tens 
of millivolts per meter just outside the polar cap and corresponds to Earthward con- 
vection of plasma. Field strengths are less than 1 mV/m. 

c) The co-rotation field (small) inside L = 4 to 6 corresponding to the rigid 
rotation with the Earth of the plasma in this volume.. 



Figure 2-16. Typical variation of the ionospheric conductivities 
with height for a nighttime ionosphere 12-38]. 

The total field is highly variable during substorms. It is generally stronger 
in the polar cap and plasma sheet during storms. There are various analytic models 
of the field giving both spatial and temporal dependence 12-5]. Reference 2-3 shows 
schematic diagrams of the various components. 
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2.8 CHARGED PARTICLES 


Ions and electrons occur in the magnetosphere with energies from approximately 
1 eV to many MeV, and dividing them into groups is somewhat artificial. However, it 
is convenient to consider the lower energy particles that may be characterized by 
number density and temperature (sometimes several temperatiu’es) separate from higher 
energy particles that have non-thermal distributions and are best characterized by 
intensity or flux. 

The effects of particles on a spacecraft can be extremely complicated and extended 
discussion is beyond the scope of this handbook. They are given briefly as follows: 

a) Energetic particles cause ionization to their penetration depth with damage 
a possible consequence. 

b) Plasma density and temperature (average energy) determine the potential 
to which the spacecraft charges (typically kT negative) , the current that may be 
drawn to or from a surface, and the sheath and wake dimensions (paragraph 2.9). 

c) Electromagnetic waves with frequency near or less than the electron plasma 
frequency 


or f — 9000 / N Hz , (N = electrons /cm ^) 
e e 




pe 


4ttN e 
e 

£ m 
. o e 


will propagate only in certain modes and directions if at all. 


2.8.1 Low Energy Plasma — Ionosphere 


The ionosphere extends from some 60-km altitude out to tbe plasmapause, the 
boundary of the plasmasphere on the L = 4 to 6 shell. The ionization is created in 
the upper atmosphere as the result of solar radiation and cosmic rays interacting with 
the atmosphere to produce a plasma usually consisting of equal concentrations of 
positively and negatively charged particles. The plasma remains in a state of dynamic 
equilibrium because charged particles lost through recombination are balanced by 
newly-created particles. Variations from geographical anomalies in the Earth’s mag- 
netic field exist (Fig. 2-17). Because the prime source of energy for the ionosphere 
is the Sun, the parameters defining the ionosphere vary diurnally, seasonally, and 
with the solar cycle (Figs. 2-18 and 2-19). At high latitudes and in the plasmasphere, 
variations are caused by magnetic storms . 


2. 8. 1.1 The Ionosphere Below 1000 km 

Three regions are defined to describe prominent features of the lower iono- 
sphere. In order of increasing altitude and increasing ion concentration they are the 
D, E, and F regions (Fig. 2-18) . Large diurnal effects occur in the ionosphere. At 
nighttime, the D region virtually disappears , and a depression in electron concentra- 
tion occurs between the E and F regions. 
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2-17. Daytime equal-electron-density contours (Units: 10 /cm ) [2-39] 



Figure 2-18. Normal electron distribution, daytime [2-40], 




Figure 2-19. Normal electron distribution, nighttime [2-40]. 


The D region, the lowest ionospheric region, has an approximate altitude range 
from 60 to 85 km. The predominant ionizing agent is the Lyman alpha radiation of 
hydrogen (1216 X) with cosmic rays contributing at lower altitudes. Nitric oxide 
(NO^) and oxygen (O 2 ''') are the principal ionic constituents, together with heavy ion 

complexes involving water vapor (Figs. 2-20 and 2-21). The D region has the lowest 

3 

electron density concentration of the three regions with a maximum of 10 electrons/ 

3 

cm near 80 km. 

The approximate altitude range of the E region is from 85 to 140 km. Ultra- 
violet and soft X-radiation are the principal ionizing agents in this region. Precipi- 
tating auroral particles are important at high latitudes. The electron concentration 

5 3 

in the E region ranges from approximately 10 electrons/cm during solar minimvim to 
a value approximately 50 percent larger during solar maximum. 

The F region has an approximate altitude range from 140 to 1000 km and in day 
time has two divisions, F^^ and F 2 . The solar spectrum from 200 to 900 X is the 

principal ionizing agent in this region. The two divisions occur because the degree 
of ionization and recombination falls off with decreasing altitude at different rates. 

The F^ region is associated with the ion production peak occurring in the vicinity of 

150 km. The Fj^ region disappears at night as the concentration of electrons 

decreases above the E region. The F 2 region is associated with the peak in the 

electron density distribution which varies with time of day, season, solar cycle, and 
latitude. It usually lies within the region from 200 to 400 km. The predominant ions 
near the low -altitude boundary are NO"^ and 02 "^; a gradual transition occurs with 

increasing altitude, until O'*' becomes the principal ion at the upper boundary. 
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Figure 2-20. Ionic composition of the solar minimum 
da37time winter ionosphere [2-14]. 
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2-21. Relative concentration of ions for solar minimum and a quiet 
summertime, (Nighttime summer profiles are also appropriate 
for winter daytime.) [2-14] 
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The section of ionosphere above the peak of the F2 region is commonly referred 

to as the topside ionosphere. Information on the topside ionosphere has been obtained 
mainly through satellite measurements and the incoherent backscatter radar technique 
[2-41]. 


2. 8. 1.2 The Ionosphere Above 1000 km 

The ionosphere above 1000 km includes the upper portion of the topside iono- 
sphere and the plasmasphere (or outer ionosphere) extending to altitudes of the order 
of 3 to 4.5 Earth radii. From this altitude to about 10 R is the region of transition 
to the interplanetary medium. 

Most information on the outer ionosphere has come through experimental research 
with the aid .of whistlers (RF electromagnetic signals generated by some lightning dis- 
charges in tropospheric thunderstorms) although in situ measurements by satellites 
have also provided useful resiilts. The data indicate a steady decrease in electron 
concentration with increasing altitude , and a rapid and abrupt concentration decrease 
of three to four orders of magnitude at 15,000 to 25,000 km above the magnetic 
equator. This decrease in electron concentration is called the "knee" or plasmapause 
12-42]. 

The predominant ions in the outer ionosphere vary with altitude. In the lowest 
part of the region, atomic oxygen (O"*") predominates, whereas atomic hydrogen (H"^) is 
the main constituent above 1000 km (Fig. 2-21). The electron and ion densities 
balance to give charge neutrality. 


2. 8. 1.3 Ionospheric Temperature 

The temperature profile of the ionosphere above the E region has been obtained 
by sounding rockets, satellites, and backscatter sounders.. The temperature of the 
ionosphere depends on the balance between the various heating and cooling processes 
there and the diurnal and seasonal variation of this balance. The daytime tempera- 
ture of the ionospheric electrons in the F^^ region is generally larger than the ambient 

gas temperature by a factor of two. This results from the relative efficiency of the 
ambient electrons in removing excess kinetic energy from the energetic photoelectrons 
that are continually released on the dayside and the relatively small number of elec- 
trons that share this energy. Below about 250 km the resulting hot electron gas is 
cooled mainly by collisions with neutral particles at lower altitudes and by collisions 
with ions at higher altitudes. The temperature of the ions approaches that of the 
electrons above 600 km; both exceed the neutral gas temperature. 

The model presented herein is based on experimental results for altitudes below 
1000 km. 

Temperature measurements made by the IMP 1 Spacecraft in the region from 2 
to 5 Rg show that the electron temperature, T^, is approximately proportional to the 

square of the radial distance from the Earth. Here, a value of 2,4 for the power of 
the radial distance is used to join lower altitude daytime temperature to data obtained 
from the OGO-1 Satellite. A value of 1.90 was used as the power to match nighttime 
values of T^ with the low temperature readings. The electron temperature profiles 

for daytime summer solar minimum conditions at altitudes below 1000 km are presented 
in Figure 2-22; the ion temperature profiles are presented in Figure 2-23. 
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Figure 2-22. Electron temperature, thermal velocity, and 
satellite /thermal velocity ratio [2-14]. 
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Figure 2-23. Ion temperature, thermal velocity, and 
satellite /thermal velocity ratio 12-14]. 




Above 1000 km, the ion temperature is taken as equal to the electron tempera- 
ture because the ion temperatures theoretically approach the electron temperature 
values. Figure 2-24 gives the electron temperature profiles during periods of minimum 
and maximum solar activity. 
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Figure 2-24. Model of electron temperature profiles [2-14]. 

The flux of the ionospheric particles to a spacecraft is important in calculating 
sheath and charging effects., , The fluxes determine the maximum currents of either 
sign that may be drawn. An isolated body goes to a sufficiently negative potential 
to reduce the electron current to equal the ion current (paragraph 2.9). Particle 
fluxes for electrons and ions below 1000 km are presented in Figure 2-25 and fluxes 
above 1000 km are presented in. Figure 2-26. For electrons, the nximber of particles 
flowing through a unit area is based on the electron thermal velocity, ; for ions, 

^ e 

the number is based on the satellite orbital velocity, V (circular orbit), because it 
is generally greater than the ion thermal velocity. 
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The topside ionosphere and plasmasphere are also populated by Van AUen belt 
partj^cles (paragraph 2. 8.2). 


2.8. 1.4 Plasma Sheet 

Beyond the plasmapause lies the plasma sheet, sometimes called the auroral 
domain because it extends at high latitudes to the diffuse auroral zone (Figs. 2-1 and 
2-2). The plasma density is low, but the average energy is higher than in the 
ionosphere-plasmasphere, so it is convenient to give the flux of particles as in Figure 
2-27, instead of density. 

The environment at geosynchronous orbit, 6.6 radius, is of special interest 

because many satellite systems operate there. At this distance the environment varies 
with both local time and magnetic activity (substorm or world-wide storm), and a 
satellite may find itself on either side of the plasmapause or even outside the magneto- 
pause on the day side as these boundary surfaces move. 



Figure 2-27. Composite electron and proton spectra in the inner 
magnetosphere and plasma sheet 12-43,2-44]. 


2.8.2 Energetic Particles 

The energetic particles found in the magnetosphere are the geomagnetically 
trapped or Van Allen radiation, together with the still more energetic cosmic rays 
that enter from outside. 
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2. 8. 2.1 Trapped Particles 

The trapped particles, accelerated within the magnetosphere, move as follows: 

a) They rotate arotand the field at the cyclotron frequency. 


00 = eB7mc (f^^ = 2.8 [MHz7Gauss] • B [gauss]) 


b) The (guiding) center of gyration moves in latitude along the magnetic field 
lines until the mirror point is reached (the point of maximum |b|). The particle is 
reflected and drifts to the mirror point in the outer hemisphere. (If a mirror point 
were in the atmosphere the particle would be lost through collisions.) 

c) - The line on which the guiding center moves drifts in longitude, electrons 
east and positive particles west, always crossing the equator at the same distance (in 
the dipole approximation). Thus a particle’s guiding center moves on a closed surface 
or L-shell. Similar motion occurs in the real geomagnetic field, with guiding centers 
following irregxilarities. 

In this motion kinetic energy is conserved so long as B is static. Temporal 
variations of the field accelerate and decelerate particles and also move them from one 
location to another. Another quantity rigorously conserved is directional intensity, 

2 

(particles per cm sec sterad keV). If we follow a group of particles along their 
trajectory, the directional intensity is constant wherever tliey go, in agreement with 
Liouville's theorem. 

This motion may be described by the adiabatic invariants which are "almost 
constants" of the motion. The first invariant is 


mv 


i 


B 


const — mv 


(#) . 


since v = const , 


with a the angle between velocity and field. Since a = 90 deg at the mirror points. 


. 2 

sin a 
B 


B 


mirror 


gives the mirroring field. The second and third invariants constrain the longitudinal 
drift around the Earth [2-28]. 

Particle lifetime ranges from hours (even less than a full revolution around the 
Earth) to years. Radiation may be durably trapped, if of proper energy , inside L = 
10. Beyond that, drift in longitude would carry it through the magnetopause on the 
day side. Because of the trapped particles’ motion, the distribution of omnidirec- 
tional flux as a function of energy in a single meridian plane defines the omnidirec- 
tional flux at all longitudes. Similarly, the directional intensity as a function of energy 
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and pitch angle along a single radius in the equatorial plane determines the intensity 
everywhere. 

Models of the trapped radiation have been compiled by Vette and others at the 
National Space Science Data Center 12-10,2-43,2-45,2-46,2-47,2-48,2-49,2-50]. These 
models provide numerical and graphical data and computer programs which are avail- 
able from the NSSDC, NASAVGoddard Space Flight Center to calculate flux at any 
position. The "fluence" (time integrated flux) for any orbit may be calciilated from 
these data. 

Some representative omnidirectional fluxes in a meridian plane, shown in Figures 
2-28 through 2-31, are from these models. Figure 2-32 shows differential energy 
spectra. 



Figure 2-28. R-X; map of AE-4 electron fluxes for Epoch 1964, 
threshold energy 0.5 MeV. Solar minimum. 
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Earth Radii 







Figure 2-31. AP8MIN R-A plot of constant intensity flux contours 

with an energy of > 1.0 MeV. 



C (MeV) 


Figure 2-32. Model electron spectra together with experimental data 
for geosynchronous orbit [2-51]. 







Chan et al. [2-43] have calculated fluxes of electrons and protons with E < 

100 keV encountered in five regions along the ISEE satellite orbit (21 apogee, 280 

km perigee, 28.5 deg inclination). Figures 2-33 and 2-34 show fluxes, spin averaged. 

The intensity received along a satellite orbit may be integrated over an orbital 
period or other interval. Figure 2-35 shows the orbit averaged spectra for ISEE, and 
Figures 2-36 and 2-37 show the omnidirectional electron and proton fluxes received in 
low altitude orbits per day [2-46]. The asymmetry of the geomagnetic field brings 
the radiation belts lowest over the South Atlantic anomaly as illustrated in Figure 
2-38. Vehicles at 100- to 400-km altitude will also encounter auroral radiation in the 
auroral zone at 60 to 70 deg geomagnetic latitude. The diffuse auroral radiation is an 
extension of the plasma sheet and has similar directional intensities. In the evening 
sector discrete, arcs occur with electron spectra probably resulting from localized 
electric fields parallel to B. Figure 2-39 was based on rocket measurements in aurora 
at approximately 300-km ^titude. Proton fluxes are less than 1 percent of the elec- 
tron flux. These electrons produce locally enhanced electron density in the ionosphere 
as they collide with the neutral atmosphere. Figure 2-40 shows auroral electron 
density as measured by the Chatanika incoherent scatter radar. 

At geosynchronous orbit fluxes are less energetic 12-52]. Figure 2-41 shows 
the flux there integrated over a day (Fig. 2-32). 


2.8.3 Cosmic Radiation 


Cosmic ray particles enter the magnetosphere from outside it and are of two 

types: 

a) Galactic cosmic rays — steady mean flux intensity from outside the solar 
system 

b) Solar cosmic rays — come in bursts from the Srm. 

The particles move in the geomagnetic field without significant energy change, but a 
particle of given magnetic rigidity, R, (R = P7Z momentum per unit charge) is pre- 
vented from reaching certain locations. The minimum rigidity that can reach a given 
location depends on the direction of arrival at that location. In the dipole approxima- 
tion to the field 


R 


60 cos^ A 


min 


[1 + (1 - cos^ A cos y)^]^ 


GV 


where 


r = radial distance in Earth radii 
A = latitude 

Y = half angle of allowed cone of arrival direction about the normal to the 
meridian plane. (Directions closer to the normal are allowed.) 
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DIFFERENTIAL FLUX (el«ctrons/(cma.a.sr.eV)] 



Figtire 2- 33. Composite electron spectra [ 2- 43] . 








Figtxre 2-35. Orbit -averaged differential spectra, ISEE [2-14]. 
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Figxare 2-38. Proton flux densities at am altitude of 296 km. 


INTENSITY (ELECTRONS/CM^ SEC STER KEY) 


0818:36.6 - 0818:40.7 UT 
o:0°-5° xi60°-65° 



Figixre 2-39. Auroral electron intensity (center of the arc) [2-53]. 
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DISTANCE NORTH OF CHATANIKA — km 


Figxire 2-40. Ionospheric electron density in an aurora [2-54] 




DIFFERENTIAL FLUX (particlM/em^ d«y M«V) 



ENERGY (M«V) 

Figure 2-41. Electron and proton differential energy spectra [2-14]. 
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One GV ==10 V; one GeV = 10 eV. For vertical arrival at Earth's surface 
r = 1, Y - 90 heg so 


R . = 15 cos A [GV] 

rain 


The rigidity of a particle is related to its energy [2-55] 
E = J (RZ)^ + (mc^)^ = total energy, GeV 

where 

2 

me = rest mass energy in GeV 
Z = atomic number 
R = magnetic rigidity. 


The external field alters R . so that in the polar cap, above about 65 deg latitude, 

min i- JT > o 

the cutoff is nearly zero. 


2.8.4 Galactic Cosmic Rays 


So far as we know these particles fill our portion of the galaxy uniformly and 
isotropically outside the solar system. They consist of approximately 85 percent protons, 
14- percent alpha particles, and approximately 1-percent lieavier nuclei. In addition, 
the electron and positron flux amounts to approximately 1 percent of the prcton flux. 

20 

Nuclei occur with energies at least up to 10 eV. The differential energy spectrum 

of all particles above a few GeV is proportional to 1/E^ with y approximately 2.6 
(Figs. 2-42 and 2-43). 


2. 8. 4.1 Effect of Solar Activity on Galactic Cosmic Rays at the Earth 

The intensity of particles with energy below approximately 10 GeV that reach 
Earth's orbit is reduced below the interstellar level by interaction with the solar wind 
and interplanetary magnetic field. It is thought that a cavity approximately 100 AU 
in radius from the Sun is affected. The result is that in the active part of the 11 
year solar cycle the galactic cosmic ray intensity at Earth is reduced. The effect on 
the energy spectrimi is shown in Figures 2-42 and 2-43. 

Figure 2-44 shows the dose rate versus latitude. Note that lowest energy par- 
ticles are reduced the most. One effect of this is that as cosmic rays penetrate the 
atmosphere , or a solid vehicle , the intensity and dose rate for particles of solar maxi- 

2 

mum increase up to approximately 40 g/cm shielding thickness. This is because the 
high energy particles produce many secondaries of various kinds (Fig. 2-45) [2-56]. 
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10 '^ 1 10 ^ 10 ^ 10 ^ 10 ^ 10 ^ 10 ® 10 ^ 


PROTON ENERGY (BeV) 

Figure 2-42. Galactic radiation integral energy spectrum. The flux 
is the number above the given energy [2-14]. 



Kinetic Energy, MeV 

Figure 2-43. Differential flux of galactic cosmic ray protons at 
solar maximum (1969) and minimum (1965) 12-56]. 




■"■”0 10 20 30 40 50 60 70 80 90 
GEOMAGNETIC LATITUDE (deg) 

Figure 2-44. Cosmic ray dose rate as a function of latitude at 
solar minimum and solar maximum [2-14]. 



Figtu’e 2-45. Ionization rate in the atmosphere versus depth [2-14] 


2.8.5 Solar Cosmic Rays 

Solar cosmic rays are ejected from the Sun in solar flare events and consist of 
protons, yphas, and electrons. They are generally of lower energy than galactic 
cosmic rays. Below 1 GeV kinetic energy and down to 1 MeV the integral flux of 

5 

solar protons is 10 greater than galactic particles, but is less above a few GeV. 
Table 2-4 shows some model spectra. 


TABLE 2-4. MODEL SOLAR COSMIC RAY SPECTRA [2-14] 



2 

T = kinetic energy = E - nrc 


P = RZ = momentum 


The stream of solar cosmic rays reaches and envelopes Earth minutes after a 
solar flare and reaches peak intensity in a few hours. It then decays in 1 to 2 days. 
Solar cosmic ray bursts are sporadic but are most frequent near solar maximum. 

Figure 2-46 shows dose rates from both solar and galactic particles [2-57]. 2 

Unshielded, the dose from solar particles is greater. However, a few tens of mg/cm^ 
of shielding excludes the effects of solar cosmic rays and galactic cosmic rays 
dominate [2-58,2-59]. 



Figure 2-46. Dose from solar cosmic ray events (vertical lines) and from galactic 
cosmic rays (hatched horizontal lines). Tops of lines T 30 MeV, 
triangles, T 100 MeV [2-57]. 


2-54 





2. 9 SPACECRAFT CHARGING 

In orbit, or in deep space, vehicles can develop an electric potential up to tens 
of thousands of volts relative to the ambient plasma, and large differential potentials 
can also occur. One of the consequences is electrical discharge or arcing which can 
damage surface structures and electronic systems [2-60]. Garrett recently published 
a comprehensive review [ 2- 61] . 

Garrett's review paper [2-60] gives the basic spacecraft charging equation in 
terms of number or current densities : 


J^CV) - (Jj(V) ^Jgg(V) ^Jgj(V) 

^ Jrlc(V) = 


where 


V = satellite potential 

Jg = ambient electron flux to a point on satellite surface 

Jj = ambient ion flux to a point on satellite surface 

Jgj, = secondary electron flux du« to Jg 

Jgl = secondary ion flux due to Jj 

Jrqe = backscattered electrons due to 

JpH = photoelectron flux 

J„ = active current sources such as charged particle beams or ion thrusters 

Jrlq - resistive, capacitive, and inductive coupling currents between a point 
and adjoining surface . 


Analysis of each of these terms is complex and depends upon the spacecraft 
configuration and materials, whether it is in sunlight or shadow, its altitude, and 
variable factors such as the flux of high energy solar particles and magnetic storm 
activity. 

Use of conducting surfaces where possible and adequate grounding techniques 
can significantly reduce differential charging, which is usually a more serious problem 
than is the development of a high sat ellite-to- space potential. Active techniques such 
as low-current electron and ion emitters have been used to clamp the satellite potential 
to the surrounding plasma. 
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SECTION 3. THE MOON 


3.1 INTRODUCTION 

We on Earth often see ourselves as isolated and independent of our neighbors 
in the solar system. But our life is dominated, even made possible, by the Sun, and 
it is regulated by the periodicity of the Moon. The months of our year are measured 
by the regular motions of the Moon about the Earth; our tides rise and fall primarily 
as a result of a gravitational tug-of-war between the Moon and the Earth. The 
influence of the Moon on our cultural development is woven throughout our mythologies, 
art, and literature. Indeed, the Moon is our planetary companion. 

Scientific interest in the Moon dates from the earliest recorded history (Refer- 
ence 3-1 gives a review of early lunar studies) , but it was not until the advent of 
the "space age" that we could examine the Moon firsthand. Nearly half of the Moon’s 
surface had never been seen. In 1959, Luna 3, launched by the U.S.S.R., obtained 
the first photographs of the "far side." Our understanding of the Moon has changed 
radically during the last 20 years as a multitude of spacecraft have passed, orbited, 
impacted, and landed. Although the ultimate origin of the Moon is stiU elusive, we 
are familiar with the major processes and events which shaped its evolution over the 
past 4 billion years 1 3- 2, 3- 3]. Because it has experienced little geologic activity for 
about 3 billion years, we may gain insight into the very early history of the Earth 
and other planets and the processes which shaped them by studying the ancient lunar 
surface . 

This chapter is not intended to review exhaustively the history and development 
of the Moon but, rather, to summarize its major physical properties, to give engineer- 
ing information for vehicle planning purposes , and to provide an entrance into the 
vast literature concerning the Moon. Taylor's 1975 review [3-4] is worthy of particu- 
lar attention. Other recent reviews include References 3-1, 3-5, 3-6, and 3-7. 

Results of lunar studies are published annually in the Proceedings of the Lunar and 
Planetary Science Conference. 


3.2 DYNAMIC PROPERTIES 


The lunar dynamic properties are listed in Table 3-1. 

TABLE 3-1. DYNAMIC PROPERTIES 


Parameter 

Value 

Period ; sidereal 

synodic (new Moon to new Moon) 

27. 322 days 
29. 531 days 

Eccentricity (geocentric orbit) 

0.0549 

Semi-major axis (mean Moon-Earth distance between body centers) 
perigee distance 
apogee distance 

384,402 km 
363,300 km 
405,500 km 

Mean orbital velocity 

1. 03 km sec ^ 

Inclination of lunar orbit to ecliptic 

5° 7’ 

Inclination of lunar equator to ecliptic 

1032' 

1 
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3.2.1 Comments 


The values in Table 3-1 were taken from References 3-6, 3-8, and 3-9, A 
lunar ephemeris is published regularly in the American Ephemeris and Nautical 
Almanac [3-10], The Earth and Moon make one revolution about their common center 
of mass in 27 days 7 hr 43 min 11.5 sec (the sidereal period). Tidal friction caused 
in part by the Earth's oceans is gradually slowing this revolution and causing the 
Moon to slowly spiral outward. The Moon is in synchronous rotation, keeping one 
face continually towards the Earth. Actually, nearly 60 percent of the lunar surface 
can be seen from the Earth because of the changing geometry between an observer 
and the Moon as it follows a slightly elliptic orbit inclined to the Earth’s equator. 

In addition, there is a small physical libration of dynamical origin, a deviation from 
uniform rotation which is of the order of several minutes of arc [3-11]. Reviews of 
the orbital/ tidal evolution of the Moon are contained in References 3-12 through 3-16. 
Recent laser ranging has permitted very refined orbital determination [3-17]. 

The large size of the Moon relative to its primary is unusual in our solar sys- 
tem and endows the Earth-Moon with the dynamic properties of a double -planet sys- 
tem. In addition, the Sun strongly perturbs the Moon's orbit. The solar perturba- 
tion results in significant departures from simple Keplerian motion. (Brown’s classic 
text [3-18] gives a discussion of lunar orbit theory and Deprit [3-19] gives a recent 
reAdew . ) 


3.3 PHYSICAL DATA 

3.3.1 Mass, Figure, and Other Physical Properties 

Table 3-2 summarizes the Moon's physical data. 


TABLE 3-2. PHYSICAL DATA 


Parameter 

Value 

Normalized (Earth = 1) 
or Comment 

Mean radius (R) 
Principal axes a 
b 
c 

1738 km (nominal) 
1738. 6 km 
1738.2 km 
1738. 5 km 

0.272 

Mass (M) 

99 

7.353 X 10 kg 

0.0123 

Mean density 

3.342 g cm'^ 

0.605 

Surface area 

37.9 X 10® km^ 

0.074 

Volume 

2.20 X 10^® km® 

0.02 

Surface gravity 

162.3 cm sec ^ 

0.165 

Escape velocity 

2. 37 km sec ^ 

0.212 

Gravitational 

parameter (GM) 

4902 km^ sec ^ 

0.0123 

Moment of inertia 
factor (I/Mr2) 

0. 391 ± 0.002 

(homogeneous = 0. 4) 

Bond Albedo 

0.067 

estimated [3-20] 

Visual magnitude 

-12.71 ± 0.06 

full moon 

Surface temperature 

102°K - 384»K 

[3-21] 

Atmospheric density 

5 3 

2 X 10 molecules/cm 

surface, night [3-22] 

-Bulk conductivity 

_ 5 

<10 ohm/m 
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Striking among the Moon's bulk properties is its mean density (3.34 g/cm ), 
which is lower than all of the inner rocky planets but greater than most other natural 

3 

satellites. Phobos and Deimos, the moons of Mars, have densities less than 2 g/cm ; 

3 

the moons of Jupiter, except for lo, have densities less than about 3 g/cm . The 
density of the Moon implies that it is composed predominantly of silicate minerals, and 
an iron-rich core, if there is one, must be very small. This deduction is consistent 
with the coefficient of the moment of inertia for the Moon (approximately 0.395 [3-15]), 
which implies that the mass distribution within the Moon is nearly homogeneous. (For 
a homogeneous body the coefficient would have the value 0.4.) 


3. 3. 1.1 Lunar Mapping Coordinate System 

Map coordinates on the Moon are referred to the directions of the principal axes 
of inertia. Latitude is measured from the equator positive north (towards Mare 
Imbrium); longitude is measured positive eastward (towards Mare Crisium) (Fig. 3-1). 
A small crater, Mdsting A, has traditionally been used as a reference point on the 
lunar surface, since there is no distinctive surface feature exactly at the point of 
intersection of the equator and the prime meridian. The coordinates of the center of 
Mdsting A are [3-6] : 

Longitude -5° 9' 50" ± 5" 

Latitude -3° 10' 47" ± 4" 


3.3.2 Gravity Field 


The gravity field of the Moon has been extensively studied by lunar orbiters. 
The expansion for the gravitational potential in spherical harmonics is: 



oo n ^ 

1 + Z Z (|-) p™ (sin <|))(Cj^jjj cosmA a- sinmX) 
n=2 m=0 


where cj) , A, and r are Moon-centered spherical coordinates, latitude, longitude, and 
radial distance. The m=0 terms, zranal harmonics, are written J^ (J^ = ~^n0’ 

convention). The user of expansions of this type must be careful to note that both 
normalized and unnormalized functions appear in the literature. The coefficients 
differ by a normalizing factor. Table 3-3, adapted from Liu and Laing [3-24], gives 
unnormalized coefficients to order and degree (3,3). (For values to J(15), C(8,8), 
and S(8,8) review Reference 3-14. The table presented here is truncated.) Accu- 
racy, estimated to be approximately 10 percent for 3 declines markedly for coeffi- 
cients of high order and degree. 

For all but the most specialized prirposes, only the J 2 and C 22 terms need be 
included, giving the simplified expression: 


1 J- 


(^) 


(1 - 3 sin^(})) + 3 C 22 cos^cj) cos2A 
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Figure 3-1. Simplified geologic map of the Moon showing the Apollo and Luna landing sites. 
The highland basins are divided into two groups based on their apparent ages - 
Imbrian (younger) and Nectarian (older). Mare basalts are shown in black, and a 
possible high-velocity zone at 50-km depth is indicated by stippling 13-23]. 





TABLE 3-3. SPHERICAL HARMONIC COEFFICIENTS OF LUNAR 
GRAVITY FIELD (TRUNCATED) [3-24] 


Harmonic 

Unnormalized Coefficient 

J(2) 

0.1996 X lo"^ 

C(2,l) 

0.8171 X 10'^ 

S(2,l) 

-0.7213 X 10"^ 

C(2,2) 

0.2359 X 10’^ 

S(2,2) 

0.4538 X 10'^ 


Uncertainty 

0.20 X 10"^ 

0.24 X 10"^ 
0,71 X 10'^ 
0. 53 X 10”^ 

0.60 X 10'^ 


J(3) 

C(3,l) 

S(3,l) 

C(3,2) 

S(3,2) 

C(3,3) 

S(3,3) 


0. 5878 

X 

10' 5 

0.29 

X 

10-5 

0. 3001 

X 

lO"^ 

0.27 

X 

10-5 

0.1421 

X 

10-5 

0.32 

X 

10-5 

0.4698 

X 

10-5 

0.28 

X 

10-5 

0. 5748 

X 

10-5 

0.17 

X 

10-5 

0.4847 

X 

10-5 

0.22 

X 

10" 5 

-0.2919 

X 

10-5 

0.13 

X 

10-5 


J(4) 


-0.1195 X 10 


0.17 X 10 


J(5) 


0. 4544 X 10'^ 


0.26 X 10'^ 


Local variations in gravity have been detected that in many cases correlate 
with surface features [3-25]. Areas with anomalously high- gravitational acceleration 
have been called "mascons," short for mass concentrations, on the assumption that 
the variations result from high-density materials near the surface. Mascons generally 
occur at large impact basins filled with younger basaltic lavas. The high values of 
gravitational acceleration are ascribed to a superisostatic basalt fill and7or relief on 
the crust -mantle boundary [3-26]. Geophysical models of the Moon incorporating 
gravitational, as well as topographic, seismic, and librational data have been used by 
Bills and Ferrari [3-27] to estimate that the crustal thiciness varies from 30 to 
110 hm. 


3.3.3 Magnetic Field 

The Moon does not at present possess an overall dipole magnetic field. How- 
ever, surface magnetic fields of varying direction and amplitude have been detected 
and mapped. Strength of surface fields is in the range: 6 to 313 gammas 
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(1 gamma = 10 gauss = 10 tesla). Strongest fields are associated with the ancient 
highland areas (paragraph 3.5. 1.2). 

Since the McK)n today has no dipolar magnetic field, the discovery of natural 
remnant magnetism (NRM) in samples returned from the Moon was unexpected. NRM 
is commonly acquired by rocks as they cool from high temperatures through a critical 
Curie temperature in the presence of an external magnetic field. Metallic iron, 
responsible for most of the lunar NRM, has a Curie temperature of about 750°C. 

Lunar magnetism and its origin are discussed in References 3-28 through 3-33. 

If the primordial field was produced by a dynamo within a fluid metallic core, it 
apparently switched off as the Moon cooled and convection in the core slowed. 
Estimates of the strength of the primordial magnetic field necessary to create the 
observed NRM values vary widely — from 600 to 160,000 gammas. Other suggestions 
to explain the observed fossil magnetism include (1) close approach to the Earth and 
its magnetic field, (2) enhanced solar or solar wind magnetic fields, (3) permanent 
magnetization acquired during the period of accretion from a strong magnetic field in 
the solar nebula, and (4) short-lived magnetic fields during meteorite or comet impact. 


3.4 INTERIOR 

3.4.1 Structure >• 

Studies of lunar rocks and experiments conducted by the Apollo science teams 
have given us a general idea of the structure of the lunar interior. ,The major units 
are (1) the crust (outer 30 to 110 km), (2) the mantle (extending to about 1000-km 
depth, and possibly (3) a central core (600 to 800 km in radius) (Fig. 3-2). 

Many of the major geologic processes that modify the surface of a planet are 
driven by forces arising within the body. Volcanic activity, tectonic movements 
(folding and faulting) , and a dipole magnetic field are examples of internally driven 
processes. Therefore, understanding the interior structure of the Moon would reveal 
much about its present and past geologic state. 

A wide variety of observations and measurements demonstrate that the interior 
of the Moon possesses a layered internal structure. The average density of rocks at 

3 

the surface is approximately 3.1 to 3.3 g7cm , while the bulk density of the Moon is 

3.34 g/cm . Although the difference is slight, it does require some density increase 
with depth, which is consistent with the value noted for the moment of inertia coeffi- 
cient 13-15]. (On Earth the density contrast is much larger l surface rocks have a 

3 3 

density of near 2.7 g7cm , while the bulk density of the planet is 5.5 g/cm .) 


3.4.2 The Core 


The nature, or even the presence, of a distinct central core is still uncertain. 
However, a zone of 600 to 800 km in radius does not transmit seismic S -waves (shear 
waves), and it is known that S-waves cannot be transmitted through liquids [3-34]. 
So perhaps a zone of partial melt exists at this radius — the lunar equivalent of an 
asthenosphere. Temperatures, calculated from electrical conductivity measurements, 
indicate that 1600°C may be approached at this depth 13-35]. The moment of inertia 
data permit an Fe-FeS core of 700-km radius or an iron core of 400-km radius [3-15]. 
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Figure 3-2. Major units of the lunar interior. 


If the Moon has no core, the asthenosphere is probably primitive, undifferentiated 
material that is heated by radioactive decay of K, U, and Th. 


3,4.3 The Mantle 

A sharp increase in the velocity of seismic waves marks the crust -mantle 
boundary at a depth of approximately 60 km. Seismic velocities in the lunar mantle 
(approximately 8 km7sec) exceed terrestrial values and indicate a composition rich 
in olivine and pyroxene (Mg and Fe silicates) with a density of approximately 3.4 

O 

g7cm , Petrologic studies of the lunar basalts also indicate that tkey were generated 
from such rocks at mantle depth 1 3- 4], The lunar crust and mantle comprise a 
mechanically rigid outer shell, often called the lithosphere. The Moon's lithosphere 
is much thicker than Earth's (approximately 1000 km versus 100 km). At the base 
of the lithosphere a prominent zone of moonquake activity has been detected, though 
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the periodicity of the quakes indicates that they are triggered by Earth-tides and are 
not the result of thermal convection as on the Earth [3-36]. Stress accumulation 
above a partially molten zone associated with the core may be responsible. 


3.4.4 The Crust 

The lunar seismic network established dm’ing the Apollo program provided 
detailed information about the crustal structure of the Moon in the Oceanus Procellarum 
region. Interpretations of seismic velocity profiles of the lunar crust indicate that it 
possesses a layered structure as well [3-35]. A relatively sharp discontinuity exists 
at a depth of approximately 25 km. The nature of the discontinuity is debated, but 
it appears to separate fragmented rock from underlying solid rock. Seismic velocities 
in the outer portions of the crust are low (2 to 4 km/sec) and gradually increase 
with depth. A layer of fine-grained breccia and broken rock fragments 1- to 2-km 
thick (the megaregolith) may be responsible for this part of the profile. The shattered 
rocks are probably bedrock smashed by repeated meteorite impact. The seismic veloc- 
ities continue to increase to approximately 6 km /sec at 25 km. These are similar to 
velocities measured for lunar basalts; therefore, this may be the depth of basaltic 
lavas in Oceanus Procellarum [3-35]. Alternatively, and probably more consistent 
with the surface geology, the seismic profile could be the result of a thin (approxi- 
mately 1 km) layer of basalt overlain by ejecta associated with large impact basins 
[3-37]. In either case, the increase in seismic velocity with depth can reasonably 
be attributed to the decreasing abundance of shock-induced microfractures [3-38]. 

From 25 to 60 km (the base of the crust in this region) seismic velocities are relatively 
uniform at approximately 7 km /sec and are similar to velocities determined for anor- 
thosites (a calcium aluminum silicate feldspar) and gabbros (a ferro- magnesian mineral) 
like those exposed in the lunar highlands [3-35]. The uniformity of the seismic 
signals relative to the upper section of fractured crust suggests that shock-induced 
cracks were never formed below 25 km. Alternatively, cracks may have been 
annealed by a thermal event after their formation. 


3.5 SURFACE OF THE MOON 

The surface of the Moon has two major regions with distinctive geologic features 
and evolutionary histories (1) the relatively smooth, dark areas which Galileo first 
called "maria" and (2) the densely cratered rugged highlands (uplands), originally 
called "terrae." The highlands occupy about 83 percent of the Moon's surface and 
generally have a higher elevation, as much as 5 km above the mean radius. In places 
the maria lie about 5 km below the mean radius and are concentrated on the uear 
side of the Moon (Fig. 3-1), a fact that still lacks an adequate explanation. Compre- 
hensive reviews of the geology of the Moon's surface can be found in References 3-1 
and 3-5. Much of the information regarding the lunar surface is the result of six 
manned landings on the Moon. Summaries of the geology of each of the Apollo land- 
ing sites are contained in a series of NASA Special Publications (Apollo 11, SP-214; 
Apollo 12, SP-235; Apollo 14, SP-272; Apollo 15, SP-289; Apollo 16, SP-315; and 
Apollo 17, SP-330). The Apollo missions were preceded by a series of unmanned 
Surveyor missions (NASA SP-184, 1969). Those missions, and the preceding Ranger 
program, returned considerable scientific information about the nature of the Moon's 
surface that is now often ignored because of the volume and quality of the informa- 
tion obtained by the Apollo missions. 
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3. 5. 1 Surface Morphology 


3.5. 1.1 Impact Craters 

The main external geologic process modifying the surface of the Moon is mete- 
orite impact. Craters range in size from minute pits only microns in diameter to 
gigantic basins hundreds of kilometers across which are the major land form on the 
Moon. Figures 3- 3a through 3- 3c show some of the variation in lunar craters. 

An often used classification for craters distinguishes simple craters fix)m complex 
craters and basins. Simple impact craters are the smallest, most abundant type and 
range from a few micrometers to approximately 15 km in diameter (Fig. 3- 3a). They 
are generally bowl-shaped depressions. Fresh simple craters have bright rays and 
rubbly ejecta blankets composed of the material excavated by the impact and locally 
"churned" material. Simple craters may reach depths of 2 to 3 km. Complex craters, 
which are larger than 10 to 20 km in diameter, are relatively shallower (i.e., lower 
depth-to-diaroeter ratio) than small craters. Most complex craters have scalloped or 
terraced walls, central peaks, and broad, flat floors (Fig. 3- 3b). The transition 
from simple to complex craters depends on several factors, including the impact 
energy and the strength of the material in which the crater was formed. 

Larger complex craters have distinctive central peaks, and the still larger basins 
display multi-ring features such as those around an Orientale basin (Fig. 3- 3c). 
Although these large basins are not as numerous as the smaller craters, their forma- 
tion dramatically altered the appearance of the Moon, creating huge depressions where, 
in some cases, lavas accumulated later. 


3. 5. 1.2 The Highlands 

The lunar terrae, or highlands, are the high albedo or bright regions of the 
Moon. They are primarily old surfaces that became extensively cratered during early 
lunar history when the flux of meteorites was high; they may record the final stages 
of accretion. 

The most abundant rocks in the highlands are anorthositic gabbros, composed 
principally of plagioclase , olivine, and pyroxene. Almost all the samples have been 
extensively brecciated, metamorphosed, and chemically contaminated by repeated 
meteorite impact, so that few, if any, pristine samples of the primordial crust appear 
in the lunar sample collection. Nevertheless, detailed studies of these ancient rocks 
have demonstrated the importance of magmatic differentiation in their formation. 
Apparently the Moon heated tremendously druring accretion as the kinetic energy of 
impacting bodies was converted to heat. Additional heat may have been generated by 
radioactive decay within the Moon. Large parts of the Moon melted; perhaps even the 
outer several hundred kilometers became molten. As this "magma ocean" cooled, 
minerals crystallized and separated according to their densities. Dense minerals, rich 
in Fe and Mg, sank, adding to the bulk of the mantle; less dense minerals, like 
plagioclase, floated to the surface, forming anorthositic cumulates now exposed in the 
battered highlands [3-39,3-40]. The earliest known volcanic rocks — varieties of 
basalt — were erupted shortly after the crust formed. 
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Figure 3- 3a. Simple crater. Notice its bright ejecta deposit 
and rocky rim (AS 17-2301). 
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Complex crater with terraced walls and 
central peak (AS 16-19580). 















Figure 3- 3c. Orientale basin, the youngest multi ’•ringed basin 
on the Moon. It is only partially filled by mare 
lavas (LO-IV-194M). 












3. 5. 1. 3 The Maria 


The lunar maria are concentrated in topographic basins on the near side of the 
Moon. The average normal albedo of the maria is a rather low 7 percent. Some maria 
occur within large circular impact basins such as Crisium, Serenitatis, and Imbrium, 
whereas others, like Oceanus Procellarum, occupy irregfular depressions. The lunar 
maria are vast plains composed of basaltic lava flows that were erupted after the 
highlands and the impact basins formed [3-41]. The basalt in the middle of the 
basins may be several kilometers thick but is still only a small part of the crust that 
is approximately 60-km thick. 

Lunar basalts are chemically and mineralogically similar to their terrestrial 
counterparts (paragraph 3.5.4), differing only in their high abundance of Ca-plagio- 
clase and titanium and the total absence of hydrous minerals. Experimental studies 
of lunar basalts demonstrate that they resiilted from partial melting of rocks from 150- 
to 450-km deep within the lunar mantle [3-35]. These magmas rose through the 
fractured crust and were modified only slightly before they were extruded onto the 
surface . 

Lunar lava appears to have flowed from a variety of vents and fissures, but 
there are no large composite cones and shield volcanoes like those on Earth and Mars. 
Small domes several, kilometers across may have been the sources for some flows [3-42]. 
In addition voluminous flows issued from large, irregular depressions associated with 
meandering lava channels known as sinuous riUes [3-43,3-44]. Hadley Rille, visited 
during the Apollo 15 mission, a typical rille, is 135-km long and approximately 1-km wide. 
Many large lava flows have no apparent vents and may have issued from long fissures 
like the terrestrial ones that fed the Columbia Plateau basalts in Washington State [3-42]. 

Other features characteristic of the maria are narrow ridges, several kilometers 
wide, with sinuous outlines that extend discontinuously for great distances and, in 
some cases, transect highland surfaces as well [3-45], (Fig. 3-4). "Wrinkle" ridges, 
as they are sometimes called , are believed to be formed by compression as mare 
basins subsided, adjusting to the load of accumulating basaltic lavas. Others have 
suggested that the mare ridges reflect a change in the style of tectonic deformation 
as the Moon began to cool and contract after an early period of slight planetary 
expansion [3-46]. Linear rilles , narrow, fault-bounded valleys, criss-cross some of 
the maria as well. These rilles appear to be more than 3.5 billion years old and may 
be the restilt of local adjustments to the weight of lava or may have been due to 
slight planetary expansion in an initially warm Moon [3-47]. 


3. 5. 1.4 The Regolith 

The surface of the Moon is strongly brecciated or fragmented. This mantle of 
weakly coherent debris, called the regolith, consists of shocked fragments of rocks, 
minerals, and glass spherules formed by meteorite impact. The thickness of the 
regolith is highly variable and depends on the age of the vinderlying bedrock and 
the location with respect to nearby craters and their ejecta blankets. In general, 
the maria are covered by 3 to 16 m of regolith, whereas the older highlands have 
developed a "soil" at least 10-m thick [3-48], The thickness of the regolith at the 
Apollo 17 landing site (near a highlands -mare boundary) ranged from 6 to 37 m [3-35]. 
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Figxire 3-4. Wrinkle ridges and lunar rilles are prominent features in 
this photograph of the margin of Mare Serenitatis. The linear 
rilles are on an older surface that was partially 
buried before the ridges formed (AS-17-601). 
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Figure 3-5. Contour map of Alphonsus compiled from oblique Apollo 16 
mapping camera photographs [4-49]. 
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3.5.3 Statistics of Lunar Surface Features 

For the design of lunar roving vehicles, as well as estimates of ages of observed 
features and studies of morphological changes , it is important to have statistics of 
terrain features such as size distributions of craters and slope -probability distribu- 
tions, which are given in the following paragraphs. For further details as might be 
needed for lunar trafficability studies consult References 3-51 through 3-53. 

3. 5.3.1 Crater Size Distribution 

Lunar crater cumulative size distributions follow the equation N = K reason- 
ably well. The limiting steady-state distribution for lunar plains has a= -2, while 
some actual distributions are better represented with a — -3. Reference 3-52 gives 
the following expressions : 

a) Rough mare : 

N = 10" (D < 100 m) 

N - 10 d"^ (D > 100 m) 

b) Smooth mare: 

N = 10" (D < 40 m) 

N = (100 m > D > 40 m) 

N = (200 m > D > 100 m) 

N = 10 D"^ (D > 200 m) 

c) Some uplands: 

N = 10“^D'^ (D < 1000 m) 

N = 10^ D"^ (D > 1000 m) 

d) Slopes 17 degrees and larger: 

N is negligible. 

e) Heavily cratered surfaces : 

N = 10" (D < 10,000 m) 

f) All craters on "steady state" surfaces (fresh, young, mature, and old 
craters) : 

-1 -9 

N = 10 D 

g) All craters on "steady state" surfaces that are fFesh: 

N = 3.1 X 10"^D"^ 
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h) All craters on "steady state" surfaces that are fresh and young: 

N = 2.5 X 

i) All craters on "steady state" surfaces that are fresh, young, and mature: 

N = 5 X 

2 

where: N is the cumulative frequency of craters/m 

D is crater diameters in meters. 

If the exponent a= -2, then the area covered by craters with diameters between 
any two logarithmic intervals is the same. Thus, craters with diameters between 100 
and 200 m would occupy the same total area as those between 10 and 20 m. Figure 
3-6 shows recently compiled crater distribution data. 



CRATER DIAMETER D (km) 

Figure 3-6. Absolute cumulative size -frequency distribution 
of crater populations [3-54]. 
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3. 5. 3. 2 Slope Distributions 

Slope distributions have been derived from photo grammetric analysis [3-49] and 
from bi static radar [3-55,3-56]. (In this context, "bistatic radar" refers to the tech- 
nique of illuminating a spot on the lunar surface by a radar signal transmitted from 
an orbiter and receiving the scattered signal on the Earth.) The two techniques 
compare very well except where many slopes exceed 20 deg, a cutoff point for 
bistatic radar. The radar method appears to underestimate the roughness of upland 
terrain at small scale lengths. 

Distributions generally are approximately semi- logarithmic (Figs. 3- 7a and 3- 7b) 
[3-49]. Cumulative slope probability distributions, if assruned to be Gaussian, show 
standard deviations that are typically; 

For Maria approximately 3.7 deg 

For Highlands 4.5 to 6 deg and higher. 

However, some actual distributions depart from the Gaussian form and show structures 
characterized as "complex" (bumps and wiggles) or "tailed" (more slopes at larger 
angles than a Gaussian) (Fig. 3- 7c). 



SLOPE ANGLE, IN DEGREES SLOPE ANGLE, IN DEGREES 

(a) Lunar cratered plains. (b) Uplands. 

Figure 3-7. Slope probability distributions by photo grammetry of selected liinar 

plains and upland regions [3-49]. 
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(c) Uplands -complex. 
Figure 3-7. (Concluded) 


3. 5. 3. 3 Block Distributions 

Blbck fields present obstacles to lunar roving vehicles. In an unpublished 
report, Moore [3-51] presents detailed counts; his data are summarized in Figures 
3-8 and 3-9 [3-53] . 


3.5. 3.4 Radar Studies 

Radar signals transmitted from Earth have been used to probe the lunar sur- 
face. Tyler has suggested that the strength of the wavelength dependence in the 
reflected signal correlates with the size of craters near the upper limit of a steady 
state crater distribution on the maria and can therefore be used to determine not 
only surface roughness but also the age of the surface [3-57]. Review References 
3-58 and 3-59. 


3.5.4 Chemistry of Lunar Surface 

Table 3-4 summarizes the major chemical elements found on the lunar surface. 
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Fraction of area covered by blocks larger than size 



Figure 3-8. Block distributions - craters, intercrater areas, 
and smooth Maria [3-51]. 
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Cumulative frequency (|0'^ blocks/meter^ ) 
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TABLE 3-4. CHEMISTRY OF MAJOR LUNAR ROCK TYPES [3-4] 



Highland Rocks 

Mare Basalts j 

Anorthositic 

Gabbro 

Gabbroic 

Anorthosite 

Olivine 
Basalt (A 12) 

Green 

Glass (A15) 

Si 02 

44.5 

44.5 

45.0 

45. 6 

Ti 0 2 

0.39 

0.35 

2. 90 

0.29 

AI 2 O 3 

26.0 

31.0 

8.59 

7. 64 

FeO 

, 5.77 

3.46 

21.0 

19.7 

MnO 

-- 

-- 

0.28 

0.21 

MgO 

8 . 05 

3.38 

11.6 

16. 6 

CaO 

14. 9 

17.3 

9.42 

8.72 

Na 20 

0.25 

0. 12 

0.23 

0. 12 

Kp 

-- 

-- 

0.064 

0.02 

^2^5 

-- 

-- 

0.07 

-- 

^"2"^3 

0.06 

0.04 

0.55 

0.41 

Total 

99. 9 

100.2 

99. 77 

99.4 


The chemical composition of the regolith is dominated by the composition of the 
underlying bedrock, indicating that it is derived from local sources. Although com- 
ponents from more distant sources have also been identified, about half the regolith 
appears to come from within approximately 3 km of the point of collection [3-60]. 
There is no major transfer of material across the highland-mare boundaries, and 
vertical mixing is limited [3-61], Although the upper 0.5 mm is estimated to have 
turned over around 100 times in 1 million years, there is only a 50-percent chance 
that material only 1-cm deep will have been disturbed. (For another estimate review 
Reference 3-62.) 

Turkevich [3-63] summarizes the elemental composition of the lunar surface in 
Table 3-5. 


3.5.5 Lunar Soil Particle Size Distribution 

Grain size distributions in lunar soils are summarized in Figure 3-10 [3-64]. 

The lunar regolith generally consists of particles less than 1 mm in diameter, 
about the size of silt to fine sand, with a few scattered rocks. The size distribution 
of fine particles (<1 mm in diameter) is approximately log-normal, but a slight 
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TABLE 3-5. ELEMENTAL COMPOSITION OF LUNAR SURFACE 


Element 

Average Lunar Surface 
Percent of Atoms 

Comments 

0 

61.0 


Na 

0.4 


Mg 

4.3 


A1 

9. 5 

more A1 in Terrae than 
in Maria 

Si 

16.3 


Ca + K 

6.0 


Ti 

0.3 

very little Ti in Terrae 

Fe 

2.3 

3 times as much Fe in 
Maria as in Terrae 


bimodal tendency has been noted [3-65]. The bimodality may be due to the formation 
of glassy agglutinates (60 to 500 ym) and the accumulation of very fine ejecta 
(approximately 31 ym) from distant craters. 

Lunar core samples show the regolith to consist of discrete layers deposited as 
ejecta blankets from small craters. 


PARTICLE SIZE, MILLIMETERS 



SAND-SIZE 

IcoarseImediumI fine 


SILT-SIZE 

COARSeImEDIUmI FINE 


PERCENT 
COARSER 
BY 


Figure 3-10. Submillimeter grain size distribution for Apollo 11, 
12, 14, and 15 soil samples [3-64]. 
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3.5.6 Mechanical Properties of Lunar. Surface Material 


The mechanical properties of lunar soils resemble those of terrestrial soils of 
comparable particle size, such as silty fine sand, even though composition may be 
dissimilar [ 3- 66] . 

Table 3-6 lists some mechanical properties of lunar surface material. 


TABLE 3-6. MECHANICAL PROPERTIES OF LUNAR SURFACE 


Lunar Soil 
Parameter 

Value 

Reference or Comment 

Density 

_3 

near 1 . 0 g cm 
1.5 to 2.0 

surface 

10- to 20-cm depth [3-66] 

Angle of internal 
friction 

30 to 50 deg 

higher for lower porosities 
[3-66] 

Mean porosity 

4. 3 ± 2.8 percent 

all Apollo sites for upper 
15 cm [3-67] 

Cohesion 

0.03 to 0.3 N cm“^ 

increase in cohesion for 
density increase from 0 . 99 

to 1.87 g cm'^ [3-68] 

Bearing strength 

0.02 to 0.04 N cm ^ 
30 to 100 N cm'^ 

_ 3 

density : 1 . 1 5 g cm 

density: 1.9 g cm ^ [3-69] 


Most of these measurements, obtained by the Apollo science teams, were antici- 
pated by the Surveyor missions [3-70]. 


3.5.7 Electrical Properties of Lunar Surface Material 

Electrical properties of lunar rocks are generally similar to Earth minerals 
(although a recent laboratory study found anomalously high conductivities [3-71]. 

Data on dielectric constant and loss tangent for lunar soils and rocks are given 
in Figure 3-11 [3-72,3-73]. 


3.5.8 Optical Properties of the Lunar Surface 

The normal albedo is 0.051 (Oceanus Procellarum) to 0.176 (floor of Aristarchus). 

(Note: The Bond albedo, the ratio of the light reflected into all directions to 

the incident light, an important physical quantity, is difficult to measure. Usually, 
normal albedo, the ratio at zero phase angle of the brightness of an area to the 
brightness of a Lambert disk, is all that can be specified.) 
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a. Dielectric constant as a 
function of wavelength. 
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b. Dielectric constant as a 
function of density. 
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c. Dielectric constant (solid 
circles, +) and loss tangent 
(open circles , x) for 
rock sample. 


d. Dielectric constant (solid 
circles, +) and loss tangent 
(open circles , x) for 
soil sample. 


Figure 3-11. Dielectric constant and loss tangent for lunar materials. 
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Hapke [3-20] summarizes the optical pi^perties of the lunar surface as follows: 


a) Low albedo . The Moon reflects approximately 7 percent of incident light. 
The brightest areas exceed the darkest by 3.5:1. 

b) Reddish color . The albedo increases monotonically with wavelength in the 
visible. 

c) Full Moon . There is a .strong peaking of brightness for every area at the 
time of a full Moon. 

d) Low polarization . Near quadrature, polarization is approximately 8 percent. 
Maximum variation place to place is about 3.5:1. 

e) The preceding ’’laws" apply for all known locations on the lunar surface. 

(Reference 3-20 gives a detailed discussion of phase functions, polarization, and other 
optical properties.) 


3.5.9 Temperatures of the Lunar Surface 

Since there is no moderating atmosphere, surface temperatures are extreme. 

At the Apollo 17 site the temperature ranged from 384°K during the day to 102°K 
just before sunrise [3-21]. A detailed study of infrared radiation at 11 ym from the 
eclipsed lunar surface revealed over a thoiosand small anomalously warm regions [3-74]. 
There is an approximate anticorrelation between infrared eclipse temperature and 
albedo. Local high temperature anomalies have been interpreted as due to the pres- 
ence of an excess of bare rocks on the surface [3-75]. 

_ 1 

Reference [3-53] quotes a thermal inertia (kpc) ® of approximately 800 near the 

2 I 

lunar equator, with a range of possible values 25 to 1400 cm s^K/cal. 


3. 6 ATMOSPHERE 

There is no continuous "lower atmosphere;" the Moon only has an exosphere. 
The "upper atmosphere" of the Moon, which extends to the surface, is a nearly 
collisionless gas. There are no "aeolean" effects. Measurements of atmospheric 
density were performed by the Apollo science team and indicate that the surface 

5 3 

level nighttime atmosphere contains approximately 2 x 10 moleculesicm [3-22]. 


3. 6. 1 Composition 

Table 3-7 summarizes the composition of the lunar atmosphere [3-76], 


3. 6.2 Sources of the Atmosphere 


The sources of the lunar atmosphere include (1) diffusive andVor sputter 
release of implanted solar wind noble gases and possibly some molecular gases formed 
by reduction of lunar minerals by solar wind hydrogen; (2) radioactive decay of K, 
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TABLE 3-7. LUNAR ATMOSPHERE 


Gas 

3 

Concentration (molecules/cm ) 

Day 

Night 1 

H 2 

A 

X 

0 

CO 

A 

X 

0 

“^He 

2 X 10^ 

4 

4 X 10 ^ 

20„ 

Ne 

— 

<10^ 

36. 

Ar 

— 

<3 X 10^ 

^®Ar 

— 

4 

4 X 10^ 


40 222 

U, and Th within the interior, producing Ar, Ra, and He; (3) meteorite and 
comet impact vaporization, providing water vapor and CO 2 ; (4) possibly some small 

contributions by degassing from impact events or gas venting; and (5) contamination 
by the spacecraft making the measurements [3-4]. To avoid contamination, mass 
spectrometer measurements were conducted at the Apollo 17 site at night when most 
gases would be condensed [3-22]. 

The lunar atmosphere loses helium at approximately the same rate as solar wind 
helium is intercepted [3-77]. Helitxn concentration goes through a maximum during 
the lunar night. Argon is adsorbed on soil grains at night, producing a nighttime 
minimum. Loss mechanisms include (1) thermal evaporation of gas molecules, (2) 
photo-ionization, and (3) gravitational (solar and terrestrial) perturbation of atoms 
in the atmosphere. The lunar atmosphere is in effect a reservoir providing tem- 
porary storage for free atoms that eventually escape into space. It also provides a 
channel for lateral transport of condensable gases to high-latitude cold traps where 
large deposits of water ice and other solidified compounds have probably accumulated 
over the past 2 billion years [3-78]. 


3.7 MAGNETOSPHERE 

The Moon has no magnetosphere per se. The solar wind impinges on the lunar 
surface; there is no bow shock. 


3.8 SPACE ENVIRONMENT 
3.8.1 Radiant Energy Fluxes 

The spectral distribution of sunlight is approximately blackbody at 5760°K. 

2 -2 

Energy flux at the Earth-Moon system is 1.97 cal/(cm min) or 1371 ± 5 W m . 

Other than direct sunlight , the most important source of radiant flux at the 
Moon is earth shine (sunlight reflected from the Earth) and thermal radiation from 
the Earth's surface and atmosphere. For a discussion, review Section 2.0 on the 
Earth's upper atmosphere. 


c 
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3.8.2 Charged Particle Flux 


The Moon does not have a magnetic field strong enough to prevent the direct 
impingement of ionized particles on its surface. There are three sources for the 
charged particles that have been detected at the surface of the Moon (1) the solar 
wind, (2) solar cosmic rays, and (3) high-energy galactic cosmic rays originating 
outside the solar system. 

Table 3-8 sxammarizes the charged particle environment at the Moon's surface [3-4]. 


TABLE 3-8. CHARGED PARTICLE ENVIRONMENT AT LUNAR SURFACE 


Source 

MeV/nucleon 

Proton Flux 
(Ncm 3 sec ^) 

Penetration 

Depth 

(cm) 

Solar wind 

10-3 

O 

00 

10"® 

Solar cosmic rays 

1 to 10^ 

10^ 

1 to 10" 3 

Galactic cosmic rays 

2 4 

10 to 10 

1 

1 to lo3 


3. 8. 2.1 Solar Wind Plasma 

The solar wind plasma consists mostly of H and He nuclei and, as noted earlier, 
is the dominant source of the lunar atmosphere. Surface features exposed to the 
solar wind with a velocity of approximately 400 km /sec are gradually smoothed by a 
process known as sputtering. This process poses little hazard to equipment on the 
lunar surface, however, because the erosion lifetime of a rock with a diameter of only 

- 6 “i 

10 m is estimated to be 10 years [3-79]. 


3.8. 2.2 Solar Cosmic Rays 

Solar cosmic rays are ejected from the Sun during the violent explosions asso- 
ciated with solar flares. Long-period variations in solar cosmic rays correlate with the 
well-known 11- year cycle of sunspot activity. The nuclei of light atoms (H, He) are 
dominant, but heavier nuclei (Ca, Fe) also occur. Since there is no blanketing 
atmosphere, these high-energy particles reach the surface unhindered, penetrating 
typical lunar surface materials to a depth of 1 cm [3-4]. Studies of tracks (revealed 
by etching rock samples with acid) and the isotopic changes produced by this bom- 
bardment have been used to determine deposition and erosion rates on the lunar 
surface, as well as to determine the nature and extent of mixing within the regolith. 
Preserved in these lunar rocks is a record of the Sun's activity, and studies indicate 
that the solar cosmic ray flux has been fairly constant over the past 10 million years 
[3-80] , although short cessations on the order of 50 years in length may have 
occurred [3-81]. 


3. 8. 2. 3 Galactic Cosmic Rays 

The most energetic charged particles reaching the Moon's surface come from 
outside the solar system in the form of galactic cosmic rays. Galactic cosmic ray 
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particles (atomic nuclei) have lunar penetration depths that exceed 1 m, but the 
flux rate is substantially lower than for solar cosmic rays. 


3.8.3 Meteoroids 

The Moon is bombarded with meteoroids and micrometeoroids. With no protec- 
tion from an atmosphere, impact velocities range from 2.4 to 72 kmjsec, even for 
small objects [3-82]. Micrometeoroids compromise three classes (1) cometary debris, 
(2) interstellar grains, and (3) Ixmar ejecta. 


3. 8. 3.1 Meteoroid Hazard 

The hazard posed to astronauts on the surface by meteoroid bombardment is 

2 

small [3-83]. The impact rate per cm appears to be between 1.1 and 50 craters per 
million years for pits greater than 500 ym in diameter [3-84]. 


3. 8. 3. 2 Flux of Larger Meteoroids 

The direct measurement of the meteoroid flux has been accomplished for par- 
ticles larger than about 100 g by using the lunar seismic network. Between 70 and 
150 events were recorded each year until the experiment was discontinued. The 
present-day flux has been described by the equation: 


log N = -1.62 - 1.16 log m 


2 

where N is the number of impacting bodies per km per year with a mass greater 
than m [3-36]. The distribution appears to be the result of two populations of 
impacting bodies, one that varies little throughout the year and larger objects that 
intersect the Moon's orbit during April, May, and June, producing a higher than 
average flux. A correlation between known meteor showers and the flux rate has 
been proposed [3-85]. (For more details on meteoroids, review Section 2.6.) 


3. 8. 3. 3 Flux of Micrometeoroids 

For micrometeoroid flux in space in the vicinity of the Moon, review Section 2.0 
on the Earth's upper atmosphere. For the flux of particles ejected from micromete- 
oroid impacts on the lunar surface, review paragrraph 4.1.3. 4 of Reference 3-53. 


3.9 SUPPLEMENTARY MATERIAL 
3.9.1 Origin of the Moon 

In the excitement of the early days of space exploration many scientists believed 
that the origin of the Moon would soon be understood. However, the lunar samples 
revealed that even the Moon, as small as it is, has experienced a complex history 
that reprocessed all of its original material, obscuring its primordial nature. A better 
understanding of the Moon's evolution after its origin about 4.6 billion years ago may 
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come from additional studies of the data gathered by the Apollo and Luna missions, 
but it appears that further spacecraft landings, new missions such as lunar polar 
orbit, or new analysis techniques will be required. 

Three theories concerning lunar origins are currently popular. The least exotic 
proposal suggests that the Moon and the Earth accreted from a nebula of gas and 
dust that surrounded the primordial Sun in mudh the same relative positions that they 
occupy today. Another theory postulates that the Moon formed in a different part of 
the solar nebula but was later captured by the Earth. The Moon may have survived 
intact or it may have fragmented and then reaccumulated in Earth orbit. A third 
hypothesis suggests the Moon separated from a partly differentiated Earth during a 
rapid rotational instability shortly after accretion [3-4,3-40,3-86], 


3. 9. 2 Where Do We Go From Here? 

As our nearest planetary neighbor, the Moon should continue to receive scientific 
attention in the coming decades. Possible futxare missions include a polar orbiter 
equipped with imaging and remote sensing devices to explore portions of the Moon not 
yet observed in detail. Unmanned sample return missions (such as those sponsored 
by the U.S.S.R.) could play an important role in future exploration, while avoiding 
the risk and expense of manned missions. Visionaries urge the eventual colonization 
of the Moon and have offered plans using its natural resources as partial justification 
for such ventures [3-87,3-88]. \ 


3.9.3 Chronology of Lunar Exploration 

A complete chronology of lunar missions is found in Table 3-9. 
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TABLE 3-9. CHRONOLOGY OF LUNAR MISSIONS 


Name 

Country 

Launch Date 

Objective 

Result 

Pioneer 0 

USA 

August 17, 1958 

Lunar orbit 

Launch failure 

Pioneer 1 

USA 

October 11, 1958 

Lunar orbit 

Launch failure 

Pioneer 2 

USA 

November 8, 1958 

Lunar orbit 

Launch failure 

Pioneer 3 

USA 

December 6, 1958 

Lunar impact 

Launch failure 

Luna 1 

USSR 

January 2, 1959 

Lunar impact 

MisSed Moon 

Pioneer 4 

USA 

March 3, 1959 

Lunar impact 

Missed Moon 

Lunar 12 

USSR 

September 12, 1959 

Lunar impact 

First lunar impact 

Pioneer (P-1) 

USA 

September 24, 1959 

Lunar orbit 

Exploded on launching pad 

Luna 3 

USSR 

October 4, 1959 

Lunar flyby 

First far side photography 

Pioneer (P-3) 

USA 

November 26, 1959 

Lunar orbit 

Launch failure 

Pioneer (P-30) 

USA 

September 25, 1960 

Lunar orbit 

Launch failure 

Pioneer (P-31) 

USA 

December 15, 1960 

Lunar orbit 

Launch failure 

Ranger 3 

USA 

January 26, 1962 

Lunar impact 

Missed Moon 

Ranger 4 

USA 

April 23, 1962 

Lunar impact 

Spacecraft failure 

Ranger 5 

USA 

October 18, 1962 

Lunar impact 

Missed Moon 

U-8* 

USSR 

January 14, 1963 

Lunar lander 

Failed to leave Earth orbit 

Luna 4 

USSR 

April 2, 1963 

Lunar lander 

Missed Moon 

Ranger 6 

USA 

January 30, 1964 

Lunar impact 

Impact, but camera failed 

Ranger 7 

USA 

July 28, 1964 

Lunar impact 

Success 

Ranger 8 

USA 

February 17, 1965 

Lunar impact 

Success 

Cosmos 60 

USSR 

March 12, 1965 

Lunar landing 

Failed to leave Earth orbit 

Ranger 9 

USA 

March 21, 1965 

Lunar impact 

Success 

Luna 5 

USSR 

May 9, 1965 

Lunar landing 

Crashed on Moon 

Luna 9 

USSR 

January 31, 1966 

Lunar landing 

First soft landing 

Cosmos III 

USSR 

March 1, 1966 

Lunar orbit 

Failed to leave Earth orbit 


*Unannounced 
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TABLE 3- 9. (Continued) 


Name 

Country 

Launch Date 

Objective 

Result 

Luna 10 

USSR 

March 31, 1966 

Lunar orbit 

First lunar orbiter 

Surveyor 1 

USA 

May 30, 1966 

Lunar landing 

Success 

Explorer 33 

USA 

July 1, 1966 

Lunar orbit 

Failed to achieve lunar orbit 

Lunar Orbiter 1 

USA 

August 10, 1966 

Limar orbit 

Success 

Luna 11 

USSR 

August 24, 1966 

Lunar orbit 

Success; no imagery 

Siirveyor 2 

USA 

September 24, 1966 

Lunar landing 

Crashed on the Moon 

Luna 12 

USSR 

October 22, 1966 

Lunar orbit 

Success 

Lunar Orbiter 2 

USA 

November 6, 1966 

Lunar orbit 

Success 

Luna 13 

USSR 

December 21, 1966 

Lunar landing 

Success 

Lunar Orbiter 3 

USA 

February 4 , 1 967 

Lunar orbit 

Success 

Surveyor 3 

USA 

April 17, 1967 

Lunar landing 

Success 

Lunar Orbiter 4 

USA 

May 4, 1967 

Liinar orbit 

Success 

SiiTveyor 4 

USA 

July 14, 1967 

Lunar landing 

Communications lost 

Explorer 35 

USA 

July 19, 1967 

Lunar orbit 

Success 

Lunar Orbiter 5 

USA 

August 2, 1967 

Lunar orbit 

Success 

Surveyor 5 

USA 

September 8, 1967 

Lunar landing 

Success; soil analysis 

Surveyor 6 

USA 

November 7, 1967 

Lunar landing 

Success; soil analysis 

Surveyor 7 

USA 

January 7, 1968 

Lunar landing 

Success; soil analysis 

Zond 4 

USSR 

March 2, 1968 

Circumlunar 

Engineering Test 

Luna 14 

USSR 

April 7, 1968 

Lunar orbit 

Success 

Zond 5 

USSR 

S'ept ember 14, 1968 

Circumlunar 

Success ; first recovery 

Zond 6 

USSR 

November 10, 1968 

Circumlunar 

Success 

Apollo 8 

USA 

December 21, 1968 

Circumlunar 

First manned lunar orbit 

Apollo 10 

USA 

May 18, 1969 

Circumlunar 

Success; manned 

Luna 15 

USSR 

July 13, 1969 

Lunar landing 

Crashed on Moon 

- - .i 
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TABLE 3-9. (Concluded) 


Name 

Country 

Launch Date 

Objective 

Result 

Apollo 11 

USA 

July 16, 1969 

Lunar landing 

First manned lunar landing 

Zond 7 

USSR 

August 7, 1969 

Circumlunar 

Success 

Cosmos 300 

USSR 

September 23, 1969 

Lunar landing(?) 

Failed to leave Earth orbit 

Cosmos 305 

USSR 

October 22, 1969 

Lunar landing(?) 

Failed to leave Earth orbit 

Apollo 12 

USA 

November 14, 1969 

Lunar landing 

Success ; manned landing 

Apollo 13 

USA 

April 11, 1970 

Lunar landing 

Landing cancelled en- route 

Luna 16 

USSR 

September 22, 1970 

Lunar landing 

First unmanned sample return 

Zond 8 

USSR 

October 20, 1970 

Circumlunar 

Success 

Luna 17 

USSR 

November 10, 1970 

Lunar landing 

First "rover" 

Apollo 14 

USA 

January 31, 1971 

Lunar landing 

Success; manned landing 

ApoUo 15 

USA 

July 26, 1971 

Lunar landing 

Success ; manned landing 

Luna 18 

USSR 

September 2, 1971 

Lunar landing 

Communications lost 

Luna 19 

USSR 

September 28, 1971 

Lunar orbit 

Success 

Luna 20 

USSR 

February 25, 1972 

Lunar landing 

Returned samples 

Apollo 16 

USA 

April 16, 1972 

Lunar landing 

Success ; manned landing 

Apollo 17 

USA 

December 7, 1972 

Lunar landing 

Success ; manned landing 

Luna 21 

USSR 

January 8, 1973 

Lunar landing 

Success 

Explorer 49 

USA 

June 10, 1973 

Lunar orbit 

Success 

Luna 22 

USSR 

March 21, 1974 

Lunar orbit 

Success 

Luna 23 

USSR 

October 28, 1974 

Lunar landing 

Crashed 

Luna 24 

USSR 

August 1976 

Sample return 

Success 


CO 
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SECTION 4 MERCURY 


4.1 INTRODUCTION 

Mercury, the innermost planet, orbits the Sun at 0.4 AU. In size it is between 
the Moon and Mars. The planet is difficult to obsirve telescopically from the Earth, 
since it is never more than 28 deg from the Sun. This means that observations must 
be made either at twilight or at dawn, and the observer must look close to the hori- 
zon through a thick layer of the Earth’s atmosphere. Thirteen times a century 
Mercury transits across the solar disk. (The next two transits will be in November 
1986 and November 1999.) 

Astronomers in the past attempted to draw sketch maps of albedo features, but 
the first detailed observations came during the late 1960s when radar signals were 
scattered from the planet [4-1]. Analysis of the radar return sigfnals revealed a 
rough surface and led to the first correct determination of the rotation rate. Then 
in 1974 and 1975 NASA's Mariner 10 Spacecraft provided the first close look during 
its triple flyby encounter with the planet. Mariner 10 followed a trajectory from the 
Earth that allowed it to skim past Mercury March 29, 1974, then loop around the Sun 
making another rendezvous the following September, repeating the process for the 
final encounter March 16, 1975. (In the literature the three Mariner 10 encounters 
are designated Mercury I, II, and III.) Pictures from Mariner 10 showed a surface 
closely resembling that of the Moon. No atmosphere was found other than a minute 
amount of helium and a trace of hydrogen. A surprise was the discovery of an 
intrinsic planetary magtietic field. Figure 4-1 illustrates the Mariner 10 closest 
approach points for the three encounters [4-2]. 



Figure 4-1. Mercury flyby points. 


Because of its short orbital period, Mercury makes its closest approach to Earth 
(inferior conjunction) three times a year. The flexibility in trajectory and launch 
parameters allowed by the planet’s orbit makes it an attractive planetary probe target 
in spite of communications problems posed by its small angular separation from the 
Sun, as viewed from the Earth and frequent solar occultations , as its orbit takes it 
behind the Sun. Important scientific questions abound that deal not only with 
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the planet itself, but with the early history of the solar system, how the solar wind 
interacts with planetary magnetospheres, planetary magnetic dynamos, and many other 
matters of a fundamental nature. 


4.2 DYNAMIC PROPERTIES 

The dynamic properties of Mercury are given in Table 4-1. 


TABLE 4-1. DYNAMIC PROPERTIES [4-3] 


Parameter 

Value 

Comments 

Semi-major axis 

5.79 X lo'^ km 

0.3871 AU 

Perihelion distance 

4. 60 X lo'^ km 

0.3075 AU 

Aphelion distance 

6. 98 X 10^ km 

0.4667 AU 

Eccentricity 

0.20563 


Orbit inclination 

7. 004 deg 


Mean orbital velocity 

47. 87 km sec ^ 


Sidereal year 

87. 969 days 


Sidereal day 

58.646 ± 0.0005 days 

[4-4] 

Synodic period 

115. 88 days 

time between 
inferior con- 
junctions 

Maximum angle from 
Earth- Sun line 

28 deg . 


Obliquity of rotatior 
axis 

within 1 deg of normal 
to orbit plane 

[4-4] 

Earth-Mercury distance 



Maximiam 

20.25 X 10^ km 


Minimum 

8. 94 X lo"^ km 



4.2.1 Comments 


Prior to ground-based radar observations made in the late 1960s, astronomers 
believed that Mercury always kept the same face toward the Sun. However, radar 
scattering data revealed that the planet turns on its axis with a period of 58.65 days. 
This measurement was refined by analysis of the Mariner 10 pictures, which also 
enabled the obliquity of the rotation axis to be determined. It was found to be 


4-2 









within 1 deg of normal to the planet's orbit plane [4-4]. The rotational period is 

4 

exactly two-thirds of the orbital period — at least within one part in 10 . This rela- 
tionship arises from a resonance which, in the case of Mercury, shoidd be stable [4-5]. 

The eccentricity of Mercury's orbit (0.205) and the orbital inclination (7 deg) 
are both the largest for any planet except for Pluto. The peculiarity of the orbit has 
led to a number of analyses of its ultimate stability and to speculations about its origin 
[4-6]. Mercury's distance to the Sun varies between 46 and 70 million kilometers. 

If the planet had an atmosphere this would lead to very large seasonal changes caused 
by variation in solar radiation flux. 

Precise observations over many decades have shown that the longitude of peri- 
helion passage is processing by about 45 sec of arc per century more than can be 
accounted for classically, taking into account the perturbing effects of all the other 
planets. This precession was predicted by Einstdn's General Theory of Relativity 
and is an important piece of confirming evidence for that theory. 


4.3 PHYSICAL DATA 

4. 3. 1 Mass, Figure, and Other Physical Properties 
Table 4-2 summarizes Mercury's physical data. 


TABLE 4-2. PHYSICAL DATA [4-3] 


Parameter 

V alue 

Comments 

Radius (mean equatorial) 

2439 km 

0. 38 Earth 

Mass 

3.302 X kg 

0.055 Earth 

Mean Density 

_ 3 

5. 44 g cm 

Earth = 5. 52 
Mars = 3. 98 
Moon = 3.34 

Gravity at surface 

_ 9 

370 cm sec 

Earth = 980 
Mars = 372 
Moon - 162 

Escape velocity 

4. 25 km sec ^ 

Earth = 11.2 

(Solar mass) /(Mass of 
planet) 

6023700 ± 300 

[4-7] 

GM 

22031. 8 ± 1.0 km^ sec“^ 

[4-7] 

^2 

(8 ± 6) X 10’^ 

[4-7] 

Magnetic dipole moment 

5 X 10^^ G cm^, 

see text 

Normal albedo 

0.125 

average over visible 

Surface temperature extremes 

-173°C to 427°C 


Surface atmospheric pressure 

<2 X 10"^” mbar 

[4-8] 

Flux of Solar Radiation 

14490 W m 2 (perihelion) 
6290 W m (aphelion) 

- 







4.3.2 Comments 


The large mean density of Mercury is striking considering the small size of the 
planet and the fact that it more closely resembles the Moon and other satellites than 
it does the other planets. The implication is that it possesses an anomalously large 
concentration of iron, very likely in the form of a differentiated core some 70 to 75 
percent of the total radius. The composition of Mercury is a matter of great interest 
to planetologists because of what it may reveal about the origin and early history of 
the solar system. (It is interesting that the surface gravity of Mercury and Mars are 
nearly equal because of Mercury’s larger density and smaller radius.) 


4.3.3 Mapping 

Antoniadi and other astronomers sketched maps of the light and dark areas on 
the planet’s surface from visual observations. ^ Names suggested by Antoniadi have 
been retained by the International Astronomical Union (lAU) for albedo features [4-2]. 
The relationship between albedo variations and surface topography as shown by Mariner 
10 pictures is different for Mercury than for the Moon, where dark. areas correspond 
to the Maria. On Mercury no great albedo difference exists between cratered terrain 
and plains; regions of high albedo appear to correspond to ray systems. 

Many of the large craters have been given names — usually famous authors, 
artists, and musicians (Fig. 4-2) [4-2]. 

The prime meridian is defined by the lAU dynamically as containing the subsolar 
point at the first perihelion passage of 1950 (J.D. 2433292.63). However, the Mariner 
10 system, tied to the planet’s surface, takes the 20 deg longitude line as passing 
through a small crater designated Hun Kal., Davies and Batson [4-9] discuss this 
and point out that at present there is a discrepancy of about 0.49 deg (20 km) 
between the two systems. 


4.3.4 Gravity Field 


The gravity field of Mercury was estimated by Esposito, Anderson, and Ng [4-7] 
from Doppler and round-trip time delay radio data during the first and third Mariner 
10 encounters with the planet. The results are far inferior to what could be obtained 
from an orbiter, but an estimate of J was obtained: 


J 2 = (8 ± 6) X 10 


Terms in the spherical harmonic expansion through ^<^ 2 ^ 



were carried in 


the analysis, but reliable estimates could not be obtained. They remark that taking 
a few point masses, rather than higher order terms, dramatically reduced the Doppler 
residuals, suggesting that Mercury may possess local gravity anomalies analogous to 
the lunar mascons. More detailed information about Merciiry's gravitational field would 
be extremely useful in studies of the interior structure of the planet — for example, 
to decide whether or not it has a differentiated core, as other evidence strongly 
suggests. 


1. It should be noted, however, that the features sketched by early astronomers 
bear little relationship to the albedo features now known. 
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4,3.5 Magnetic Field 


Mariner lO's discovery of a magnetic field intrinsic to the planet (that is, not 
induced by interaction with the solar wind) was unexpected [4-10]. The maximum 
field intensity measured (encounter Merciiry III) was 400 gammas, a factor of 20 
larger than the ambient solar field [4-11]. Estimates of the dipole intensity range 
22 3 

from 2 to 6 X 10 G cm [4-12], depending upon assumptions regarding solar wind 
conditions at the time of the Mariner 10 encounters. Two mechanisms have been 
suggested (1) a planetary dynamo produced by the core, assuming that it is flmd 
[4-13] or (2) a frozen-in crustal field, the remnant of ah early dynamo or external 
field [4-14]. 

Attempts have been made to calculate quadrupole and octupole components from 
the data, but this analysis is still controversial [4-12]. The derived dipole has a 

22 3 

strength about 5 x 10 G cm , inclined about 11 deg to the axis of rotation, with 
the same polarity as the Earth's field [4-3]. 


4. 4 INTERIOR 
4.4.1 Core 

- 3 

Mercury's high density (5.44 g cm ) is similar to the Earth (5.52) and Venus 
(5.25), but considerably higher than other comparable bodies in the solar system, the 
Moon (3.34) and Mars (3.98). The implication is that the planet must be largely 
metallic with up to 70 percent of its mass in an iron-nickel core with overlying silicate 
mantle and crust. If fully differentiated, the core is estimated to extend to about 
75 percent of the planetary radius or 42 percent of its volume [4-3]. The last figure 
is huge compared to values for the Earth (16 percent), Venus (12 percent), and other 
planets. Mercury's intrinsic magnetic field may imply a differentiated metallic core 
[4-15]. Mercury's assumed core is larger than the total volume of the Moon. 


4.4.2 Mantle and Crust 

Overlying the core. Mercury has a silicate mantle. The surface appears to be 
covered by a fine-grained fragmental layer quite similar to the Moon's regolith.. 


4. 5 SURFACE 

The resiilts of the Mariner 10 imaging experiments were reported by Murray and 
others [4-16] and the planetological and geological implications have been extensively 
reviewed [4-3,4-17]. Here we attempt only to indicate a few of the major current 
ideas regarding Mercury's surface. 


4.5.1 The Mariner 10 Imaging Experiments 

During its three encounters. Mariner 10 collected more than 2700 pictures 
ranging in resolution from 4 km down to a best resolution of 100 m — approximately 
the same as the best resolution of the lunar surface obtainable with terrestrial tele- 
scopes. In all, about 45 percent of Mercury's surface was observed. Stereoscopic 
coverage was obtained for parts of the southern hemisphere., (The high-solar 
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elevation angles for some of the imagery make analysis and interpretation of topog- 
raphy and morphology difficult.) 

NASA's Atlas of Mercury [4-21 gives pictorial coverage and diseussion. Figure 
4-3, adapted from Reference 4-18, compares coverage and ground resolution for 
Mercury from Mariner 10 with the Moon and Mars. 



Figure 4-3. Ground resolution (km). 


4.5.2 Surface Morphology 

The appearance of Mercury is very similar to the Moon. Impact craters are 
numerous, with crater diameters covering a full spectrum from lunar-basin sizes like 
the Caloris Basin at 1300 km, down to the 100-m resolution limit of the imaging camera. 
There is no reason not to suspect a fiill range of small crater sizes as found on the 
Moon. The presence of a regolith (fine well-churned soil) shows that Mercury, like 
the Moon, has been bombarded by micrometeoroids throughout its history. In the 
absence of an atmosphere, micrometeoroids reach the surface, which they have grad- 
ually degraded and weathered. The degree of degradation covers a wide range and 
depends on the age of the crater and subsequent meteoric flux. But, in addition to 
weathered craters, many fresh sharp-rimmed craters are seen which still possess 
bright ray systems of ejecta, indicative of their relatively recent origin. A geologic 
map of Mercury was published by Trask and Guest [4-19]. 

Three main types of surface [4-17] are described in the following paragraphs. 


4. 5.2.1 Heavily Cratered Terrain 

The heavily cratered terrain is made up of large craters (thirty to hundreds of 
kilometers in diameter) so closely packed together that they frequently overlap. 

Inside these craters are often found smooth plains that are only lightly cratered and, 
therefore , are considered to be younger than the intercrater plains surrounding the 
large craters. Heavily cratered terrain resembles the Lunar highlands. 
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4. 5.2.2 Intercrater Plains 


The intercrater plains surround and merge with the large craters of the heavily 
cratered terrain. These plains, which cover much of the planet, have rolling slopes 
and are densely covered with small craters (5 to 10 km) which frequently are elongated 
and occur in groups or chains and, therefore, are thought to be secondary in origin. 
The intercrater plains are generally thought to be among the oldest surface features 
on the planet since they apparently were formed at a period of high meteor flux. 

They do not appear to be derived from the ejecta of the major (basin-forming) meteor 
impacts, since Mercury, unlike the Moon, has relatively few large basins. They may 
have formed later from volcanic extrusions at a time between the formation of the 
heavily cratered terrain and the smooth plains. (It is, however, not known for certain 
that such a period of volcanic activity actually occurred.) 


4. 5. 2. 3 Smooth Plains 

The smooth plains are relatively flat areas with little evidence of major crater- 
ing. They may, therefore, be among the youngest features on the planet's surface. 
They resemble the lunar Maria and have sinuous ridges similar to features seen both 
on the Moon and Mars. However, unlike the lunar Maria which contrasts distinctively 
with the highlands because of their lower albedo, the smooth plains are similar in 
brightness to both the heavily cratered terrain and the intercrater plains. 

The smooth plains often exist within and about the Mercurian basins, such as 
Caloris and the various basins in the north polar region, and in the floors of other 
large craters. The smooth plain material is thought to be derived either from volcanic 
magma, impact ejecta, or from melt. 


4. 5. 2. 4 Other Surface Features: Chaotic Terrain and Lobate Scarps 

At the antipodal point to the Caloris Basin, the largest impact structure seen on 
Mercury, is a particiilarly hummocky and disrupted region. Prominent hills and valleys 
are seen, with many broken crater rims. This region is generally interpreted as 
having been produced by seismic waves generated by the Caloris impact and focused 
at the antipodal region. 

Widespread over Mercury's surface are series of lobate scarps thought to be 
thrust or high- angle reverse faults from a global crustal shortening. Their origin 
has been extensively discussed [4-20,4-21]. The accepted Interpretation is that they 
are due to a combination of early tidal de-spinning of the planet and contraction asso- 
ciated with cooling of the planet's large core and thin mantle. (The thermal history 
of Mercury has been investigated by Toksoz, Hsui, and Johnston [4-22] . ) 


4. 5. 2.5 Other Comments on Surface Characteristics [4-3] 

a) The intercrater plains are the predominant surface feature, covering much 
of the planet so far surveyed. 

b) Albedos of the younger plains and the older features are not greatly 
different. Only a slight correlation with surface morphology has been established. 
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c) The Mercurian craters do not expel ejecta blankets or secondary missiles to 
as great distances as on the Moon, since the surface gravity is stronger. 

d) The widespread distribution of scarps seen on Mercury are not observed 
on the Moon. 


4. 5.2.6 Surface Feature Statistics 

Strom [4-3] notes that the crater di am eterV density distribution of craters in the 
Mercurean highlands is very similar to the distribution in the lunar highlands and the 
heavily cratered regions of Mars. For planning purposes, therefore, it would be 
reasonable to adopt for Mercury the statistical distributions for craters and blocks 
that have been proposed for the Moon for regions that appear similar. 

One difference is that ejecta from Mercurean craters cover a smaller area than 
those from lunar craters, due to the larger gravity [4-23]. Similarly, for Mercuty 
the areal density of secondary craters peaks closer to the primary crater than is the 
case for the Moon. 


4. 5.2.7 Volcanism 

There is no present day volcanic activity. The absence of volatiles in the 
atmosphere can be cited, along with the lack of pictorial evidence from Mariner 10. 
Certain surface features have been interpreted as impl3dng volcanism in the past [4-24] , 
although the issue is still considered open [4-25]. 


4. 5.2.8 Surface Chemistry 

There is very little information available on the chemistry of Mercury's surface. 
McCord and Clark [4-26] detected the Fe"^ ion from ground-based observations at 
0,89 microns. (This is an example of the continthng importance of telescopic observa- 
tions of the planet.) Although the surface appears similar to the Moon, color differ- 
ences suggest differences in Fe and Ti abundances [4-27]. Requirements of vapor- 
phase deposition from micrometeoroid impacts sets the FeO content of the surface 
layer between 3 and 6 percent. The absence of argon in the atmosphere could be 

40 

interpreted as indicating a deficiency of potassium in surface materials, since Ar 

40 

results from the radioactive decay of K. Such a deficiency was predicted, but it 
is premature to set a definite limit [4-8]. 


4. 5, 2. 9 Optical Properties 

The optical properties of Mercury are given in Table 4-3, 

Bright craters and their ray systems are bluer than their surroundings by 
about 12 percent. 

Mercury spectral reflectance increases linearly from 330 to 1060 nm. (terres- 
trial observations) ; the reflectance spectra for Mercury closely resemble those for 
mature lunar soils [4-28]. 
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TABLE 4-3. ALBEDO 


Region 

Typical Albedo 
at 486 nm [4-27] 

Bright craters 

0.35 

Smooth plains 

0.12 

Intercrater plains 

0.14 

Hummocky plains 

0.13 

Bright rays 

0. 22 

Flat bottom craters 

0.15 

Secondary crater fields 

0.13 

Normal albedo 

0.125 [4-3] 


The bolometric emissivity for the planet, assumed to be independent of wave- 
length and directionality, was reported by Hansen [4-29] to be 0.82 ± 0.06. (This 
determination was based on an ingenious application of ground-based observations.) 


4.5,2.10 Thermal Properties: Temperatures 

At the equator near local midnight , a minimum temperature of 100°K was mea- 
sured by the Mariner 10 infrared radiometer. At the subsolar point near perihelion 
the temperature reaches about 700°K, in excess of the melting point of lead. Because 
of the peculiar rotational and orbital relations, together with the large orbital eccen- 
tricity, mid-day lasts significantly longer at 0 to 180 deg longitudes, than at the 90 
and 270 deg longitudes. Hence, at low latitudes, mean temperature and temperature 

extremes vary with longitude. The surface temperature at the poles remains well 

_ 1 

below 0°C [4-30]. A typical thermal inertia, (kpc) was found to be 


600 cm^ s^ K/cal 


from the cooling curve. Areas were found where this parameter is as much as twice 
the quoted value. Strom [4-3] interprets these findings as consistent with a layer 
of insulating silicate dust similar to the lunar regolith, with large-scale regions that 
are more compacted, or where rock outcrops or boulders are not blanketed by dust. 
The thermal surface properties closely resemble those of the Moon suggesting similar 
soil properties [4-31]. 


4. 6 ATMOSPHERE 

Prior to Mariner 10, although ground-based observations had failed to detect 
an atmosphere for Mercury, a surface pressure of as much as 0.1 mbar could not be 
ruled out. The pressure is now known to be at least nine orders of magnitude less 
than that. Mariner 10 carried two spectrometric instruments, one for remote 
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atmospheric sensing of radiation below 95 nm during solar eclipse by the planet and 
the other for observation of emissions at discrete locations in the spectrum [4-32]. 
The result was that He and H, and probably O, were detected. For other gases that 
were not detected, upper limits were estimated. Even at the surface gas molecules 
very rarely collide, so the Mercurean "atmosphere” is really an exosphere. It is not 
possible or useful to compute model atmospheres for the planet. 


4.6.1 Pressure 


The total pressure at the surface was determined to be less than 2 x 10 


-10 


mbar. 


4.6.2 Composition 


Only helium and hydrogen (and O with less certainty) were detected directly. 

No volatiles were detected, nor was argon, a popular candidate because of its genera- 

40 

tion by radioactive decay of K. (Review the discussion in Reference 4-8.) Upper 
limits were estimated for Ne, Ar, Xe, C, CO 2 , and H 2 O. Possible constituents 

and partial pressures at the subsolar point [4-33] are given in Table 4-4. 


TABLE 4-4. ATMOSPHERIC PARAMETERS 


Possible 

Sub solar Point Density 

Partial Pressure 

Constituent 

(atoms/ cm^) 

(mbar) 

He 

4500 

5 X 10"^^ 

H 

8 (thermal component) 

82 (non-therm al component) 


Ne 


< 1.2 X 10'^® 

Ar 


< 6.9 X 10 "^^ 

0 


< 5.6 X 10 “^^ 

C 


< 1.5 X 10 


No measurable absorption was obtained in the 30 to 95 nm wavelength region 
for sunlight when Mariner 10 passed into the umbra of the planet. Most gases have 
large absorption cross-sections in this region and, thus, limits can be placed on the 
densities of many gases of interest near the surface. He, Ar, O, CO 2 , H 2 O , O 2’ N 2 . 

and H 2 densities at 8-km altitude are all set at less than 2 to 4 x 10 ^ molecules/ 

cm^ [4-33]. 

4.6.3 Origin 

Helium atoms escape thermally from the planet and must, therefore, be made up 
either by accretion from the solar wind or through radioactive decay of uranium and 
thorixim in the crust. 
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Hydrogen may be supplied from outgassing of the crust in the form of water 
vapor or hydrocarbons, or from the solar wind. A likely mechanism is photodissocia- 
tion of H^O (concentration below the detectability threshold). The source of the 

assumed HgO vapor cxjuld be outgassing, interaction of solar wind protons with the 

soil , or meteoroid flux [ 4- 8] . 


4.7 IONOSPHERE 

No ionosphere was detected by Mariner 10 radio occultation experiments [4-34], 
Upper limits : 

3 

4000 electrons7cm night side 

3 

1500 electrons /cm dayside. 

(The difference is due to h-f signal fluctuations probably not planet -related.) 


4.8 MAGNETOSPHERE 

Experiments aboard Mariner 10 showed the existence of a planetary magnetic 
field; a permanent magnetic dipole about 0.06 percent of the terrestrial strength 
[4-13]. Though weak, this field is sufficient to give rise to a magnetosphere. The 
plasma science experiment, along with the magnetometer, identified a magnetopause, 
magnetosheath, and bow-shock associated with the magnetosphere (Fig. 4-4 [4-3]). 
The presence of a magnetotail is generally accepted, although Mariner lO's trajectory 
did not allow it to collect data in this region. Time -depen dent processes similar to 
the magnetospheric substorms on Earth may have been observed [4-35]. 

The magnetosphere of Mercury is relatively much smaller than Earth's and only 
extends to about 1.5 planetary radii whereas the terrestrial magnetosphere extends 
to over 10 Earth radii. Moreover, the size and shape of the magnetosphere were 
found to vary markedly with the flux of incident solar wind. The magnetosphere 
may even be compressed to the surface at the planet at times [4-12], The weakness 
of the magnetosphere and the ease with which it distorts rule out trapped particle 
regions such as Earth's Van Allen Belts. 


4. 9 SPACE IN THE VICINITY OF MERCURY - FIELDS AND FLUXES 

A. Solar radiative flux 

Mean solar flux at 0.38 AU: 9260. W m 

NOTE: Except in the visible, solar fluxes are highly variable [4-36,4-37]. 

B. Solar magnetic field 

Twenty to sixty gammas at Mercury orbit during Mariner 10 encounters [4-12]. 
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C. Solar wind (typical values) 


velocity 

density 


550 km7sec 

- 3 

14 particles cm (Mercury I [4-12]) 

MERCURY 



Figure 4-4. Comparison of the magnetic fields of 
Mercury and Earth [4-3]. 
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SECTIONS. VENUS 


5.1 INTRODUCTION 

The second planet out from the Sun is Venus, named after the Roman goddess 
of love and beauty. Venus is similar to the Earth in size and density, but there the 
similarity stops. It rotates extremely slowly in the retrograde direction, one revolu- 
tion in 243 days. A thick atmosphere, largely carbon dioxide, blankets the planet. 

The surface atmospheric pressure is 90 times that of the Earth and the temperature 
exceeds 700°K. Because of uninterrupted cloud cover the surface can only be "seen" 
by using radio wavelengths in which the clouds are not opaque. 

From the Earth, Venus appears through a telescope to be a bright object, 
essentially featureless in visible light, with a slightly yellowish color. Infrared 
observations also fail to reveal features, but in the ultraviolet, dark markings, often 
of a characteristic horizontal Y shape, appear in the clouds that enshroud the planet 
(Fig. 5-12). 

Venus is often called the morning or evening star, since when it is prominent 
in the sky it is always seen near the rising or setting Sun. It is the third brightest 
object in the sky, behind the Sun and the Moon, attaining a visible magnitude of -4.3. 
Observed from the Earth, Venus occasionally transits across the face of the Sun, 
though less frequently than Mercury. The last transits were in 1874 and 1882, and the 
next will be June 7, 2004, and June 5, 2012. 

Venus has recently been the destination of a number of planetary probes. The 
U.S.S.R. has concentrated its planetary program on Venus; between 1965 and 1979 
12 Venera spacecraft impacted, soft-landed, or flew by the planet. The U.S. pro- 
gram has included flybys in 1962 by Mariner 2, in 1967 by Mariner 5, in 1974 by 
Mariner 10 and, most recently, the extremely successful Pioneer Venus Project. 

Pioneer Venus consisted of two spacecraft, the Orbiter and the Multiprobe. 

Both were launched in the summer of 1978 and reached the planet in December, 1978. 
The Orbiter was placed in a highly eccentric orbit about the planet. Perigee altitude 
was initially controlled to stay within the upper atmosphere, at 160 km, with apoapsis 
66,900 km, period 24.03 hr. Perigee control was halted during the summer of 1980, 
retaining sufficient propellant to allow attitude control for many more years. (Data 
continues to be received from the Orbiter.) Solar perturbations will eventually cause 
perigee to increase to a maximum of about 2200 km in 1986, with corresponding 
decrease in apogee. Then the perigee will begin to decrease and, eventually, in 
1992, the Orbiter will burn up in the venusian atmosphere [5-1]. 

The Multiprobe consisted of a main bus similar to the Orbiter, together with 
four separate entry probes : the large Sounder Probe and three other identical small 

probes (North Probe, Day Probe, and Night Probe). All of these vehicles entered 
the atmosphere of Venus on December 9, 1978. The four probes entered first, 
staggered over a 10-min interval; followed by the bus entry 90 min later. The small 
probes, falling freely, reached the surface at different locations; the bus burnt up 
at about 110 km. The probes transmitted information for about 55 min until impact, 
except for blackouts lasting about 1 min. The Day Probe survived for 67 min after 
impact. Figure 5-1 shows the locations of the probe landing points in the lAU coor- 
dinate system [5-1]. The landing points for various Venera missions are also shown, 
along with tics at the top that indicate the longitudes of the subsolar points at the 
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Figxire 5-1. Vemis landing coordinates, lAU, 


time of landings. The U.S.S.R. also took advantage of the 1978 launch opportunity, 
and the Venera 11 and 12 probes softlanded on Venus on December 25 and December 
21, 1978. 


5.1.1 General References 


A. May and June, 1977, issues of Space Science Reviews [5-2,5- 3] gave 
planning data for Pioneer Venus a year before launch. 

B. Special issues of Science, February and Jialy, 1979, gave preliminary 
results from Pioneer Venus [5- 4, 5- 5]. 

C. The most comprehensive collection of papers on Pioneer Venus experimental 
results and their implications is in the special edition of the Journal of Geophysical 
Research [5-6]. This section draws heavily on this source. 

D. U.S.S.R. Venera results are summarized by Marov [5-7]. Preliminary 
results for Venera 11 and 12 are given in COSPAR 1980 [5-8]. A comprehensive 
account of Russian and U.S. results through 1979 and most of 1980 appears in the 
review by Moroz [5-9]. 

E. The recent Scientific American article ’’The Atmosphere of Venus” [5-10]. 

F. Venus, eds. D. M. Hunten, L. Colin, T. M. Donahue, and V. I. Moroz, 
University of Arizona Press, Tucson (expected publication date: late 1982). 

Other than for historical interest, much pre- Pioneer Venus information about 
Venus must be considered obsolete. Exceptions are orbital data. Earth- based radar 
imaging, and some photometric and spectroscopic measiirements. 


5.2 DYNAMIC PROPERTIES 

The dynamic properties of Venus are given in Table 5-1. 
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TABLE 5-1. DYNAMIC PROPERTIES [5-1,5-11,5-12,5-13,5-14] 


Semi -major axis 

1.082 X 10^ km (0. 7233AU) 

Perihelion Distance^ 

1.075 X 10^ km (0.71844AU) 

Aphelion Distance 

1.089 X 10® km (0.72824AU) 

Eccentricity 

0.0068 

Inclination of orbit 

3.393 deg 

Mean Orbital Velocity: 
Linear 
Angular 

-1 

35.01 km sec 
1.602 deg day ^ 

Sidereal Year 

224. 701 day 

Sidereal Day 

243.0 ± 0.1 day (retrograde) 

Local (Solar) Day 

116.75 days 

Inclination of equator to 
the orbit 

2.2 deg (cf. Earth 23.5 deg) 

Synodic Year 

583. 92 days (time between 
inferior conjunctions) 

Earth-Planet Distance : 
Maximum 
Minimum 

2.5861 X 10® km (1.7292AU) 
0.4327 X 10^ km (0.2893AU) 


5.2.1 Comments 

At closest approach, just over 43 million kilometers, Venus is our nearest 
planetary neighbor. Mars never comes closer then 56 million kilometers. 

The spin- axis of the planet is very near perpendicular to its orbit plane, 
unlike the Earth. The spin-period, besides being very long, is retrograde (clock- 
wise , looking down on the Solar System — in the opposite sense from the orbital 
motion). Venus' rotation is almost (but not qmte) in resonance with the orbit of 
the Earth, so that nearly the same face is presented at each inferior conjunction. 

Of the planets, Venus and Neptune have the most nearly circular orbits, with 
orbital eccentrides less than 0.01. 


5.3 PHYSICAL DATA 

Table 5-2 summarizes the physical data of Venus. 
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TABLE 5-2. PHYSICAL DATA [5-9,5-12,5-14,5-15] 


Radius , equatorial 

6051.3 

(sO.95 Earth) 

Mass 

4. 870 X 10^^ kg 

(=0.8150 Earth) 

Mean Density 

_ 3 

5. 24 g cm 

(Earth s 5.41) 

Surface Gravity 

- 2 

887 cm sec 

(Earth = 980) 

Escape Velocity at the 
Surface 

10. 4 km sec ^ 

(Earth = 11.2) 

Magnetic Dipole Moment 

Probably zero 

Paragraph 5.3.3 

Albedo, Average Over the 
Visible Spectrum 

0.80 ± 0.02 


Surface Atmospheric Pressure 

93 ± 1. 5 bar 

( 0 - 95 , X Earth) 

Surface Temperature 

s740®K 


Oblateness (J 2 ) 

6 X 10'® 

(Earth = 1083 x lO"®) 

Solar Flux 

2620 + 35 W m"^ 

Range in orbit 

Bond Albedo 

0.75 ± 0.03 



5.3.1 Comments 

The atmosphere of Venus is quite unlike the Earth's atmosphere. Its unusual 
properties are discussed at len^h in paragraph 5. 6. 

The limit quoted for Venus' magnetic dipole is a limit only and not to be con- 
strued to imply that Venus has any intrinsic magnetic moment [5-16]. 

A striking feature of the planet is its extreme sphericity. 

An albedo of 0.75 is required for balance between absorbed and emitted energy 
and is thought to be acceptable within realistic error bounds. 

5. 3. 2 Gravity Field 

The global gravity field was estimated from long-term orbital perturbations of 
tile Pioneer Venus Orbiter [5-17]. The oblateness (Jg = -Cgg) is estimated to be about 

6 X 10 . Estimates of unnormalized coefficients to order and degree 6 are given in 

Table 5-3. 

For any practical orbital calculations, only J„ (= -C„„) and C„„ would be used 
(paragraph 3.3.2). ^ 
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TABLE 5-3. SIXTH -DEGREE AND SIXTH-ORDER UNNORMALIZED 
COEFFICIENT ESTIMATES AND THEIR UNCERTAINTIES [5-17] 


/ 

m 

G/,1 


■S/m 


2 

0 

-0.597E-5* 

0.32£-5 



2 

1 

0.640£-8 

0.11 £-5 

-0.299£-6 

0.13£-5 

2 

2 

-0.332S-6 

0.81 £-6 

-0.1 74 £-5 

0.74£-6 

3 

0 

0.779£-5 

0.73£-5 



3 

1 , 

0.191£-5 

0.26£-5 

0.985 £-6 

0.24£-5 

3 

'2'"‘ 

-0,175£-6 

0.51 £-6 

0.438£-6 

0.56£-6 

3 

3 

-0.150£-7 

0.1 5 £-6 

: 0.4l3£-7 

0.1 3 £-6 

4 

0 

-O.I6.5£-5 

0.42£-5 



4 

1 

O.I9X£ (. 

O.I4£-5 

0.1 27 £-6 

0. 14£-5 

4 

2 

0.42l£-<> 

0.33£-6 

0.220£-(i 

0.33£-6 

4 

1 

-0.5Wt£ -7 

0.54£-7 

-0.234£-7 

0.56£-7 

4 

4 

-0,l79£-7 

0.1 5 £-7 

0.1 73 £-8 

0. 14£-7 

5 

0 

-0.691 £-6 

0.32£-5 



5 

1 

-0.379£-7 

0.84£-6 

0.404£-8 

QME-ti 

5 

2 

O.I56£-7 

0.16£-6 

0.925 £-7 

0.l6£-6 

5 

3 

0.1l5£-7 

0.30£-7 

-0.541 £-8 

0.29£-7 

5 

4 

0.743 £-K 

0.52£-X 

0.296£-8 

O.-Mf.-S 

5 

5 

0. !02£-X 

0.13£-X 

-0.163 £-8 

0 14£-8 


0 fJ', 

-0.424£-6 

0.24£-5 



6 

1 

-0.l80£-7 

0.54£-6 

0.546£-8 

0.54£-6 

6 

2 

0.711 £-8 

0.85£-7 

0.1 59 £-7 

0.85£-7 

6 

3 

0.152£-9 

0.14£-7 

0.225 £-8 

0.14£-7 

6 

4 

0.404£-ll 

0.33£-8 

0.1 79 £-8 

0.23 £-8 

6 f 

.■•5> .. 

0.305 £-10 

0.4 1 £-9 

-0.1 38 £-9 

0.38 £-9 

6 

6 

-0. 120£-9 

0.91 £-10 

-O.I09£-9 

0.89 £-10 


• Read -().5.9,7/i'-r5 as -0.597 x 10 ' 


5.3.3 Magnetic Field [5-16] 

Venus lacks a dipole magnetic field that might be associated with a planetary 
dynamo. Reference 16 reports an upper limit for an intrinsic planetary dipole of 
1 6 3 

4. 3 X 10 G m , based on measurements made by the Pioneer Venus orbiter. The 
observations are consistent with zero intrinsic planetary moment. 


5. 3. 4 Mapping of the Planet 

The lAU has established a coordinate grid for Venus. By agreement, the lAU 
approve all feature names. The prime meridian is placed at the point that faces the 
Earth at inferior conj imction ; the bright feature, the summit of Maxwell Montes, is at 
longitude 2.2 deg. Figure 5-2 shows a contour map of Venus with identification of 
the major named features [5-18]. 


5.4 INTERIOR 
5.4.1 Core 

Venus has a, mean density similar to the Earth and Mercury. Its relatively 
high value leads to the expectation of a compacted core. 
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Figure 5-2. Contour map of Venus. 

Generalized topographic map of Venus showing elevations of surface features. Contour interval is 
1.0 km. The highest point is the summit of Maxwell Montes (11.1 km) [5-18]. 



Although the absence of a planetary magnetic field of interior origin counts 
against a liquid core, the near resonance between the spin and orbit of Venus and 
the Earth's orbit suggests that the frictional damping of a liquid core existed at some 
stage in the planet's history. Further evidence is the existence ^f basaltic rocks at 
the landing sites of some of the later Venera landers 15-7]. 


5,4.2 Crust and Mantle 


The presence of basaltic rocks at some surface locations ipdicates that the outer 
layers of Venus at some time in the past underwent volcanic episodes. 

Data from ground-based radar observations suggested an offset of the center- 
of-mass of Venus of as much as 1.5 km toward the Earth at inferior conjunction. 
However, analysis of Pioneer Venus data indicates less than a 0.5-km offset [5-15]. 

Many gravity anomalies have been noted although they are not nearly as large 
as those for the Moon or Mars. They correlate with topography [5-19] with important 
exceptions. The most prominent highland area, Ishtar Terra, does not show a 
corresponding gravity anomaly. Aphrodite Terra and Beta Regio, other highland 
regions, however , do correspond to gravity anomalies, suggesting basic differences 
between their geology and Ishtar's [5-18], 

Some linear features have been interpreted as due to local or regional tectonic 
activity [5-15]. 


5.5 SURFACE PROPERTIES 

Although the surface of Venus is obscured in visible, ultraviolet, anu i..„-ared 
wavelengths, the cloud layers are transparent at radio wavelengths and since 1961 
many radar studies have been made [5-20]. The first spacecraft close-up radar obser- 
vations were made by Venera 9 and 10 in 1975 [5-7], using a bistatic system. Since 
1978 the Pioneer Venus Orbiter's monostatic radar imaging system has been in opera- 
tion, returning altimetry and imaging data to the Earth [5-21]. Global topography of 
Venus has become available from PV Orbiter imaging data. The Goldstone terrestrial 
radar proAudes better imaging resolution (several tens of km), though limited 
coverage [18]. 

In situ measurements of venusian surface characteristics were made by Venera 
soft-landers. Venera 9 and 10 each sent back one TV picture (Fig. 5-5). 


5.5.1 Surface Morphology 

Venus has a flat surface compared to the Earth, the Moon, or Mars. Ninety 
percent of the surface lies within a height interval of 3 km. Extremes: the highest 

point, Maxwell Montes, is 11.1 km above the planetary datum; the lowest, a trench 
named Diana Chasma, is 2.0 km below datum. Total relief, 13 km, is two-thirds that 
on Earth (20 km). Venusian highlands are relatively small in area and low in eleva- 
tion. The two major highland regions, Ishtar and Aphrodite Terrae, resemble con- 
tinents on the Earth in area. 
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Analysis of the Pioneer Venus Orbiter radar data (wavelength: 17 cm) shows 
that the largest percentage of Venus' surface has meter-scale rms slopes in the 
range of 1 to 3 deg, much smoother than the Moon or Mars. 


5. 5.1.1 Rolling Plains Province 

The rolling plains, the white areas in Figure 5-3, make up about 65 percent of . 
the surface. They are at elevation 0 to 2 km above datum, with low slopes (1.5 to 
3 deg), and a number of craters. Volcanic features appear to be present. The 
rolling plains may be an ancient surface that represents the original venusian crustal 
differentiate. Venera 8 obtained a "granitic" composition from gamma ray spectrometry 
data [5-7] . 


5, 5.1. 2 Highlands 

The highlands make up less than 8 percent of the imaged surface of Venus, but 
their radar brightness and elevation are striking. They occur in three areas (1) 
Ishtar Terra, centered at latitude 65 deg N, longitude 350 deg; (2) Aphrodite Terra, 
on and south of the equator between lon^tude 60 and 205 deg; and (3) Beta Regio, 
between latitude 15 and 40 deg N and longitude 275 and 296 deg. Aphrodite is 

roughly the size of Africa; Ishtar the size of Australia. Ishtar contains the central 

mountainous area of Maxwell Montes. 

Maxwell Montes is a very bright radar feature, implying very steep debris -laden 
slopes. There is no corresponding gravity anomaly, a condition that implies isostatic 
compensation. Maxwell may be more silicic than Beta, which is probably basaltic. 

Complex terrain is found in the Eastern Ishtar Terra region, with ridges, 
troughs, and relatively less boulders and other debris than on Maxwell's slopes. 
Aphrodite Terra also shows complex topography and a bright radar return, though 
somewhat less pronounced than Maxwell. 

Beta Regio is dominated by the irregular shield-shaped features Theia Mons and 
Rhea Mons, near which Venera 9 and 10 landed. These featiues are thought to be 

volcanic, and Beta may be one of the youngest volcanic regions on the planet. 


5. 5. 1. 3 Lowlands 

The most extensive lowland basin on the planet is Atalanta Planitia, centered at 
latitude 65 deg N, longitude 165 deg. It is about the size of the Gulf of Mexico and 
is about 1.5 km below datian. It appears very smooth, like the lunar maria, and 
shows no apparent cratering. Thus, it may represent a basaltic: lava flow. (Venera 
landers derived a basaltic composition for volcanic rocks near Beta, which lends 
credance to this hypothesis [5-7]). 


5. 5. 1.4 Tectonic Features 

A number of linear features have been observed, associated with Alpha Regio, 
Maxwell Montes, Beta Regio, and regions east of Ishtar and Aphrodite Terrae. 
Classed as ridge-and-trough systems, lineaments and scarps, it is possible to inter- 
pret them as tectonic faults. (The suggested tectonic activity is local or regional, 
not global.) 
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Figure 5-3. Venus topographic provinces. 

Map showing distribution of topographic provinces on Venus: rolling plains (no pattern), 

0(6051.0) to 2 km (6053.0); highlands (hatching), higher than 2 km (6053.0-6062.1); and 
lowlands (dotted pattern), lower than 0 km (6049.0-6051.0). The numbers in ^parentheses are 
planetary radii in kilometers [5-18]. 








5. 5. 1. 5 Volcanic Features 


The resolution available so far does not allow for positive distinction between 
volcanic and impact features. However, the' morphology of regions such as Beta 
Regio and differences between crater densities and preservation in different regions 
indicates volcanic activity and, presumably, associated tectonic activity. 


5. 5.1.6 Craters 

The smallest detectable crater is about 20 km in dimneter. Crater statistics 
are given in Figure 5-4 [5-18]. 



Figure 5-4. Venusian crater densities. 


A tentative conclusion from the crater statistics is that the venusian surface is 
relatively young. However, better data await improved surface imaging resolution, 
either from a planetary orbiter or by using synthetic aperture radar from the Earth 
[5-15]. 


5.5.2 Surface Characteristics 


In situ observations have only been made by the Venera so ft -landers. The 
Venera TV pictures (Fig. 5-5) show a rocky desert-hke terrain. Weathering appears 
to be extreme, presumably due to high- temperature sintering, corrosion, and aeolian 
processes. Wind velocitie^s of up to several meters per second are plenty high enough 
to move surface material [5-221, 
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Figure 5-5, Phototelepanoramas of the Venus surface from Venera 9 
(top) and Venera 10 (bottom) [5-7], 


The Pioneer Venus Day Probe raised a dust cloud upon impact — it appeared to 
be composed of particles in the 10- to 100- mi cron diameter size range [5-23], 


5.5.3 Chemical Composition — Mineralogy 

Gamma-ray spectrometers on Venera 8, 9, and 10 landers measured the surface 
abundances of Uranium, Thorium, and Potassium (Table 5-4), The inference is that 
rocks at the Venera 8 landing site resemble terrestrial granites, while at Venera 9 
and 10 sites basalts predominate. 


TABLE 5-4. SURFACE ABUNDANCES OF U, Th, AND K [5-7] 



Uranium 


Thorium 


Potassium 

Venera 8 

2.2 ± 0. 7 ppm 

6.5 

±2.0 ppm 

4. 9 

±1.3 percent 

Venera 9 

0.5 ± 0. 2 ppm 

4.0 

± 0. 25 ppm 

0.9 

±0.2 percent 

Venera 10 

0.7 ± 0.2 ppm 

1.1 

±0.2 ppm 

0. 3 

±0.1 percent 

Typical Terrestrial Rocks; 





Granite 

9.04 


21. 9 


3.24 

Basalt 

0.86 


2.1 


0.76 


5.5.4 Physical Properties 

The density of material at the Venera 10 site was determined to be about 2.7 

_ 3 

g cm for the top few decimeters. This value is about what would be expected for 
basaltic rocks and is not in conflict with estimates previously made from radar 
reflectivities [5-7]. 
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5. 6 ATMOSPHERE 


Observations of Venus from the Earth show only the top of a thick cloud, layer, 
whose nature was one of the key scientific questions leading to exploration of the 
planet. The extremely high- atmospheric temperature and pressure at the surface 
(740®K, 92 atm) and the fact that CO 2 is the major constituent, were, known prior to 

Pioneer Venus. Terrestrial observations had led to the identification of cloud particles 
as micron- si zed sulfuric acid droplets [5-24]. However, the Pioneer Venus mission pro- 
duced many surprises. One was the discovery of solid crystals' in the clouds in addi- 
tion to sulfuric acid droplets. The structure of the lower sitmbMp here was studied by 
instruments carried on the probes, as well as indirectly from the Orbiter. The struc- 
ture of the upper atmosphere was derived from the atmospheric drag on the Orbiter. 

This led to one of the most important and perplexing results to come out of the 

Pioneer Venus mission: the very low temperature, near 100°K, determined for the 

upper atmosphere on the nightside (Contrast the nightside exospheric temperature on 
the Earth of greater than 800°K.) 

The three major altitude regions of the atmosphere of Venus will be discussed 
separately below. They are (1) the Lower Atmosphere, below 65 km, which includes 
the cloud layer from 25 to 65 km; (2) the Middle Atmosphere from ‘65 to 120 km; and 

(3) the Upper Atmosphere, above 120 km. The information reporj^ed was derived 

almost entirely from the Pioneer Venus mission. A recent general review can be 
found in Reference 10. 


5.6.1 Lower Atmosphere 

State properties of the Lower Atmosphere were measured by instruments carried 
on the four Probes and by Venera probes and landers. The information obtained is 
given in the following paragraphs. 

Temperature and pressure data from the four probes were quite similar. Figure 
5-6 [5-25] shows contrast between the probes (as compared to data from the Sounder 
Probe). 

ALTITUCC (km) 



Figure 5-6. Temperature soundings of the Venus atmosphere. 
The altitude scale is for the SoTinder Probe. 
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5. 6, 1.1 Models of the Lower Atmosphere 


Table 5-5 txsntains thermodynamic pix>perties of the lower atmosphere of Venus 
for latitudes <45 deg. 

TABLE 5-5. THERMODYNAMIC PROPERTIES OF THE LOWER 
ATMOSPHERE (<45 deg) [5-9] 


Altitude 

(km)'*’ 

Average model 


Maximal model 


Minimal model 


T(K) 

P(alni) p(kgm'*) 

riK) 

Plaimi /Mkgm ’) 

r(K) 

Platml plkgm’i 

-6 

780 



783 



777 



-4 

765 

115.3 

79.9 

768 

117.3 

81.3 

762 

1 12.9 

78 .’ 

-2 

750 

102.0 

72.1 

753 

104.2 

73.4 

747 

99.8 

■'O.K 

0 

735 

90.0 

64.8 

738 

92.0 

66.1 

732 

89,0 

63.'' 

2 

720 

79.2 

58.3 

723 

81.0 

59.4 

717 

■'7,4 

57.3 

4 

705 

69.6 

52.3 

708 

71.2 

33.3 

702 

67.9 

51.3 

6 

680 

60.9 

46.8 

692 

62.3 

47.7 

686 

39.4 

45 6 

8 

. 674 

33.2 

41.8 

678 

.34.4 

42 ft 

670 

’ 1.8 

41.0 

10 

658 

4ft,2 

.37.3 

661 

47.4 

38.0 

653 

43.1 

.36.3 

12 

642 

40.1 

33.1 

645 

41.1 

33.S 

639 

39 1 

32.4 

14 

626 

234.7 

29.4 

629 

35.6 

.30.0 

623 

' 

28 8 “ 

Ifi ' 

' 609 

29.9 

26.0 

612 

30.6 

26.3 

606 


25.4 

IK 

59.3 

25.6 

22.9 

596 

26.3 

23.4 

390 

24,9 

22.4 

20 

576 

21.9 

20.1 

579 

22-5 

20.6 

373 

21.1 

19.6 

TO 

560 

IK. A 

17.6 

.56.3 

19 1 

18.0 

337 

18.0 

17.2 

24 

543 

13 7 

15.4 

.346 

16 2 

13 7 

’40 

13.2 

13,0 

2 b 

52ft 

1 1.3 

13.4 

529 

1 < ft 

13.7 

323 

12.3 

13 0 

28 

509 

111 

11. ft 

512 

11.4 

1 1.8 

306 

10.7 

11.2 

?0 

492 

9.23 

9.95 

495 

9. 84 

10.2 

489 

8.92 

9 o’ 

32 

47ft 

7.63 

8.51 

479 

7.90 

8,73 

473 

7.36 

8.24 

34 

4ft0 

ft.27 

7.23 

463 

6.49 

7.44 

457 

6.04 

•’ o| 

3ft 

444 

3.12 

ft.l 1 

447 

3 31 

6.29 

441 

- 4.92 

5 d2 

38 

428 

4.15 

5.13 

431 

4.31 

5.?'* 

423 


4 Dh 

40 

413 

3.. 34 

4.28 

416 

.3.47 

4.42 

J ID 

3,20 

4.14 

4 2 

4illl 

2.67 

3.33 

403 

2.78 

3.65 

5*l7 

2.55 

3.40 

44 

387 

2.11 

2. .89 

390 

2.21 

2.99 

i 

2.02 

2. ■’9 

4ft 

374 

i.ftft 

2.3ft 

378 

1.74 

2.43 


1,39 

•>17 

4H 

3ft 1 

i.30 

l.ftl 

3ft8 

1 .36 

1.96 

^’3 

1.23 

1.84 

^0 

347 

l.uo 

1 .I.! 

^S7 

1.06 

1.57 

t ID 

■' 43R-1 

1.46 

^2 

328 

7.ft4E-l 

1.24 

338 

8.I1E-1 

1.28 

i") 

■ 13 E-1 

1.1 ft 

'4 

3(I)S 

8.7.1E-1 

9.87 E-1 

319 

6.I4E-1 

1.03 

'Dll 

23E-I 

9.43E-1 

^6 

2 . 8 ft 

4.22E-1 

7.73E-1 

um 

4..58E-1 

8.11 E-1 

' D 

' 90E-1 

7.36E-I 

38 

2^4 

3.05E-1 

3.91 E-1 


3..36E-1 

ft.2.3E-l 


79E-1 

’ 3KE-1 

ftO 

2 ft 1 

2.17E-I 

4.4IE-I 

*>77 

2 42E-1 

4,70F-I, 


. 9ft E-1 

4.I2F-I 

m 2 


I.32E-I 

.3. 2 3 E-1 

2(>2 

1 .72E- 1 

3.30I--I 

:-:o 

1 3’E-l 

2.9SE-I 

b4 

245 

1.05E-1 

2.29E-I 

2^7 

I.21E-I 

2.30R-I 

' 'ft 

ft 21 E -2 

2.08E-1 

flft 

241 

7.24E-2 

1.38E-1 

252 

8.49E-2 

1.79E-1 


6.24E-2 

1.43E-1 


2.3ft 

4.95E-2 

1.1 IE-1 

248 

3.90E-2 

I.23E-1 

■' T 7 

4.20E-2 

D 'up . 2 

'() 

232 

3.3ftE-2 

7 t»8E-2 

243 

4.08E-2 

8.91R-2 

' 3 ■> 

2.80E-2 


*'5 

220 

1 .20E-2 

2.'OE-2 

2 ii» 

1.39E-2 

.3.32E-2 

210 

9 93E-3 

2.41R-2 

Stt 

-OS 

4.C)7E-3 

I.II3E-2 

2! 8 

3.44R-.3 

1 32R-2 

2D() 

.■'.09E-3 

8.1DH-5 


192 

1 .27E-3 

3.30E-3 

202 

l.,SOR-3 

4.7IR-2 

‘ 


2 .h''h-.l 

'Ml 

1 « 


1.071-.-3 

s o 2 

3 33R-4 

l,33l'-i 

1 ‘ "ft 

2,’,sR-4 

■■.84F-4 

•t' 

! s 2 

i.02E-4 

■Vo U--4 

i'»ii 

2.0 H-: -4 

4:.=i3l-4 

' ' 1 ■ft 

o I'l H-3 

2.! II--4 

1 1 M 1 

1 73 

2.72E-5 

8.23E-3 

! 0 (l 

4.S7R-4 

1.33E-4 

! dD 

l..’8E-5 

5,2ftF.-3 


Above reference le\ci R =60>’ km. » 
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5. 6. 1= 2 Vertical Stability 

Estimates of atmospheric vertical stability are difficiilt since they depend upon 
the z-derivative of the temperature. The atmosphere was found to be stably strati- 
fied for all altitiudes down to about 28 km, except for a nearly adiabatic region in the 
clouds at 50 to 55 km. The lapse rate is nearly adiabatic again from 20 to 28 km, 
and the Pioneer Venus probes measured stable stratification again from 20 km down to 
10 km. Their estimates do not extend below this level 15-25], Venera 10 data confirm 
these findings and extend stability estimates down to the surface, showing a neutral 
stability layer below 5 km. Detailed information is given in Reference 5-25. 


5. 6.1.3 Chemical Composition 

The mean moleciilar weight for all altitudes below 80 km is about 43.44 [5-26]. 

Gas chromatographic data on atmospheric composition were obtained at 51.6, 
41.7, and 21,6 km altitudes. The data for 21.6 km are given in Table 5-6 [5-27]. 


TABLE 5-6. COMPOSITION OF ATMOSPHERE AT 21.6 km [5-27] 


Gas 

Volume Mixing Ratio (Percent) 

Comments 

CO 2 

96.4 ± 1.0 


^2 

3.41 ± 0.01 


H 2 O 

0.135 

Upper limit ; the value is 
uncertain. H 2 SO^ in the 



clouds is a sink for H 2 O [5-27]. 

^2 

16. 0 ± 7.4 X 10 "^ 

This is the 41-km value 
(upper Umit) [5-9]. 

Ar 

67.2 + 2.3 X 10 "^ 

Note isotope ratio in following 
table . 

CO 

19. 9 ± 3.12 X 10'^ 

Mixing ratio falls off rapidly 
with altitude. 

Ne 

4. 31 ± 4 X 10 "'^ 


SO 2 

185 (+350 - 155) X 10'^ 

Only sulfur compound actually 
detected . 


Using the mixing ratios measured directly, together with the assumption of 
equilibrium conditions , model curves were calculated for two values so f water vapor 
as shown in Figure 5-7 [5-27]. 

Table 5-7 gives valdes for atmos|)heric constituents based on Venera and earlier 
spectroscopic data [5-28,5-29], 
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(a) (b) 

Figure 5-7. Models of the composition of the Venus atmosphere as a function 
of altitude, (a) Water = 0.135 percent (one observation by Pioneer 
Venus) and (b) water = 50 ppm (reported Venera value). 


TABLE 5-7. ATMOSPHERIC COMPOSITION BASED ON VENERA 

AND EARLIER DATA [5-28] 


C>(is 

o 

o 

N5 

^'2 

SO 2 

H2O 

S.2 

HCl 

HF 

CO 

cos 

H2S 

»2 

Mixing' 

Rutio 

0. 9G 

3.4 10'“ 

1.3 ■ lO"* 

2 . lO-* 

1 • 10"^ 

1 ■ 10« 

00 

0 

1.5 10'^ 

2 • 10‘® 

3 ■ lO’’' 



The Neutral Mass Spectrometer on the Pioneer Venus Sounder Probe provided 
independent data on mixing ratios for minor constituents. Agreement was found with 
the Gas Chromatograph data except for possible discrepancies for water vapor. Upper 
limits were put on the carbonyl sulfide (COS) and sulfur dioxide (SOg) mixing ratios, 
data of importance in atmospheric chemistry [5-30]. 


Mixing Ratio, ppm 


SO 2 

below 

55 km 
24 km 

< 10 
<300 

cos 

above 

below 

24 km 
24 km 

< 3 
<500 


5. 6.1.4 Isotopic Composition [5-30] 

A comparison of the isotopic ratios of some gases in the Venus atmosphere to 
the Earth's atmosphere is presented in Table 5-8. 
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TABLE 5-8. ISOTOPIC RATIOS OF SOME SUBSTANCES FOUND 
IN THE VENUS ATMOSPHERE 


Gas 

Venus Atmosphere 
Isotopic Ratio 

Earth Atmosphere 
Isotopic Ratio 

’He/"He 

<3 X 10~* 

1.4 X 10-‘ 

“Ne/^®Ne 

0.07 ± 0.02 

0,09' 

^“Ne/’‘Ar 

0.3 ± 0.2 

0.58 

”Ar/“Ar 

0.18 ± 0.02 

0.187 

““Ar/^Ar 

1.03 ±0.04 

296 


< 1.19 X 10-^ 

1.11 X 10-^ 

I 8 Q/I 60 

2.0 ± 0.1 X 10"^ 

2.04 X 10"’ 


22 20 

The neon isotopic ratio, Ne/ Ne, is about 25 percent lower than the value 
for the Earth, but is close to that of the solar wind. 

40 

The major Venus-Earth difference is for the argon isotopes, Ar, which is 
40 36 

radiogenic (from K decay), and Ar, which is non- radiogenic. Allowing for the 
differences in total atmospheric mass, the quantity of radiogenic argon in the Venus 
atmosphere is about one-third that of the Earth, while the Venus atmosphere contains 
about twice as much non- radiogenic argon as the Earth. 


5. 6.1.5 Winds in the Lower Atmosphere 

The zonal winds in the lower atmosphere blow from east toward the west in aU 
locations, that is, in the direction the planet rotates. Tj^pical wind speeds are listed 
in Table 5-9 [5-25]. 


TABLE 5-9. TYPICAL WIND SPEED PROFILE 


Altitude 

Wind Speed 

Surface 

0. 3 to 1.0 m sec ^ 


-1 

7 km 

1 m sec 


-1 

28 km 

30 m sec 


-1 

65 km 

100 m sec 


In the altitude range 28 to 65 km, zonal winds increase from the surface to 

100 m sec ^ near cloud top, which is the 4-to-5 day "superrotation" observed from 
the Earth. The mechanism responsible for driving the zonal circulation is not yet 
wen understood [5-1]. Figures 5-8 through 5-10 show the zonal winds: speed, mass- 

flow, and energy density from probe data [5-25]. Figure 5-11 shows sketches 
of possible meridional flows and eddies. In addition to zonal flow, 
Venus has meridional winds whose speed is an order of magnitude less^ 

At cloud top, the meridional wind speeds are from a few to 10 m sec~ 
in the poleward direction. This corresponds to the top of a Hadley cell 
circulation system, with equatorward flow at about 50- to 55-km altitude 
[5-25] . 
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Figure 5-8. Zonal winds. 


Figure 5-9. Zonal mass flow. 



KINETIC ENERGY DENSITY 
(103 J/m3) 

Figxore 5-10. Zonal kinetic energy density. 


5. 6. 1. 6 Clouds 

The cloud layer on Venus continuously obscures the surface. It is featureless 
except in the ultraviolet region of the spectrum, where low- contrast dark markings, 
often with a characteristic horizontal- Y shape, frequently appear and persist for hours 
or days. (Review paragraph 5. 6. 2.2 for a disctission of UV absorbers .) . These 
recurrent features are shown diagrammatically in Figure 5-12, along with the descrip- 
tive names associated with their structure in the literature [5-31], 
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Figiire 5-11. Meridional flows and eddies [5-25]. 



Figijre 5-12. Schematic diagram defining the basic types of cloud features 
observed in Venus ultraviolet images. The two views depicted here 
typically occur 2 days apart and represent the maximum and minimum 
tilt configurations , respectively. Tick marks on each circle indicate 
20 and 50 deg of latitude in each hemisphere. 


The clouds play an important role in the radiation budget of the planet (para- 
graph 5.6.4). Prior to the Pioneer Venus mission, careful analysis of ground-based 
observations led to the remarkable conclusion that the clouds are made up of single - 
mode droplets of stilfuric acid with radius 1.05 microns [5-24]. The Sounder Pix>be 
instruments confirmed the presence of such droplets, but also found populations of 
smaller and larger particles. Data was collected both from the Orbiter and by direct 
measurements by the particle size spectrometer on the Sounder Probe during entry. 

The great opacity of the cloud deck, which has an optical depth in. excess of 
30, results mainly from its very great thickness. Visibility within the clouds is 
nowhere less than several kilometers [5-32]. 
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Haze exists throughout the atmosphere except in the lowest 25 km, which is 
remarkably clear. Haze extends upwards to altitudes of at least 90 km [5-32]. 

The particles that form the clouds may be divided into three clases: Mode 1, 

submicron aerosol "haze" particles; Mode 2, micron-sized H2SO^ droplets; and Mode 3, 

larger crystals. The Mode 3 crystals are known to be larger solid particles. Their 
form and composition is uncertain, though the current best guess is that they are a 
chloride [5-28,5-33] and are plate-like in crystalline form. They may account for the 
major fraction of the total columnar mass loading of cloud particles (approximately 

1. 8 X 10 ^ g cm ^) . 

Earth and Venus cloud system properties are compared in Table 5-10. 

TABLE 5-10. COMPARATIVE VENUS AND EARTH CLOUD 
SYSTEM PROPERTIES [5-33] 


Property 

Earth 

Venus 

Percent coverage 

40 

100 

Average optical depth 

5-7 

25^ 

Maximum optical depth 

300-400 

40 

Composition 

solid and liquid HjO 

H2SO4 droplets, crystals, 
plus contaminants 

Number density 

100-1000 cm“’ (liquid) 

50-300 cm"’ (liquid) 

0. 1-50 cm"’ (ice) 

10-50 cm*^ (crystals) 

Average mass loading 

0.3-0.5 g m"’ 

0.01-0.02 gm"’ 

(Mass Density) 


0.1-0.2gm-^ 

Maximum mass loading 

10-20 g m"* 

(Mass Density) 

Distribution function 

normal-log normal 
bimodal (ice) 

Multimodal 

Typical background aerosol 

1 cm“^ 

100-200 cm-’ 

at cloud base 0.5 microns 

Condensation process 

homomolecular 

heteromolecular 

Average pfecipitable mass 

0.03-0.05 mm 

0. 1-0.2 mm 

Mean scattering size (diameter) 

10;im 

2-4 ;im 

(In the visible) 

Mean mass size (diameter) 

30 jam 

lOfim 

Dominant optical cloud form 

stratiform 

stratiform 

Dominant mass cloud form 

cumulus 

stratiform 

Potential latent instability 

high 

low 

Temporal variability 

high 

slight 

Dominant cloud atmosphere 

latent heat 

radiation 

beat exchange process 


Table 5-11 summarizes the properties of the various regions [5-33]. For radia- 
tive properties review paragraph 5.6.4. Over the poles the cloud tops are 6 to 9 km 
lower [ 5- 34] . 


5. 6. 1. 7 Lightning 

Whistler mode electromagnetic noise bursts were detected by the plasma wave 
instrument aboard the Pioneer Venus Orbiter [5-35] . Presumably these were produced 
by lightning discharges in the lower atmosphere , though no direct confirmation has 
been made. Several of the Venera spacecraft also detected electromagnetic noise 
signals that could be attributed to lightning [5-36]. An optical search for lightning 
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TABLE 5-11. SUMMARY OF VENUS CLOUD PROPERTIES 


Region 

Altitude, 

km 

Temperature, 

“K 

Optical 

ciepth 

T 

Average 
Number 
Density 
N, cm"’* 

Mean 

Diameter, 

pm 

Compositionf 

Upper hazet 

70-90 

225"- 190“ 

0 . 2 - 1.0 

500 

0.4 

H 2 SO 4 + contaminants 

Upper cloud 

56.5-70 

286"-225" 

6 . 0 - 8.0 

(1) 1500 

(2) 50 

Btmodal 
0.4 and 2.0 

H 2 SO 4 + contaminants 

Middle cloud 

50.5-56.5 

345"-286“ 

8-10 

(1) 300 

(2) 50 

(3) 10 

Trimodal 
0.3. 2.5. ' 0 

H.SO 4 + crystals 

Lower cloud 

47.5-50.5 

367"-345" 

6-12 

( 1 ) 1200 

(2) 50 

(3) 50 

Trimodal 
0.4, 2.0, 8.0 

H 2 SO 4 + crystals 

Lower haze 

31-t7.5 

482“-367" 

0 . 1 - 0.2 

( 2)20 

0.2 

H 2 SO 4 + contaminants 

Precloud layers 

46 and 47.5 

378" and 367" 

0.05 and 0. 1 

50 and 150 

Bimodal 
0.3 and 2.0 

H 2 SO 4 + contaminants 


flashes showed only that lightning activity on Venus is not greater than on the Earth. 

2 

An upper limit was set of about 30 flashes per km per year [5-37]. This does not 
appear to be enough activity to make lightning an important factor in the production 
of trace constituents in the atmosphere. 

The origin of lightning in the Venus lower atmosphere poses an interesting 
problem. Even though the mechanisms responsible for production of hghthing in 
Earth clouds are not completely understood, it is certain that microphysical processes 
involving water-ice play the key role. (The only important exception is for lightning 
associated with volcanic activity. ) On Venus , the Mode 3 particles , since they are 
crystalline, may be responsible for charge separation [5-32]. 


5. 6. 2 The Middle Atmosphere 

Above the cloud top level at 65 km to the homopause at approximately 130 km 
is the Middle Atmosphere. Temperature decreases to about 170°K between 90- and 
100-km altitude, then increases again to around 190®K at 110 km (0.002 mb). The 
altitude variation of temperattme with altitude (pressure) is strikingly different from 
the Earth as measured by infrared remote sensing from the Orbiter [5-38]. Thermo- 
dynamic properties of the middle atmosphere are listed in Table 5-12. 

Tentative chemical composition results are available from a recent Soviet model 
calculation (Fig. 5-14) [5-28]. 


5. 6. 2,1 Winds in the Middle Atmosphere 

No direct measurements of wind speeds were made above 60 km (cloud tops) , 

but strong zonal wind flows are expected. Models predict a maximum zonal wind 

- 1 -1 
speed of about 150 m sec at 70 km, decreasing to about 133 m sec at 95 km, _ and 

increasing again to superrotation speeds at 150 km [5-25]. 
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Figure 5-13. Profiles of temperature versus pressure at 30 deg N latitude 15-25,5-39]. 

TABLE 5-12. THERMODYNAMIC PROPERTIES OF THE VENUS MIDDLE 

ATMOSPHERE (SOUNDER PROBE) 


Altitude, 

km 

p, ban 

T, K 

p, kg/m’ 

116 

0.5199E-06 

190.2 

0.1405E-05 

114 

0.8255E - 06 

190.2 

0.2236E-05 

112 

0.1318E-05 

190.2 

0.3577E - 05 

110 

0.2075E - 05 

194.4 

0.5520E - 05 

108 

0.3273E-05 

196.3 

0.8630E - 05 

106 

0.5171E-05 

190.1 

0.1409E-04 

104 

0.8338E - 05 

183.0 

0.2363E - 04 

102 

0.1362E-04 

180.7 

0.391 IE -04 

100 

0.2252E-04 

174.1 

0.6718E-04 

98 

0.3773E - 04 

173.3 

0.1132E-03 

96 

0.6350E-04 

168.9 

0.1957E-03 

94 

0.1078E-03 

168.9 

0.3326E-03 

92 

0.1848E-03 

167.1 

0.5768E - 03 

90 

0.3168E-03 

165.4 

0.9999E-03 

88 

0.5469E - 03 

165.2 

0.1728E-02 

86 

0.9402E - 03 

163.8 

0.2997E - 02 

84 

0.1634E-02 

164.7 

0.5184E-02 

32 

0.2757E - 02 

181.1 

0.7951E-02 

80 

0.4449E - 02 

195.8 

0.1188E-01 

79 

0.5577E - 02 

204.1 

0.1428E-01 

78 

0.6926E-02 

212.4 

0.1704E-01 

77 

0.8553E - 02 

215.1 

0.2077E-01 

76 

0.1054E-01 

216.9 

0.2540E - 01 

75 

0.1296E-01 

221.1 

0.3064E-01 

74 

0.1587E-01 

227.3 

0.3647E - 01 

73 

0.1934E-01 

227.8 

0.4435E-01 

72 

0.2349E - 01 

233.5 

0.5256E - 01 

71 

0.2850E - 01 

234.9 

0.6339E-01 

70 

0.3440E-01 

243.5 

0.7381E-01 

69 

0.4166E-01 

230.5 

0.9444E - 01 
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Figure 5-14. Composition of Venus atmosphere between 
50 and 100 km (tentative). 


5. 6. 2.2 Ultraviolet Absorbers 

The dark patterns observed at ultraviolet wavelengths in the cloud tops have 
been shown by carefiil analysis to require the presence of at least two different UV 
absorbers in the middle atmosphere. The absorption from 2000 to 3200 A can be 
accounted for by sulfur dioxide at around 75 km within the haze layer above the main 
cloud deck. At longer wavelengths the observed dark patterns cannot be explained 
by SO 2 , since its absorption drops to too low a value. Therefore, a second broad- 
band absorber is reqioired in the first few optical depths of the upper cloud level 
[5-27]. Elemental sulfur was once suggested, but it has since been shown to be 
insufficient [ 5- 40] . Chlorine gas or chlorine compounds are candidates , possibly FeCl^ 

dissolved in the concentrated iH 2 SO^ [5-28]. Interest in UV absorbers is spurred by 

recognition that they play an important role in the variation in local energy deposition 
in the atmosphere. 


5. 6. 3 The Upper Atmosphere 

The temperature of the dayside Upper Atmosphere had been well established by 
remote sensing prior to the Pioneer Venus mission [5-12]. The day-night temperature 
structure is compared to the Earth in Figure 5-15. Venus' dayside Upper Atmosphere 
is not unlike the Earth's thermosphere, although temperatures are considerably lower. 
However, the nightside upper atmosphere temperature structure is totally unfamiliar 
and has not yet been explained. (It has been termed a "cryosphere" [5-25].) 


5. 6. 3. 1 Upper Atmosphere Temperature 

Asymptotic temperatures for great altitudes become: 

Dayside: T ^ 285°K 

Nightside: T ^ 110°K 

Global Average: 'v 199°K 
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Figure 5-15. Comparison of the vertical thermal structures of the atmospheres 

of Venus and the Earth [5-25,5-39]. 

5. 6. 3. 2 Upper Atmosphere Density and Composition 

Representative densities for the dayside and night side are [5-41]: 

- 3 - 3 

Dayside g cm Nightside g cm 


150 

km 

about 

8 

X 10-13 

24 

X less 

170 

km 

about 

7 

X lO"!'^ 

82 

X less 


The variation contrasts with a diurnal variation of only 1 percent on the Earth. 

Figure 5-16 shows the diurnal variation in atmospheric density at 160-km altitude 
as derived from Orbiter drag data [5-41]. The temperature drops by 150®K going 
across the terminators within only 30 deg of lon^tude. 



Figure 5-16. Diurnal variation of upper atmospheric density. 
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The temperat Tires shown in Figure 5-17 are also for 160 km, derived from 
Orbiter drag data. The cxirve labeled represents the assumed daily variation of 

temperatures at the 115-km level [5-41]. 



Figure 5-17. Diurnal variation of upper atmospheric temperature. 


MASS DENSITY /> (g/cm^) 

IQ-17 10“'* 10"'’ 10“’’ 10“” 



Figupe 5-18. Tentative model nimber and total mass 
densities versus altitude [5-42]. 
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According to the model resxilts, oxygen atoms take over from CO2 as the domi- 
nant constituent for all altitudes above 155 km on the dayside and at a somewhat 
lower altitude on the hightside. For all constituents except helium and hydrogen, 
dayside density exceeds nightside. 




Figure 5-19. Diurnal variaton of He, N, O, N^, CO, and CO„ 

at 157-km altitude 15-42]. 


5. 6. 3. 3 Winds in the Upper Atmosphere 

Model calculations suggest wind speeds increasing from 100 m sec ^ at 125 km to 

200 m sec ^ above 150 km [5-43]. Superrotation above 150 km is superposed on a 
general subsolar to antisolar flow and may only occvir at low altitudes [5-25]. 

5. 6. 4 Radiative Properties of the Venus Atmosphere 

Effective radiating temperature of the planet : 231 ± 7°K 

Bolometric Bond Albedo: 0.75 ± 0.02 

- conformed to flux (in or out) of 162 W m ^ [5-9]. 
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5. 6. 4. 1 The "Greenhouse" 

A research problem of fundamental importance is to understand how the surface 
of the planet is maintained at a temperature near 740°K. The only viable explanation 
appears to be that the atmosphere as a whole acts as a "greenhouse" to trap infrared 
radiation from the surface, while admitting enough solar radiation to act as the primary 
heating source. Pioneer Venus data may not be fully consistent with the greenhouse 
hypothesis [5-44] , but one-dimensional- models have been fairly successful. The role of 
various absorbers in contributing to a recent greenhouse model is illustrated in Table 
5-13 [5-45]. Ah interesting point is that in this model the clouds are not the major 
factor — they take third place behind CO 2 and water vapor. Any one -dimensional 

model must be incomplete, however, since horizontal motions must play the role of 
redistributing the nonuniform heating into the observed uniform temperature. 

TABLE 5-13. DECREASE IN SURFACE TEMPERATURE ACCOMPANYING 
THE REMOVAL OF VARIOUS SOURCES OF INFRARED 

OPACITY [5-45] 


Source Deleted 

Change in Surface Temperature, K 

HCI 

0.1 

CO 

13 

SOi 

52 

Clouds 

113 

H 2 O 

218 

CO 2 

463 


It is important to note that the exact role of water vapor in the infrared opacity 
will remain uncertain vintil discrepancies in the measurements of its concentration in 
the atmosphere are resolved. Another factor that undoubtedly contributes to the heat 
balance in the atmosphere is atmospheric circulation, but its exact role is uncertain 
at present. Reference 5-25 discusses questions of heat transport by meridional flow 
and eddy diffusion. 


5. 6. 4. 2 Optical Properties of the Venus Cloud System 

The cloud layer has a profound influence on the penetration of solar energy and 
therefore on. the input; side of the radiative balance of the atmosphere. Figure 5-20 
[5-32] shows the results of measurements made on the Sounder Probe. 


5. 7 IONOSPHERE 

Extensive ionospheric data were obtained during the Pioneer Venus mission. 
The interactions between the ionosphere and the solar wind are complex and difficult 
to summarize. For details, re-view published acco-unts such as References 5-46 
through 5-49. 


5.7,1 lonsopheric Ions - Temperature and Total Densities 

Ion densities except for H^ decrease by over an order of magnitude from the 
dayside to the nightside, while ion temperature, constant over the entire dayside, 
increases to peak at the antisolar point. Figures 5-21 and 5-22 summarize data on ion 
densities [5-48]. 
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Figure 5-2i. Ion density profiles. Figure 5-* 22. Ion temperature profiles. 


5.7.2 Ions — Composition 

The Venus ionosphere is mostly O above 200 km and mostly below this, 

down to 160 km. All ions other than and He"*^ are more abundant on the dayside, 
while the two lightest ions show asymmetric nightside bulges. Typical ion concentra- 
tions at 200 km at the Venus ionosphere are shown in Table 5-14. 

TABLE 5-14. TYPICAL lONSOPHERIC ION DENSITIES 
AT 200 km (SUBSOLAR) [5-49] 


Ion 

_ 3 

Density Number cm 

0 + 

10^ 

O2" 

4 

6 X 10^ 


20 

NO^ 

10^ 

CO2" 

2 X 10^ 

He"^ 

100 

CO"^ plus 

10^ 


8 X 10^ 

n'^ 

10^ 
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5.7.3 Ionospheric Electrons 


Electron temperatures vary by less than a factor of two from subsolar point to 
antisolar point (Fig. 5-23). Below, 180 km the electron temperature rises from dayside 
to nightside, while above this level the converse is true. Electron densities drop by 
two orders of magnitude from subsolar to antisolar point. Table 5-15 [5-48] summarizes 
the results of model calculations based upon Pioneer Venus data. The cut-off altitude 
in the table was taken as a representative ionopause location, though the actual iono- 
pause was found to be extremely variable (paragraph 5.7.4). 



FtFCTRON TFMPFRATURF (KI 


Figure 5-23. Median electron temperatures from 30 deg intervals 
in SZA showing the similarity of median electron temperature 
profiles at different solar zenith angles. 


5.7.4 Ionospheric Variability 

The extent of the ionosphere is highly dependent upon the pressure exerted by 
the solar wind. On the dayside the height of the ionopause may vary between 290 km 

- 8 “2 

at high solar wind pressure (above 4 x 10 dyn cm ) up to 100 km or more for low 

- 9 -2 

pressures (below 5 x 10 dyn cm ) [5-46]. The nightside ionopause height is even 
more variable, with low values below 200 km and high values in excess of 3500 km. 

On the nightside there is evidence of extreme spatial variation in plasma density and 
temperature. Mean ionopause heights are: 
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w TABLE 5-15. COMPUTED MODEL VALUES FOR ELECTRON TEMPERATURE AND ELECTRON NUMBER 

w DENSITY (T^ AND log^.^ (N.)) [5-48] 


Solar Zenith Angle, 


h, km 

0 

IS 

30 

45 

60 

75 

90 

105 

120 

135 

ISO 

165 

180 

160 

1136 

1157 

1210 

1262 

1279 

1249 

1197 

1163 

1180 

1259 

1383 

1502 

1552 


5.50 

5.48 

5.46 

5.42 

5.40 

5.22 

4.88 

4.49 

4.15 

3.92 

3.8! 

3.78 , 

3.78 

no 

1514 

1530 

1567 

1598 

1598 

1560 

1505 

1467 

1472 

1528 

1618 

1703 

1739 


5.30 

5.30 

5.29 

5.28 

5.25 

5.07 

4.73 

4.35 

4.01 

3.78 

3.67 

3.63 

3.62 

180 

1893 

1900 

1915 

1921 

1902 

1856 

1801 

1760 

1753 

1783 

1838 

1890 

1912 


5.18 

5.18 

5.17 

5.16 

5 13 

4.94 

4.61 

4.24 

3.91 

3.68 

3.56 

3.51 

3.50 

190 

2255 

2253 

2243 

2221 

2181 

2129 

2076 

2034 

2015 

2020 

2041 

2062 

2071 


5.10 

5.10 

5.09 

5.06 

5.01 

4.84 

4.51 

4.14 

3.82 

3.59 

3.47 

3.42 

3.41 

200 

2589 

2577 

2542 

2491 

2433 

2375 

2324 

2283 

2254 

2235 

2225 

2220 

2218 


5.05 

5.04 

5.02 

4.98 

4.91 

4.75 

4.43 

4.07 

3.75 

3.52 

3.40 

3.35 

3.34 

225 

3275 

3240 

3150 

3040 

2941 

2871 

2827 

2790 

2743 

2679 

2609 

2555 

2534 


4.96 

4.94 

4.89 

4.81 

4.71 

4.57 

4.26 

3.92 

3.62 

3.40 

3.28 

3.23 

3.22 

250 

3757 

3708 

3581 

3431 

3305 

3226 

3184 

3149 

3092 

3004 

2902 

28^ 

2789 


4.88 

4.85 

4.78 

4.68 

4.56 

4.44 

4.13 

3.82 

3.53 

3.33 

3.21 

3.17 

3.16 

275 

4083 

4027 

3883 

3712 

3570 

3482 

3436 

3400 

3337 

3239 

3124 

3032 

2997 


4.78 

4.75 

4.68 

4.57 

4.45 

4.35 

4.06 

3.74 

3.47 

3.28 

3.17 

3.14 

3.13 

300 

4304 

4246 

4099 

3923 

3773 

3675 

3621 

3576 

3509 

3409 

3295 

3204 

3169 


4.67 

4.64 

4.58 

4.48 

4.36 

437 

3.99 

3.69 

3.43 

3.24 

3.15 

3.12 

3.11 

350 


4516 

4391 

4232 

4083 

3967 

3881 

3808 

3726 

3631 

3535 

3463 

3436 



4.44 

4.40 

433 

4.23 

4.15 

3.89 

3.61 

3.37 

3.20 

3.12 

3.10 

3.10 

400 




4479 

4346 

4211 

4085 

3969 

3864 

377! 

3697 

3648 

3631 





4.23 

4.14 

4.06 

3.82 

3.55 

3.32 

3.17 

3.11 

3.09 

3.10 

450 





4598 

4450 

4281 

4115 

3976 

3876 

3817 

3788 

3779 






4.05 

3.99 

3.77 

3.51 

3.29 

3.15 

3.09 

3.09 

3.09 

500 




I 

4846 

4695 

4484 

4264 

4084 

3968 

3913 

3896 

3894 






3.96 

3.94 

3.72 

3.47 

3.26 

3.13 

3.08 

3.08 

3.08 

6(K) 






5190 

4920 

4598 

4323 

4146 

4066 









3.87 

3.64 

3.40 

3.20 

3.08 

3.04 



700 






5634 

5368 

4972 

4597 

4332 










3.85 

3.58 

3.33 

3.14 

3.03 




800 







5775 

5361 

4897 

4524 











3.53 

3.26 

3.09 

2.98 




900 








5730 

STSfl 













3.19 

3.03 





1000 








6048 














3.12 











Approximate Mean and Variability 


Subsolar Point 330 km ± 60 km 

Solar Zenith Angle 90 deg 850 km ± 300 km 

(AM Terminator) 

Solar Zenith Angle 120 deg 1200 km ± 550 km 

At the ionopause the electron density drops abruptly, from the ionospheric value 
4 3 

near 10 to less than 100 per cm . The electron temperature, however, remains 
unchanged across the ionopause (Fig. 5-23). 


5. 8 MAGNETOSPHERE 

Venus lacks a planetary dipole field, and so lacks a magnetosphere of the sort 
possessed by the Earth. However, the interaction between the solar wind and the 
Venus ionosphere produces magnetic effects. Figure 5-24 shows the average magnetic 
field from this cause as a function of radial distance from the planet center [5-50]. 

The maximum field is only a few tens of gammas, but the resulting magnetohydrody- 
namic plasma effects are very complex [5-46 through 5-49]. 



Figure 5-24. The average magnetic field strength for a field line 
intersecting the terminator at a radial distance r. 


Figure 5-25 shows schematically the form of the interaction of Venus and the 
solar wind. The "mantle" is the transition region between the ionosheath and the 
ionosphere. 


5.9 PROSPECTS FOR THE FUTURE 

Venus presents a rich world for investigation; Pioneer Venus results have just 
begun to open the planet for exploration. Futxire exploration by orbiters and soft- 
landers will provide much interesting and valuable information about our strange 
sister planet and, by analogy and contrast, will enhance our understanding of the 
Earth. 
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DISTANCE ALONG VENUS SUN lINE 

Figure 5-25. Schematic diagram of the interaction between 
the Venus ionosphere and the solar wind [5-51], 

Previous probes to the plant Venus which provided the majority of the data in 
this section are listed in Table 5-16. 
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TABLE 5-16. SPACE PROBES TO VENUS [5-9] 


Name of 
Mission 

Type of 
Mission 

Launch 

Arrival 

What was Measured 

Mariner- 2 

fly-by^ 
(35.000 km) 

26. 8. 62 

14.12.62 

Infrared and radio brightness tem- 
perature, magnetic field. 

Venera- 4 

descent 

probe 

+ fly-by 

12.6. 67 

18.10.67 

T , P , p (on altitudes 55- 25 km) , 
abundance of COg, Ng, Og, HgO , 
wind. 

La and 01 1300 %. (resonance 
emissions) . 

Mariner- 5 

fly-by 
(4000 km) 

14.6.67 

19.10.67 

Radio occultation. La. 

Venera- 5 

descent 

probe 

5.1.69 

16.5.69 

T, P (on altitudes 55-20 km), abun- 
dances of CO„, N„, 0„, H„0, wind. 

Venera- 6 

descent 

probe 

10.1. 69 

17.5.69 

The same. 

Venera- 7 

lander 

17.8.70 

15.12.70 

T (between 55 km and surface) . 

Venera- 8 

lander 

26. 1. 72 

22. 7.72 

T , P, abundance NHg, flux of solar 

scattered radiation (between 55 km 
and surface), wind, K, U, Th abun- 
dance in the ground. 

Mariner- 10 

fly-by 
(5800 km) 

8.11. 73 

5.2. 74 

UV-imaging of the clouds, IR Tg, UV- 
spectrum, radio occultation. 

Venera- 9 

lander 

+ orbiter 
(1545 km) 

8.6. 75 

22.10.75 

In the atmosphere T, P, chemical com- 
position (mass-spectrometer and 
narrow -band photometer) , solar flux, 
clouds (nephelometer and photometers) , 
wind. On the ground; imaging of the 
surface, density of ground,- K, U, Th, 
Imaging of the clouds. IR Tg narrow- 

band photometry between 0.4 and 2.5 
ym, polarimetry, spectrum of night- 
flow, La, radio occultation, solar wind, 
interaction, magnetic field. 

Venera- 10 

lander 
+ orbiter 
(1665 km) 

14.6.75 

25.10.75 

The same as Venera- 9. 

Pioneer- 

Venus 

Orbiter 
(150 km) 

20.5.78 

4.12.78 

Imaging of clouds, IR thermal sounding, 
UV-spectrum, polarimetry, radio 
occultation, upper atmosphere's in situ 
measurements drag, neutral and ionic 
composition , Tg , n^) . 

Pioneer - 
Venus 

Multiprobe 
(large probe 
+ 3 small 
probes + bus) 

8.8.78 

9.12.78 

T , P , p , solar and planetary fluxes , 
chemical composition (mass-spectrome- 
ter and gas chromatograph) , clouds 
particle size spectrometer and 
nephelometer . 

Venera- 11 

lander 

+ fly-by 
(40,000 km) 

9. 9. 78 

25.12.78 

T, P, p, chemical composition (mass- 
spectrometer) spectrum of scattered 
solar radiation, thunderstorms, clouds 
(nephelometer) . 

UV-spectrum, 

Venera- 12 

lander 

+ fly-by 
(40,000 km) 

12.9. 78 

21.12.78 

T, P, p chemical composition (mass- 
spectrometer and gas -chromatograph) , 
spectrum of scattering solar radiation, 
clouds (composition) . 

UV- Spectrum. 


a. For fly-by probes and orbiters the minimal distance to the surface is given. 
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SECTION 6. MARS 


6.1 INTRODUCTION 

Mars, the fourth planet from the Sun, was once considered to be a likely can- 
didate for extra-terrestrial life, but its atmosphere turns out to be thin and oxygen 
deficient; its surface is cold, dry, and apparently lifeless. Until the Mariner 4 flyby 
mission in 1965, ground-based observations provided the only information about the 
planet. NASA Special Publication 179, The Book of Mars [6-1], is a comprehensive 
guide to the early observations. 

Mariner 4, followed by Mariners 6 and 7 flybys in 1969, yielded only brief looks 
at the planetary surface. Then, taking advantage of the favorable 1971 launch window 
Mariner 9 was successfully put into an orbit around the planet. At the time the 
spacecraft arrived at the planet, its starface was totally obscured by a dust storm. 

The storm subsided and by the end of the mission, over 7000 pictures had been 
returned electronically to Earth. The Geology of Mars [6-2] reviews the results of 
these missions and early Soviet missions. 

Two Viking spacecraft, each consisting of an orbiter and a lander, reached Mars 
in 1976. A detailed account and analysis of the Viking missions is provided in Refer- 
ences 6-3 and 6-4. 

Knowledge of Mars has Increased enormously with planetary landers. Much more 
could be learned from samples returned from the planet; iindoubtedly , manned expe- 
ditions will eventually be sent to explore the planet and carry out observations on 
the surface and within the atmosphere. If the time ever comes when humans colonize 
the planets , Mars remains the most likely candidate. 


6.2 DYNAMIC PROPERTIES 

The dynamical properties ^of Mars are given in Table 6-1. 


6.2.1 Comments 

These data were adapted from The American Ephemeris and Nahtical Almanac 
[6-6]. The orbit eccentricity is comparatively large, which means that there is a 
significant difference between perihelion and aphelion distances. The Earth-to-Mars 
configuration varies considerably from one opposition to another, affecting oppor- 
tunities for spacecraft exploration. 

The obliquity of the axis of rotation and the eccentricity of the Mars orbit are 
believed to have varied widely over geologic time. The obliquity has ranged from 

5 

10.8 to 38 deg, with oscillatory periods as rapid as 10 years. The orbital eccentricity 
has varied between 0.004 (circular) and 0.141 [6-7]. 
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TABLE 6-1. DYNAMICAL PROPERTIES 


Parameter 

Value 

Comments 

Semi -major axis 

22. 794 X lo"^ km 

1.523691 AU 

Perihelion distance 

20. 665 X lo”^ km 

1.381398 AU 

Aphelion distance 

24. 918 X 10^ km 

1.666984 AU 

Eccentricity 

0.093387 

Earth: 0.017 

Orbit inclination 

1°50'59.28" 


Mean Orbital Velocity 

" 24.129 km sec-1 


Sidereal year 

686. 97964 days 


Sidereal day 

24^37"“22.663±0.002® 

16-5] 

Solar day (1 "sol") 

2 4^3 9"“ 35. 238® 

noon to noon 

Synodic period 

779. 94657 days 

time between Earth- 
Mars oppositions 

Obliquity of rotation 
axis 

25 deg 

Earth: 23. 5 deg 

Longitude of perihelion 
from vernal equinox 

253 deg 

Epoch 1980 

Earth-Mars opposition 
distance : max 

min 

10. 1 X 10^ km 
5.6 X 10 km g 



6.3 PHYSICAL DATA 

6.3.1 Mass, Figure, and Other Physical Properties 

Table 6-2 summarizes the physical data of Mars. 

6. 3. 2 Earth- Mars Configurations 

Table 6-3 lists the Mars opposition distances through 1990 while Figure 6-1 
shows some of the orbit transfer velocity data. 


6. 3, 3 Mapping 

The prime meridian for Mars passes through the small crater "Airy-0" within 
the large crater Airy at 5.19 deg S [6-10]. A very precise determination was made 
of the angle between Airy-0 and the vernal eqvhnox, allowing the longitudes of the 
Viking Landers to be determined accurately. 
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TABLE 6-2. PHYSICAL DATA 


Parameter 

Value 

Comments 

Radius 



Equatorial 

3397.2 km 


Polar 

3375. 5 km 


Mass 

6.418 X 10 kg 

0.1074 Earth 

Mass Ratio: Sun/ Mars 

3098700. 


Mean density 

3.933 g cm"^ 

Earth: 5.41 

Gravity at surface 

372. 52 cm sec ^ 

Earth: 980. 

Escape velocity 

5. 024 km sec ^ 

Earth: 11.2 

Flattening^ 



Dynamic 

0.00522 


Optical 

0. 00612 

Earth: 0.00335 

GM 

42828.32 km^ sec ^ 


'^2 

0.00196 

Earth: 0.001083 

Moment of inertia 
parameter (I /MR 2) 

0.365 

[6- 8] 

Magnetic dipole moment 

22 3 

less than 10 G cm 

see text 

Direction of rotation 

RA 317. 340±0.003 deg 

[6-5] 

axis (1976 epoch) 

Dec 57-710±0.002 deg 


Normal albedo 

0.1 to 0.4 


Surface temperatxare 
extremes 

130 to 300®K 


Surface atmospheric 
pressure 

5. 9 to 15. 0 mb 

varies seasonally 
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TABLE 6-3. MARS OPPOSITION DISTANCES THROUGH 1990 


Date 

Distance 
(10® km) 

Mars Season 

Aerocentric 
Solar Longitude 
North 

January 22, 1978 

98 


30 

Februeiry 25, 1980 

101 


80 

April 2, 1982 

94 

Summer 

130 

May 11, 1984 

80 

Summer 

175 

July 10, 1986 

65 

Fall 

225 

September 30, 1988 

59 

Fall 

270 

November 28, 1990 

82 

Winter 

320 


(o)EARTH launch 


V 

km/sec 



(b)MARS ARRIVAL 
V 

km/sec 


fc)MARS 0ER4RTURE 


(d) EARTH ARRIVAL 



V 

km/Sec 



V 

kmAec 


EARTH LAUNCH YEAR 



Figiire 6-1. Earth-to -Mars transfer orbit velocity increments [6-9], 


Topographic charts of Mars are available from the U.S. Geological Survey, as 
are detailed charts of Viking lander areas [6-11]. 

Figures 6- 2a and 6- 2b are global maps of Mars showing the major surface 
features and landforms [6-12]. 


6. 3. 4 Gravity Field 


Spherical harmonic coefficients for the martian gravity field have been derived 
to degree and order six and twelve [6-13,6-14]. Table 6-4 gives the coefficients for 
normalized spherical harmonics. (To convert to unnormalized format review paragraph 
2.7.1.) 

Except for very specialized purposes, only and would be necessary to 
use (paragraphs 2.3.2 and 3.3.2). 
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Figure 6- 2a. Major landforms on Mars [6-12]. 
Modern nomenclature as established by the lAU. 
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Legend : pi (permanent ice) 

ep (etched plains) 
pc (cratered plains) 
hk (hummocky-knobby) 
p (plains, undivided) 
cu (cratered terrain) 


Id (layered deposits) 

V (volcanic deposits) 
he (hummocky terrain, chaotic) 
c (channel deposits) 
g (grooved terrain) 
m (mountainous terrain) 


Figure 6-2b. Physiographic provinces of Mars [6-12]. 
(Based on Mariner 9 data.) 



TABLE 6-4. NORMALI2ED SPHERICAL HARMONIC 
COEFFICIENTS 


n 

m 

Cnn, X 10’ 

Snm X 10’ 

2 

0 

-8746 ±2* 


2 

2 

-850 ± 3 

497 ±2 

3 

0 

-126 ±9 


3 

1 

34 ±3 

243 ±6 

3 

2 

-157 ±3 

80 ±3 

3 

3 

352 ±4 

250 ±4 

4 

0 

61 ±3 


4 

1 

44 ±3 

- 33 ±4 

4 

2 

-9 ±2 

-88 ± 3 

4 

3 

68 ±2 

-1 ±3 

4 

4 

-2 ±4 

-121 ±5 

5 

0 

-23 ±4 


5 

1 

-2±4 

25 ±3 

5 

2 

-44 ±3 

-7±4 

5 

3 

30 ±3 

0± 3 

5 

. 4 

-48 ±6 

-36 ±4 

5 

5 

-50 ±6 

31 ± 5 

6 

0 

20 ±6 


6 

t 

24 ±4 

-1 ±4 

6 

2 

14 ±5 

20 ±3 

6 

3 

9 ± 4 

-7±4 

6 

4 

23 ±4 

31 ± 5 

6 

5 

21 ± 5 

2±6 

6 

6 

28 ±5 

2±5 


*Cjo = (-1.9557 ± 0.0004) X 10 ’ unnormalized. 


6,3.5 Magnetic Field 

The magnetic field of Mars has not been adequately investigated. Interaction of 
the planet with the solar wind may imply the existence of an intrinsic dipole of the 
23 3 -4 

order 10 G cm , (or 10 the terrestrial value) but this is at the present limit of 
detectability [6-15]. Such a small value is surprising, if Mars is assumed to have a 
core about the size of Mercury’s (as some believe) since its rotation rate is close to 
Earth’s. If there were an interior dynamo, it shoiild produce a planetary magnetic 
moment about 1000 times greater than the limit placed by the experimental data [6-16]. 
It may be that Mars has some residual surface magnetism and lacks an operating 
dynamo, like the Earth’s Moon. 

Questions about Mars' magnetic field will only be resolved by placing magnetome- 
ters on the planet or in low orbit, or by return of rock specimens. 


6.4 INTERIOR 
6. 4. 1 Core 

- 3 

The average density of Mars (3, 9 g cm ) is considerably less than Earth’s, 

suggesting less pronounced core differentiation. Past volcanism requires a molten 

interior, at least at one time. Estimates of core size and composition limits are given 

in Reference 6-17: pure Fe: radius 0. 4 R , pure FeS : radius 0. 6 R . The center 

m m 

of figure for Mars is displaced from the center of mass by 2.5 km toward 62 deg S, 
272 deg W. The geometric flattening (0.00612) is too large and the dynamic flattening 
(0.00522) is too small for Mars to be homogeneous and in hydrostatic equilibrium. 
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6. 4. 2 Crust 


Assuming zero crustal thickness at Hellas , minimum crustal thickness is 
estimated to be 28 ± 4 km [6-18], 


6.5 SURFACE 

The global topographic contour and albedo (Fig. 6-3) map of Mars shows several 
large scale features: the Tharsis Bulge centered at 101 deg W, 14 deg S; the large 

southern basins of Hellas and Argyre at 290 deg W, 50 deg S, and 40 deg W; and the 
huge equatorial canyons comprising Valles Marineris just east of the Tharsis summit. 



Figure 6-3. Topographic profile along martian equator. Zero circle 
represents topographic datum, and other ra(Uans represent 
elevation [6-19]* (Radial scale is greatly exaggerated.) 

The distribution of craters is asymmetrical over the planet. The northern 
third is relatively lightly cratered; the southern two-thirds is heavily cratered. 
Light cratering of the north implies it is a younger region with a different crustal 
type or thickness [6-20]. 

The mean altitude above datum for the two hemispheres also is very different 
as clearly shown in Figure 6-4. Note the topographic lows surrounding the polar 
regions (Fig. 6- 2b). 

4 

I 3 

s * 

S ' 

w 0 


-*o*-«o* - 30 " 0* air SO* 90 * 

kATITUOt 

Figure 6- 4. Average elevation of Mars latitude belts [6-2] . 



6-8 




6.5.1 Surface Morphology 


Mars exhibits a great variety of surface features: craters and cratered terrain, 

channels and fluvial features, tectonic features, volcanic features, aeolian features, 
and the polar regions. In this discussion, special attention is given to the Viking 
Lander sites, since detailed photographic information is available for them. 


6.5.2 Cratered Terrain and Craters 

Mars can be transected by a plane inclined approximately 40 deg to the equator 
such that over 90 percent of the densely cratered terrain lies on the southerly hemi- 
sphere so defined [6-21]. 

The sparsely cratered plains, common in northern latitude regions, appear to 
antedate the events which bombarded the ancient densely cratered terrain. They are, 
nonetheless, fairly old (3 to 4 billion years). On the sparsely cratered plains are a 
complex variety of features that have been attributed to the interplay of impact, 
volcanic, aeolian, and ice- related processes. 

The largest impact crater on Mars is Hellas Planitia, centered at 40 deg S and 
290 deg W. The basin is 1600 km in diameter and 6-km deep, more than 4 km below 
the datum at its lowest point. Roughly circular in outline, Hellas Planitia is undoub- 
tedly an original impact structure, but it has been highly degraded. Ar gyre Planitia 
is a similar, although smaller, impact crater [6-2]. 

The majority of moderately sized craters (< 100 km) have ejecta blankets dis- 
tinctively different from those on the Moon or Mercury. They show one or more 
layers of ejecta, each with a distinct outer ridge or rampart, and were given the 
name "rampart craters." Their characteristic ejecta configuration appears to be due 
to a radial outward flow of ground-hugging debris after crater formation, perhaps 
from melting ground ice or permafrost [6-21]. 

Pedestal craters , another featiire peculiar to Mars , are found primarily in the 
northern hemisphere between 30 and 70 deg N. Typically, the crater is surrounded 
by a low, nearly circular platform or pedestal. Pedestal craters are believed to be 
the result of repeated deposition and erosion of extensive debris layers by wind. 

Such events appear to have occurred widely in this region of the northern hemisphere. 
The ejecta blanket is believed to have protected the underlying layers from removal, 
resulting in the pedestals [6-21]. , 

The densely cratered hemisphere has been heavily modified into regions of 
chaotic, fretted, and knobby terrain. 

Chaotic terrain refers to areas that have apparently sunk below the surrounding 
terrain to form irregularly shaped, jostled blocks. Large channels, often character- 
istic of fluvial flow, originate at many of the chaos areas. According to current 
theory, chaotic terrain was created when the ground collapsed after melting of sub- 
surface ice. 

Fretted terrain is characteristic of the boundary between the densely cratered 
highlands and the sparsely cratered lowlands. Typically, this terrain is a complex 
interfingering of high standing remnants of the densely cratered terrain and inter- 
crater plains and, at a low elevation (usually 1 to 2 km), the more sparsely cratered 
plains. Throughout the fretted terrain, debris flow occurs at the boundary between 
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the cratered plateaus and the lower plains. The mechanisms of debris flow may be 
linked to the melting of ground ice. 

Knobby terrain represents the transition regions between the fretted terrain 
and the sparsely cratered plains to the north. Mechanisms responsible for the erosion 
and/or reformation of the northern plains seem to have reshaped the cratered plateaus 
of the fretted terrain into smaller, and form-rounded, equi dimensional hills Charac- 
teristic of the knobby terrain. 


6. 5.2.1 Crater Statistics 

The observed distribution of craters is bimodal (Fig. 6-5), with two main types:; 
larger craters (diameter > 5 km) which tend to be shallow with relatively low rims; 
and smaller craters that are generally bowl-shaped. Crater distribution varies with 
size. For very large craters (>64 km) the distribution is nearly saturated (like the 
lunar highlands). 



Figure 6-5. Crater distributions on Mars [6-2]. 


6. 5. 2.2 Block Distributions 

Block distribution data have bpen derived for the Viking 1 and 2 landing sites 
from the detailed pictures sent back. Figure 6-6 shows rocks in the neighborhood of 
Viking Lander 2. Figure 6-7 shows derived size distributions. The similarity of the 
two sites, which differ geologically, lends some measure of confidence to regional 
extrapolations [6-22]. 
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Figixre 6-6. Viking Lander 2 site 16-12] 
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6.5.3 Channels and Fluvial Features [6-21] 


The observed channels suggest formation by running water. (There is no 
connection between these features and the famous ’’canals," which proved to be an 
optical illusion.) Liquid water cannot exist on Mars at present since the atmospheric 
pressure is below the triple-point; therefore, either the channels were formed by other 
processes, such as faulting, lava flows, or wind action, or atmospheric conditions were 
very different sometime in the past. There are three main t37pes ; fretted channels, 
runoff channels, and outflow channels. 

Fretted channels are flat-floored sinuous valleys. They are mostly in fretted 
terrain which occurs in two latitudinal bands, 25 deg wide, centered on 40 deg N and 
45 deg S. 

Runoff channels typically start small and increase in size and depth downstream 
and have a characteristic V-shaped cross section. They almost exclusively occur in 
densely cratered terrain. The extensive, well-integrated tributary networks asso- 
ciated with these channels are strong evidence for water erosion in some past epoch. 

The outflow channels appear to be the result of catastrophic flooding and most 
originate in chaotic terrain. They are deepest close to their sources and extend for 
hundreds of kilometers with widths typically tens of kilometers. 


6. 5.3.1 Tectonic Features: Tharsis and Valles Marineris 

Compared to Earth, Mars, is seismically inactive at present, though there is 
ample eAudence for tectonic activity in the past. Most notable is the Tharsis Bulge, 
a roughly circular feature covering nearly a quarter of the planet’s surface, with a 
summit at an elevation over 11 km above the datum. Most of the large shield volcanoes 
are found atop Tharsis. Surrounding Tharsis is a radial pattern of faults. As only 
the oldest surfaces retain the fault pattern, Tharsis is believed to have been formed 
early in the planet’s' history. 

Valles Marineris is a collection of steep walled canyons extending east from the 
base of Tharsis for about 4500 km just south of the equator. Different canyon sys- 
tems merge and in places reach depths of over 7 km and widths of 600 km [6-21] . 

The general alignment of Valles Marineris, the presence of fault scarps, and the 
absence of erosion deposits ..strongly suggest the canyons were produced by faulting, 
probably related to the formation of the Tharsis Biilge [6-24]. 


6,5.4 Volcanic Features 

Mars has many volcanic features. Large shield volcanoes, roughly circular 
mountains with central caldera and gently sloping (<6 deg) flanks, are primarily 
clustered in regions of Tharsis and Elysium. The largest, Olympus Mons, is centered 
at 18 deg N, 133 deg W, with summit 22 km above the surrounding plains, 27 km 
above the datum. 

At the base of the central edifice is an outward facing escarpment, 550 km in 
diameter and in places up to 6-km high. In many places laya flows drape over the 
escarpment and extend far beyond it, so the mean diameter of the shield is approxi- 
mately 700 km. On the Tharsis Bulge are three other large shield volcanoes, Arisia 
Mons, Pavonis Mons, and Ascraeus Mons which resemble Olympus Mons, as well as 
terrestrial volcanic shields such as Mauna Loa. 
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On the Elysium Bulge are three volcanoes, Elysium Mons, Hecates Thalus, and 
Albor Thalus. The largest, Elysium Mons, is a shield volcano, 14-km high, 170-km 
wide, comparable to the Tharsis shields, though with significantly steeper flanks 
(10 to 12 deg). 

Numerous volcanic features, which range in size down to the limit of resolution 
of the available photography, have been identified outside the two main provinces of 
Tharsis and Elysium. 


6.5.5 Aeolian Features 


Mars wind action is well known from the vast annual dust storms. Yet their 
erosional effect appears to be small; they mostly recirculate already eroded debris 
[ 6 - 21 ]. 


Extensive dune fields exist on Mars. Around the north pole is a nearly con- 
tinuous belt of dunes which includes both transverse and barchanoid types [6-12]. 

The Same pattern persists over large areas, suggesting a constant direction and 
strength for the wind that formed them. In a few areas several wind systems are 
implied. Dime fields are also common around the south pole, but are mostly restricted 
to the interiors of large craters. The mechanisms for forming sand-size particles are 
still unknown. The windblown dust particles deposited at the Viking Lander sites 
would not form dunes. It is not known whether dune forming processes are still 
active. 

Dust storms primarily redistribxite fine debris. This redistribution seems to 
explain the seasonal color changes that were once thought to be evidence of vegeta- 
tion. (The bright regions are attributed to accumiilations of a thin layer of fine 
debris, while the dark regions are the underlying solid surface.) The wind-driven 
debris deposits are thin, only centimeters or even millimeters thick. 

A peculiar aeolian feature is the polar laminated terrain discussed in the follow- 
ing paragraph. 


6.5.6 Polar Regions 

The distinctive polar caps of Mars are easily visible from the Earth, and their 
seasonal advance and retreat have been observed for more than two centuries. The 
two caps are not symmetrical, they differ in maximum and minimum size, as well as in 
composition. The northern cap esrtends down to between 50 and 64 deg N at maxi- 
mum and retreats to 87 deg N at minimum. The. south polar cap extends down to 
40 deg S at maximum, retreats more rapidly in the spring, and nearly disappears in 
summer. Both seasonal caps are CO 2 ice; the northern residual cap is H 2 O ice, the 

residual south polar cap may be mostly COg [6-25]. Their differences have been 

attributed to a number of factors: the planet’s eccentric orbit and near coincidence 

of southern summer solstice and perihelion; high global dust-storm activity during 
southern summer, resiilting in cooler surface temperattmes and more isothermal 
atmospheric temperature profiles ; and the different elevation and surrounding 
topography of the two polar regions. 
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Sequences of layered deposits form the unique laminated terrain of the martian 
polar regions. They apparently are thick accumulations of volatiles (most likely water- 
ice [6-21,6-26] and airborne dust [6-27 through 6-29]. The layered structure is believed 
to be due to climatic changes that have resvilted either from orbit variations [6-7] or 
from variations in solar luminosity [6-29]. If the layered structure has, indeed, 
arisen from solar variations, its study cotild lead to better understanding of long-term 
climatic changes on the Earth. 

The polar laminated terrain has been sculptured into a series of stair-like 
incised valleys and gently sloping escarpments, 100- to 1000-m high, 3- to 10-km wide, 
aijid hundreds of kilometers long [6-21,6-30]. The origin of these valleys and escarp- 
ments remains controversial [6-27,6-31,6-32,6-33,6-34]. 


6.5.7 Permafrost 

There is evidence that a thick permafrost (mean annual temperature less than 
273°K) covers the entire planet and extends to a depth of about I km at the equator 
and several kilometers at the pole. Within 40 deg of the equator, however, the 
ground is permanently dehydrated to depths of at least 1 m. Existence of water-ice 
below these depths depends on the permeability of the soil [6-12,6-35]. 


6.5.8 The Viking Lander Sites 

Both Viking Lander sites are in the sparsely cratered northern lowlands. 

Viking 1 landed on Chryse Planitia, 47. 97 deg W, 22.27 deg N, at about 16:00 local 
time. Viking 2 landed on Utopia Planitia, 225.75 deg W, 47.67 deg N, at about 
10:00 local time. Both landed during local summer. The Viking 1 and Viking 2 sites 
are 1.5 km and 2.5 km below the mean pressure level of the planet, respectively 
[6-36]. 

The choice of lander sites was based upon a requirement for smooth terrain. 

The two sites are topographically similar: low rolling hiUs, mostly covered by fine- 
grain debris, and only 8 percent of each area contained rocks 3.5 cm or larger 
(Fig. 6-7). 

The surface samples analyzed from both sites are similar and probably repre- 
sentative of fine debris everywhere on the planet. The atmospheric and meteorological 
observations from the two sites are also likely to be indicative of the planet as a 
whole. 


6. 5. 9 Composition of Stirface Material 

The X-Ray Fluorescence Experiment analyzed the elemental composition of the 
loose debris at the Viking landing site [6-37]. The experiment could not detect eie- 
ipents with atomic numbers less than 12 (Mg). The determination for a typical Mars 
sample appears below , and the assumed mineralogical composition is given in Tables 
6-5 and 6-6 [6-38 through 6-40]. (Phosphorus is not listed becaiase its signal is hidden 
by S and Si, although minor amounts may be present.) 

The fine-grained debris probably consists mostly of Fe-rich clays with minor 
amounts of Meserite (MgSO^. H 2 O), calcite (CaCO^), and rutile (Ti0 2 ). The surface 
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TABLE 6-5. ELEMENTAL COMPOSITION OF THE' VIKING 1 

LANDER SITE [6-37] 


Element 

Percent by Mass 

Mg 

5. 0 ± 2. 5 

A1 

3.0 + 0. 9 

Si 

20. 9 ±2.5 

s 

3.1 ± 0.5 

Cl 

0.7 ± 0.3 

K 

<0. 25 

Ca 

4.0 ± 0. 8 

Ti 

0. 51 ± 0.2 

Fe 

12. 7 ± 2. 0 

L* 

50.1 ± 4. 3 ■ 

y** 

8 . 4 ± 7.8 

Rb 

<30 ppm 

Sr 

60 ± 30 ppm 

Y 

70 ± 30 ppm 

Zr 

<30 ppm - 


*L is the sum of all elements not directly determined. 

**If the detected elements are all present as their common oxides 
(Cl excepted) then X is the sum of components not directly 
detected, including H 2 O , NaO, CO 25 and NO^. 

TABLE 6 - 6 . CHEMISTRY OF THE VIKING 1 LANDER SITE [6-39,6-40] 


Compound 

Percent Mags 

Si 0 2 

44. 7 

AI 2 O 3 

5.7 

^®2^3 

18.2 : 

MgO 

8 . 3 

CaO 

5.6 

K2O 

<0.3 

Ti 02 

0 . 9 

SO 3 

7.7- 

Cl 

0.7 

SUM 

91.8 
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material contains significant quantities of adsorbed volatiles. When heated to 500°C, 
the sample material released about 1 percent water by weight and from 50 to 100 ppm 

-7 -1 

CO 2 [6-41]. Additionally, trace amounts of O 2 (less than 7 x 10 moles g ) and 

other gases were released upon humidification [6-42], There may be significant 
amounts of CO 2 and H 2 O adsorbed into the top regolith layers, perhaps as much as 

400 g cm~^ adsorbed CO„, which is equivalent to a surface pressure at equilibrium of 
over 100 mb [6-43]. 

The results of the Viking biological experiments support the suggestion that the 
martian top regolith contains highly oxidizing compounds [6-44], The hypothesized 
peroxides and superoxides could be produced by the ultraviolet radiation which 
reaches the martian surface Mnce there is no blocking ozone layer in the atmosphere. 


6.5.10 Optical, Thermal, and Mechanical Soil Properties 

Optical, thermal, and mechanical soil properties of the Mars top regolith are 
given in Table 6-7 while those for the Viking 1 landing site are compared in Table 6-8 
with the Moon. 


TABLE 6-7. PROPERTIES OF THE MARS TOP-REGOLITH* 


Parameter 

Nominal Value 

Thermal Conductivity 

2 to 20 X 10 ^ cal sec ^ cm ^ K ^ 

Specific Heat 

0.15 to 0.19 cal g"^ K“^ 

Emissivity 

0.90 to 0.98 

2 i -i 

Thermal Inertia 

100 to 600 cm s^K7cal: (k p c) " 

Albedo 

0.2 to 0. 4 

Bolo metric Albedo 

0.43 to 0.93* 

Dielectric Constant 

2.3 to 3. 5 

Bulk Density (porosity) 

- 3 

1 to 1 . 8 g cm 

Penetration Resistance 

0. 3 N cm ^ cm ^ 

Cohesion 

0.01 to 0 . 1 N cm'^ 

Adhesion 

10 “^ to 10'^ N cm'^ 

Coefficient of Sliding Friction 

0. 55 to 0. 65 

Density 

3. 933 g CP ^ 


^Equatorial Regions [6-45]. 
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TABLE 6-8. COMPARISON OF VIKING 1 LANDING SITE 
AND LUNAR SOIL PROPERTIES 



Mars (Viking 1) 

Lunar Depth 

Property * 

Sandy 

Flats 

Rocky 

Flats 

0 to 0. 01 m 0. 1 to 3 m 

Bulk density (porosity) 
(g/cm^) 

Particle si ze (siirface 
and near surface) 

1 to 1.6 

1.8 

1.0 to 1.3 1.5 to 2.1 

10 to 100 pm (;%)■ 

60 

30 

30 to 60 

100 to 2000 pm (%) 

10 

30 

30 to 35 

Angle of internal 
friction (deg) 

20 to 40 

40 to 45 

35 to 50 

Penetration resistance 
(dyne /cm 2 /cm) 

3 X 10^ 

6 X 10^ 

3 X 10^ 

2 

Cohesion (dyne /cm ) 

3 4 

lO*^ to 10 

10^ 

2 

Adhesion (dyne/cm ) 

10^ to 10^ 

10^ to 10^ 

Coefficient of sliding 
friction 

0. 55 to 0. 65 

0. 5 to 1 


*Soil properties deduced from Viking 1 data. Lunar soil is included for 
comparison and because the nominal engineering design soil model was based 
on the lunar soil [6-46]. 


The Viking Infra-Red Thermal Mapper experiment aboard the orbiters indicates 
a wide range of surface thermal properties and large atypical regions. Planetwide 
measurements of the infrared emission and reflection yielded information on the thermal 
inertia of the surface [6-45]. Temperatures observed during the Viking primary 
mission ranged from 130 to 290°K. 


6. 6 ATMOSPHERE 

The general temperature structure of the atmosphere is shown in Figure 6-8. 


6. 6. 1 The Lower Atmosphere, Temperature and Pressure 

The mean temperature profile for the martian lower atmosphere at summer tem- 
perate latitudes, derived from data gathered during the Viking entries, is shown in 
Figure 6-9 [6-36]. 

Model atmospheres for "warm” and "cool” conditions, ±10°K from nominal, 
coupled with high and low surface pressures, are shown in Table 6-9 [6-36J. 
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Figure 6-8. General temperature structure of the atmosphere of Mars. 



Figure 6-9. Nominal mean temperature profile for summer temperate latitudes (6-36). 
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TABLE 6-9. SUMMER. MID-LATITUDE , DAILY-MEAN ATMOSPHERE 
OF MARS (COOL AND WARM MODELS) (6-36) 



Cool, Low Pressure Model 

Warm, High Pressure Model 

z, km 

T. “k 

p, mb 

P. kg/m^ 

mm 

T. *K 

p, mb 

P, kg/m® 

P/Po 

0 

204 

5.9 

1.51 X 10*^ 

l.OOO 

224 

7,8 

1.82x 10*^ 

1.000 

4 

204 

4.03 

1.03 

0.683 

224 

5.51 

1.29 

0. 706 

8 

199 

2. 74 

7. 20 X 10*® 

0.464 

219 

4.09 

9.77 X 10*® 

0.524 

12 

191 

1.84 

5.04 

0.312 

211 

2.85 

7.07 

0. 366 

16 

185 

1.22 

3.45 

0. 207 

205 

1.97 

5.01 

0. 252 

20 

178 

7.96 X 10*‘ 

2.34 

0.135 

198 

1.34 

3.54 

0. 172 

24 

173 

5. 13 

1.55 

8.70 X 10*^ 

193 

9.03 X 10** 

2.45 

0. 116 

28 

168 

3.27 

1.02 

5.54 

188 

6.03 

1.68 

7.73 X 10*^ 

32 

163 

2.06 

6.60 X 10*'* 

3.49 

183 

3,99 

1. 14 

5. 11 

36 

158 

1.28 

4.23 

2.17 

178 

2.61 

7.67 X 10** 

3.35 

40 

152 

7.81 X 10*^ 

2.69 

1.32 

172 

1.69 

5. 13 

2. 16 

44 

148 

4. 70 

1.66 

7.96x 10*® 

168 

1.08 

3.36 

1.38 

48 

144 

2. 79 

1.01 

4.73 

164 

6.83 X 10*^ 

2. 18 

8.75 X 10*® 

52 

140 

1.64 

6.12X 10*® 

2. 78 

160 

4.28 

1.40 

5.48 

56 

137 

9.49 X 10’^ 

3.62 

1.61 

. 157 

2.65 

8.84X 10*® 

3.40 

60 

134 

5.44 

2.12 

9. 22 X 10** 

154 

1.63 

5.55 

2.09 

64 

132 

3.09 

1.22 

5.23 

152 

9.99 X 10‘® 

3.44 

1.28 

68 

130 

1.74 

7.01 X 10*® 

2.95 

ISO 

6.08 

2. 12 

7.79X 10** 

72 

129 

9. 76 X 10"^ 

3.96 

1.65 

149 

3.68 

1.29 

4.72 

76 

129 

5.47 

2. 22 

9.27 X 10*® 

149 

2. 23 

7.83 X 10*® 

2.86 

80 

129 

3.07 

1.24 

5.20 

149 

1.35 

4.75 

1.73 

84 

129 

1.72 

6.99 X 10*^ 

2.92 

149 

8.21 X 10** 

2.88 

1.05 

88 

129 

9.70 X 10'® 

3.93 

1.64 

149 

4. 99 

1.75 

6.39 X 10*® 

92 

129 

5.46 

2. 22 

9. 26 X 10*® 

149 

3.03 

1.07 

3. 89 

96 

129 

3.08 

1. 25 

5.22 

149 

1.85 

6.49 X 10*’ 

2. 37 

100 

129 

1.74 

7.06 X 10’® 

2.95 

149 

1. 13 

3.96 

1.44 


6-19 

















To represent variations in the summertime atmosphere, temperature profiles were 
displaced by ±10°K from the nominal. This variation everywhere exceeds the maxi- 
mum deviation of the Viking 1 and ViMng 2 mean curves from the nominal by at least 
a factor of 2. The range ±10 deg about the nominal is chosen as a reasonable bound 
for warm and cool models and appears suitable down to the near- surface level and at 
the surface to near the end of the northern summer season. 

The cool model has pressure falling most rapidly with altitude and, hence, is 
combined with a selected minimum value of surface pressure, 5.9 mb, to define a cool, 
low-pressure model. This model would tend to occur near 160- deg solar longitude at 
northern latitudes about 45 deg. By similar reasoning, the warm model is associated 
with a maximum level of surface pressure, 7.8 mb, to define a warm high-pressure 
model. The 7. 8-mb maximum occurs in southern early summer around 280-deg solar 
longitude and might tend to occur around 25-deg south latitude where the solar 
incidence at noon is normal. 

The selected combinations of surface pressure and temperature profile do inot 
yield the greatest range in surface density which, in the models, ranges from 0.0151 

3 3 

to 0.0182 kg7m or ±9.3 percent about the mean of 0.0166 kg/m . The maximum range 

of surface density is obtained by combining the cool temperature limit with the maxi- 
mum sijrface pressure and vice versa. This yields surface (datum level) densities from 

3 

0.0137 to 0.0200 kg/m . (The lower limit might occur near solar longitude 160 deg at 
subsolar latitudes; the upper limit, at Ls = 280 deg, 6 = 45 deg S.) However, the 
selected models do yield the maximum density range at high altitudes where entry 
vehicles must decelerate, showing a range of 5.6/1 at 100 km [6-36]. 


6. 6. 2 Lower Atmosphere Composition 

The composition of the lower atmosphere as determined by mass spectrometers 
on the Viking landers is given in Table 6-10 while Table 6-11 gives isotopic ratios. 
The major constituent is CO 2 * Its abundance varies during the martian year by 

26 percent, sUghtly changing the mixing ratios of the other gases [6-47]. 

TABLE 6-10. COMPOSITION OF THE MARTIAN LOWER ATMOSPHERE [6-48] 


Gas 

Percent Volume 

Carbon Dioxide (CO 2 ) 

95.32 percent 

Nitrogen (N 2 ) 

2. 7 percent 

Argon (Ar) 

1. 6 percent 

Oxygen (O 2 ) 

0. 13 percent 

Carbon Monoxide (CO) 

0.07 percent 

water Vapor (H 2 O) 

0.03 percent variable 

Neon (Ne) 

2. 5 ppm 

Krypton (Kr) 

0. 3 ppm 

Xenon (Xe) 

0. 08 ppm 
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TABLE 6-11. ISOTOPE RATIOS IN ATMOSPHERIC GASES (6-36) 


Ratio 

Mars 

Earth 

12c/13c 

90 

89 

IGq/ISo 

500 

499 

l^N/l^N 

165 . 

277 

^^Ar./^^Ar 

3000 

292 


2.5 

0. 97 


The traces of ozone detected by Mariner 9 are not enough to block solar \iltra- 
■violet radiation, though ozone may influence the atmospheric radiation balance [6-49]. 

A maximum of 57 ym STP over the polar regions during winter was reported [6-50]. 
Typical values range from 3 ym in summer to 16 ym in winter over the polar regions. 

The average amount of atmospheric water vapor is 0.03 percent by volume, or 
approximately 10 precipitable microns. The water vapor concentration is closely linked 
to the temperature distribution: during most of the year and at most latitudes the 

atmosphere holds all the water possible (i.e., the relative humidity is 100 percent at 
night) [6-51]. 

The isotopic ratios in atmospheric gases, compared to terrestrial values, appear 
in Table 6-11. The anomalous nitrogen and argon ratios can be used to put con- 
straints on the history of volatiles in the atmosphere of the planet [6-38,6-52,6-53, 
6-54,6-55]. 


6. 6. 3 Meteorology 

6. 6. 3. 1 Meteorological Parameters and Length of Seasons 

Tables 6-12 and 6-13 list global meteorological properties and compare the 
seasons in the northern and southern hemispheres. The most striking feature is the 
seasonal variation in atmospheric pressure due to condensation of CO 2* 


6. 6. 3. 2 Pressure Variations 

Figure 6-10 gives the pressxire at both lander sites for an entire martian year 
and clearly shows the asymmetry between the poles [6-56]. The north cap does not 
exchange as much mass as the south polar cap, primarily because of the eccentricity 
of the orbit (0.093). The near coincidence of the perihelion (aerocentric longitude, 

LS = 253 deg) and northern winter solstice (Ls = 270 deg) results in more extreme 
seasons in the southern hemisphere. Note that the Viking 2 site is nearly 1 km below 
the Viking 1 site, referred to the mean datum. 
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TABLE 6-12. METEOROLOGICAL PROPERTIES 


Parameter 

Nominal Value 

Note or Range 

Pressure 

f- 8 mb- 

7 to 9 mb 

Temperature 

215°K 

130 to 290°K 

Temperature Lapse Rate 

2°K km"^ 


Adiabatic Lapse Rate 

4. 5®K km"^ 


Scale Height 

11 km 


Tropopause Height 

40 km 


Turbopause Height 
Wind Speed 

120 km 

0 to 50 m sec ^ 

above boundary layer 

Solar Radiation 
mean 

perihelion 
ape he lion 

590.0 W m‘^ 

718.0 

493.0 

based on the solar 
constant at 1 AU of 
1371 ± 5 W m"2 

Albedo 
average 
ice caps 

Solar UV Flux 

0.2 to 0.4 
0.7 

4 -2 - 1 

10 erg cm sec 

200-300 ran, top of 
atmosphere, siibsolar 
point, perihelion. 


TABLE 6-13. DURATIONS OF SEASONS ON MARS 


Season.) 

Northern 
Hemisphere .r 

: f ‘ ■ 

Southern 

Hemisphere 

Length 

Days 

Sols 

spring 

autumn 

199 

194 

summer 

winter 

183 

178 

autumn 

spring 

, 147 

143 

winter 

summer 

158 

154 



687 

669 










Figure 6-10, Surface pressure variation over one Mars year [6-56], 
Center section: Daily mean pressures at the two landers for 700 

sols. Upper and lower sections: Standard deviations of 

pressvire within each sol. All pressures and standard 
deviations are in mb. The abscissa is time measured in 
sols. The scale labeled Ls is the aerocentric (refers to 
Mars orbit around Sun) longitude of the Sun. The dates 
of initiation of two global dust storms are marked. 

Gaps are due to lost data. 


6. 6. 3. 3 Diixrnal Surface Temperatures 

Large diurnal temperature variations were observed at both landing sites. 

Fig. 6-lla shows the air temperature for one sol (24.46 hr) at both Viking landing 
sites early in the mission dxiring martian summer, as measured 1.6vm above the sur- 
face [6-57]. The Viking orbiter also measured surface temperatures over the entire 
planet which ranged from 130 to 290°K [6-45], 

Minimum air temperature differs only slightly from minimum soil temperatures 
while maximum air temperatures are considerably lower than surface temper at vires. 

Fig, 6- lib shows a typical diurnal temperature cycle for equatorial eqviinox conditions 
[6-58]. 
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Figure 11a. Diurnal temperatxire variations at Viking sites [6-57]. 
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Figure lib. Di\irnal temperature variations for Mars air and 
surface; at the equator at equinox [6-58]. 


6. 6. 3. 4 Winds, Aerosols, and Dust Storms 

The winds at the Viking lander sites ranged from 2 to 7 m sec with strong 
diurnal and seasonal variation in both direction and magnitude. Winds above the 

boundary layer may be as high as 50 m sec [6-57]. (The daily winds result from 
meteorological disturbances and are not directly due to the pole-to-pole transfer of 
CO2.) 

Dust storms have long been observed. They tend to occur near perihelion (Ls - 
253 deg). Because of the eccentricity of the martian orbit, the insolation at perihelion 
is fully 45 percent stronger than at aphelion. Additionally, perihelion nearly coincides 
with southern summer solstice (Ls = 270 deg). The martian atmosphere tends to be in 
radiative equilibrium, and this intense period of summertime radiation coupled with 
topographic variations apparently triggers the global dust storms. Storms may last 
for several months. The opacity of the dust cloud in visible wavelengths may be 
quite high, reaching 5 optical depths in the visible [6-59]. 
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The particles in the air follow the size distribution: 


n(r) = Cr exp [-4r/r^] 


where n(r) is the number density of particles of radius r, and r , the mode radius, 

rn 

equals 0.4 microns [6-60]. The total column mass loading during a global storm is 
- 3 - 2 

typically 10 g cm . The single scattering properties of the particles are listed 
below, averaged over the solar spectrum [6-60]. 

Average Scattering Extinction 

Single Scattering Asymmetry Parameter Efficiency 

Albedo (jj cos 0 Q + 

0Xl 

0.86 0.79, 2.74 

The slightly red color of the martian sky observed by the Viking cameras is 
due to this dust in the atmosphere; even on "clear'* days the opacity exceeds 0.18 
[6-61]. 

It is interesting to note that the two global dust storms observed , identified by 
orbital imaging, resulted in only minor redistribution and accumulation of debris at the 
lander sites [6-62]. Martian dust storms are not large walls of fast moving debris, as 
often imagined. While martian wind speeds may be typically ten times those of terres- 
trial winds, the density of the martian air at the surface is 120 times less. The drag 
of the martian winds is, in fact, comparable to that of typical terrestrial winds because 
wind drag is roughly proportional to air density times wind velocity squared. However, 
the dust storms and transport of dust into the polar regions are a major factor in 
martian climate. 


6. 6. 3. 5 Clouds 

i 

Clouds, frequently seen in orbiter photographs, can generally be divided into 
two types, dust and condensate ices (water and carbon dioxide) . 

Yellow clouds are generally associated with dust storms. The storms typically 
begin near perihelion and the yellow cloud may then qmckly obscure the entire surface. 

A wide variety of white condensate cloud types occur. They are composed pri- 
marily of water, but a few have been observed at altitudes high enough (i.e., 80 km) 
to imply CO 2 composition. 

Large , well-developed spiral cloud formations seen by the Viking orbiters 
appeared to be at 65 to 81 deg N and at about 6 to 7 km ^titude [6-63]. 

Topographic clouds are produced by air flow over craters, ridges, and moun- 
tains [6-64]. The clouds associated with the Tharsis summits are particularly striking, 
covering large areas and obscuring all but the peaks. 

Frost and early morning fog layers are common [6-61]. They are often seen in 
orbiter pictures in the bottom of canyons and were detected at both lander sites 
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during the summer. Hoods have been observed over both poles in late winter with 
the north polar hood being much larger in temporal and spatial extent. 


6. 6. 3. 6 The "Blue Clearing" 

Astronomers have noted that the contrast of the martian surface decreases 
toward the blue end of the spectrum, and that occasionally a "clearing" occurs. 
Though this phenomenon is well documented, nothing that can be related to it was 
observed during spacecraft encounters. A discussion of recent ground-based obser- 
vations appear in Reference 6-65. 


6.7 UPPER ATMOSPHERE AND EXOSPHERE 

The atmosphere is well mixed below 120 km. Above this altitude (the turbo- 
pause) composition varies according to the molecular weights of individual constituents. 
Carbon dioxide remains dominant up to at least 180 km. 

The eddy diffusion coefficient at 120 km, as determined by height profiles of 
7-1 

Ar and Ng, is 5 x 10 cm sec . 

The thermal structure of the upper atmosphere is complex and variable. The 
Viking 1 and 2 entries observed average temperatures near 200®K, which is consider- 
ably colder than the 300 to 350®K temperatures observed by Mariners 6, 7, and 9 
[6-50,6-66]. 

The exosphere and the escape of gases are important factors in determining 

8 "2 “1 

atmospheric evolution. The escape rate of H 2 is of the order of 10 cm sec , 

twice the escape rate of 0„. The estimated average of 10 precipitable microns of 

^ 4, 

water vapor in the atmosphere could support this escape rate for only 10 years. 

Isotopic separation of oxygen and nitrogen may occur in the martian exosphere, which 

could account for the observed nitrogen isotopic ratios. The oxygen ratios reflect a 

large reservoir of exchangeable oi^gen (6-52) . 


6.8 IONOSPHERE 

The martian ionosphere extends upwards from 120 km; CO 2 is a dominant con- 
stituent, and free electrons are mainly produced by its photoionization. Mars does 
not have an E or F 2 ionosphere layer but does have an peak formed near unit 

optical depth for solar ultraviolet radiation in the wavelength range from 200 to 900 
Just below 130- km altitude the F^ ion concentration reaches a peak of approxi- 

mately 10 cm of which 90 percent is O 2 and 10 percent is COg • At higher 

altitudes, ions were detected with a peak concentration near 225 km of less than 

10 cm . Viking 1 measured ion temperatures of approximately 150°K near the F^ 

peak increasing to an apparent exospheric temperature of 210°K near 175 km. Above 
this altitude departures from thermal equilibrium with the neutral gas occur, and the 
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ion temperatxare increases rapidly to over 1000°K at 250 km [6-67], Figure 6-12 shows 
profiles of ion concentrations compared with theoretical chemical model results. 



Figure 6-12. A comparison of the measured ion concentration from 
Viking 1 profiles and a theoretical model ionosphere in 
chemical equilibrium [6-67]. 


6. 9 MAGNETOSPHERE 

The interaction of Mars with the solar wind appears to require a magneto spheric 
component in addition to an atmospheric component. Intriligator and Smith [6-15] 
interpret spacecraft results as indicating a situation between "Venus-hke" (effectively 
no intrinsic magnetic field, interaction directly with ionosphere) and "Earth-like" 
(interaction dominated by the intrinsic magnetic field). 


6. 10 SATELLITES 

Asaph HaU discovered and named Mars' two natural satellites Phobos (fear) and 
Deimos (terror) in 1877. Their orbital elements and physical parameters are listed in 
Table 6-14 [6-68]. 

The orbit of Phobos is just outside the Roche limit of Mars and inside a sta- 
tionary orbit (analogous to geostationary orbit) . An observer on the surface of Mars 
would see Phobos rise in the west and cross the sky to set in the east in under four 
hours; the disc woiild subtend 12 min of arc (compared to 30 min for the Earth's 
Moon). This is aliiK>st half the apparent diameter of the Sun as seen from Mars and 
the penumbral shadow of Phobos sweeps across the martian surface on every orbit. 

As is characteristic of satellites within stationary orbit , the tidal drag on Phobos is 
decreasing its orbital distance from Mars and Phobos will strike the martin surface 
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TABLE 6-14. CHARACTERISTICS OF PHOBOS AND DEIMOS (6-68) 



Phobos 

Deimos 

Orbital Elements 



Semi-major axis 

9378 km 

23459 km 


2. 76 R.. 

Mars 

^Mars 

Eccentricity 

0.015 

0.00052 

Inclination (deg) 

1.02 

1.82 

Sidereal period 

7^3 9™ 13. 85® 

30^17“54.87® 

Physical Parameters 



Longest axis 

13. 5 km 

7.5 km 

Intermediate axis 

10. 7 km 

6. 0 km 

Shortest axis 

9. 6 km 

5. 5 km 

Rotation 

synchronous 

synchronous 


„ -3 

. -3 

Density 

2.0 g cm 

1. 9 g cm 

Mass 

9.8 X 10^^ g 

2.0 X 10^^ g 

Albedo 

0.05 

0.06 

Surface gravity 

-2 

1 cm sec 

0. 5 cm sec ^ 


8 

in less than 10 years [6- 69 through 6-71]. It is possible that tidal forces win result in 
disruption first as Phobos crosses the Roche radius. The surface of Phobos is heavily 
cratered. The two largest craters, "Stickney" and "HaH," are 10 and 5 km in diame- 
ter, respectively {6-72]. 

The orbit of Deimos lies just outside the stationary altitude. Its surface is also 
heavily cratered. There is evidence of regoliths on both satellites to depths of at 
least 5 m and probably more [6-72,6-73]. 

(Photos from the Viking Missions are available from NASA.) 

6.11 SPACE IN THE VICINITY OF MARS 

Mars' nearest celestial neighbor is the Earth-Moon system. Mars is near the 
asteroid belt, and interactions may have been an important factor in the evolution of 
the planet and its satellites. 
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6.12 THE SEARCH FOR LIFE ON MARS 

In addition to the lander cameras, which would be expected to show any large 
life-forms, the Viking landers contained thi*ee experiments to search for indications of 
life. 


a) T.he gas exchange experiment , designed to detect metabolization and related 
gaseous exchange products in the presence of a nutrient solution [6-42]. 

. b) The labeled release experiment designed to detect the release of radioactiv- 
ity from labeled carbon bound in a nutrient solution [6-74]. 

c) The pyrolitic release experiment which was based on the assumption that 
any martian life would be able to incorporate radioactivity-labeled atmospheric carbon 
in the presence of sunlight (i.e. , by photosynthesis) [6-751. 

All three experiments showed signs of activity, but probably of chemical not 
biological origin [6-44,6-76,6-77,6-78,6-79]. 

In addition, the Gas- Chromatograph/Mass Spectrometer (GCMS) failed to detect 
organic material at even very low concentrations , less than parts per billion for heavy 
organics and parts per million for lighter ones [6-41]. 

Although the resiolts of the biology detection experiments were hardly encourag- 
ing, all the chemical elements necessary to support Earth-type life exist on Mars, and 
speculation about life there continues [6-14,6-74,6-80,6-81,6-82,6-83]. 

The properties of the martian environment that seem most unfavorable to life as 
we know it , are : 

a) The scarcity of water and lack of liquid water at the surface 

b) The cold temperatures and extreme temperature variations that occur both 
diurnally and seasonally 

c) The solar ultraviolet radiation between 1900 and 3000 X, which penetrates to 
the surface (Table 6-12, paragraph 6. 6. 3. 1) 

d) The presence of strong oxidants in the soil 

e) The low surface atmospheric pressure 

f) The low atmospheric concentration of N 2 and the apparent absence of any 
form of nitrogen in the top-regolith. 

The issue will not be settled definitively until further studies are made, pref- 
erably through specimens returned to Earth or manned exploration. 


6,13 SPACECRAFT MISSIONS TO MARS 

Table 6-15 is a list of Mars missions. 
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TABLE 6-15, MARS MISSIONS 


Mission 

Origin 

r" — — 

Launch Date 

Outcome 

Mars 1 

USSR 

November 1, 1962 

- Contact lost after 105 million 
kilometers 

Mariner 3 

USA 

November 5, 1964 

Failed (shroud malfunction) 

Mariner 4 

USA 

November 28, 1964 

Fly-by on July 14, 1965 

Zond 2 

USSR 

November 30, 1964 

Batteries failed on May 5, 
1965 , 

Mariner 6 

USA 

February 24, 1969 

Fly-by on July 31, 1969 

Mariner 7 

USA 

March 27, 1969 

Fly-by on August 5, 1969 

Mars 2 

USSR 

May 19, 1971 

Capsule ejection on 
November 27, 1971 

Mars 3 

USSR 

May 28, 1971 

Capsule soft landing on 
December 2, 1971 

Mariner 9 

USA 

May 30, 1971 

Orbiter since November 13, 
1971 

Mars 4 

USSR 

July 21, 1973 

Orbiter since February 10, 
1974 

Mars 5 

USSR 

July 25, 1973 

Orbiter since February 12, 
1974 

Mars 6 

USSR 

August 5, 1973 

Landed on March 12, 1974 

Mars 7 

USSR 

August 9, 1973 

Contact lost on March 9, 
1974 

Viking 1 

USA 

August 20, 1975 

Landed on July 20, 1976 

Viking 2 

USA 

September 9, 1975 

Landed on September 3, 
1976 
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APPENDIX A, NEUTRAL ATMOSPHERE MODELS 


A.l DENSITY VARIATIONS AND MODELS 

Since the advent of the first orbiting satellite in 1957, significant advances 
have been made in our knowledge of the Earth's upper atmosphere. The accepted 
model of the Earth's upper atmosphere at that time was one in which the temperature, 
pressure, density, and composition varied only with height. Predictions of futixre 
sightings (orbital positions) of the first satellites using this invariant atmosphere 
soon showed that the model was in error. Subsequently, several periodic variations 
in the density, temperature, and composition of the upper atmosphere have been 
identified. They are classified as: 

a) Variations with the solar cycle approximately 11 years 

b) Variations with the daily change in activity on the solar disc 

c) The divirnal variation 

d) Variations with geomagnetic activity 

e) The semiannual variation 

f) Seasonal-latitudinal variations of the low thermosphere 

g) Seasonal-latitudinal variations of helium at all altitudes 

h) Rapid density fluctuations probably connected with tidal and gravity waves. 

The first seven have some regularity and can be modeled with varying degrees of 
accuracy. Since temperature plays a minor role in comparison with density, current 
models have been developed to represent, insofar as practical and possible, the 
variability of the ambient mass density rather than the temperature. The models 
are based on temperature profiles (Fig. A-1) [A-1] which have been adjusted so as 
to produce the density values derived from the analyses of satellite orbital decay 
data. 


A. 1.1 Variations with Solar Activity 

The ultraviolet solar radiation that heats and causes compositional changes 
(Fig. A- 2) [A-1] in the Earth's upper atmosphere consists of two components, one 
related to active regions on the solar disc and the other to the disc itself. The 
active-region component varies from day-to-day while the disc component varies 
more slowly, presumably with the longer periodicities in the solar activity; i.e., the 
approximately 11-year solar cycle. The atmosphere has been observed to react in 
a different manner to each of these two components. Jacchia and Slowey have found 
that the disc component of the solar radiation is , for all practical purposes , linearly 
related to the 10.7-cm solar flux smoothed over six solar rotations (162 days). 
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Figiare A-1. Four temperature profiles from the 
present models [A-1]. 

The effect of solar activity on atmospheric density is shown in Figures A- 3 and 
A-4 after Jacchia [A-1], When the short -period oscillations which are caused by the 
diurnal variation are removed, there is essentially an 11-year variation that parallels 
the smoothed 10.7-cm solar flux data. 

Figure A-5, after Jacchia [A-1], shows some of the variations in greater detail. 
The variations with a 2 7- day period are caused by the Sun's rotation - the active- 
region component. The semiannual variation can also be seen. 


A. 1.2 The Diurnal Variation 


Analyses of satellite orbital decay histories have shown that upper atmospheric 
densities reach a maximum around 2 p.m. local solar time at a latitaide approximately 
equal to that of the subsolar point while the minimum occurs between 3 and 4 a.m. 
at about the same latitude in the opposite hemisphere. Consistency between temper- 
ature and density cannot be achieved on a diurnal basis in a static model; therefore, 
the temperature profile has been used as a parameter so that observed density values 
can be reproduced by the models. The temperature distribution needed to achieve 
this required reproduction of the observed density values is depicted on a global 
scale in Figure A- 6 [A-1]. Even though the global temperature distribution is an 
artifact developed solely for use in the density and composition models, some experi- 
mental results are in good agreement [A- 2 through A-5]. Thomson- scatter temper- 
ature measurements [A- 6 through A- 8], as shown in Figure A- 7 [A-1], generally 


A-2 



GEOMETRIC HEIGHT (km) 


GEOMETRIC HEIGHT (km) 



GEOMETRIC HEIGHT (km) 



Figure A- 2. Atmospheric composition for three values of the exospheric 
temperature [A-1]. (Number densities (n) are plotted against 
geopotential height. The corresponding geometric heights 
are marked at the top of the diagrams.) 






Figure A- 3. Density profiles for seven 
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Figiore A- 4. Ten-day means of the densities obtained from the drag 
of the Explorer 1 satellite compared with variations in the 
10.7-cm solar flux [A-1]. [M.J.D. is the Modified 

Julian Day (J.D. minus 2 400 000.5).] 
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Figure A- 5. Densities obtained fixjm the drag of the Explorer 1 
satellite compared with variations in the 10.7-cm solar flxix. 
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Figure A- 7. Atmospheric temperatures obtained on November 30, 1967, 
by Carru and Waldteufel [A- 7] by use of Thomson-scatter techniques, 
compared with temperatvires predicted by the present models 
for a height of 300 km (solid line). 

show that the temperature maximum occurs between 3 and 5 p.m. rather that near 
2 p.m. This controversy has not been resolved at this time; however, it appears 
as if there is a phase lag between the density maximum and the temperature maximum 
which cannot be included in the current atmospheric models. 


A. 1.3 The Semiannual Variation 

No satisfactory explanation has been found for this variation. It was initially 
assumed that this density variation could be linked with a temperature variation; 
however, data from analyses of more recently orbited satellites showed that the 
original assumption was in error. The amplitude of this density variation is strongly 
height -dependent and variable from year-to-year, with a primary minimum in July 
and principal maximum in October; however, it does not appear to be related to solar 
activity. Figure A- 8 [A-1] shows the semiannual variation as derived from the orbital 
decay analysis of the Explorer 32. 


SATELLITE 1966 44A (EXPLORER 32) 
M J O. 



figure A- 8. The semiannual density variation as derived from drag analysis 
on Satellite 1966 44 A (Explorer 32). [All other variations have been 
suppressed by using the appropriate equations. M.J.D. is the 
Modified Julian Day (J.D. minus 2 400 000.5).] 
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A. 1.4 Variations with Geomagnetic ActiAdty 


Analyses of orbital decay histories can give only a blurred picture of the com- 
plex reaction of the upper atmosphere to geomagnetic disturbances (Fig. A- 9) [A-1]. 
Blamont and Luton [A- 9] and DeVries [A- 10] have shown that the upper atmosphere 
first reacts in the auroral zones with the energy subsequently propagating toward 
the equator apparently in the form of wave-like perturbations. It appears as if the 
atmosphere reacts with a zero time delay in the auroral zones with the geomagnetic 
storm effects showing up in the equatorial zone about 6 to 9 hr later. It is very 
difficult to adequately include this effect in a sthtic model. The current models 
calculate the density variations on the basis of a global increase in exospheric temper- 
ature. Observations [A- 9, A- 10] have shown that this is not the case; however, for 
satellite lifetime prediction calcvdations , an assumed global temperature increase is 
acceptable. For some calculations, such as control dynamics analyses, control moment 
gyro analyses and aerodynamic torques, instantaneous temperature (and therefore 
density), increases for short time periods, and specific locations may be required. 
These studies may require special applications of the models. 


A. 1.5 Seasonal-Latitudinal Variations of the Lower Thermosphere 


Presently accepted models assume constant temperature and density values at 
90 km to prevent the models from becoming too complex even though large temper- 
ature and somewhat smaller density variations are known to exist. Current models 
therefore are constructed with a seasonal-latitudinal density variation which varies in 
the vertical from 0 at 90 km to a maximum at 110 and back 0 at 170 km. In the 
horizontal, the maximum occurs on December 27 at the North Pole and on June 25 
at the South Pole. These variations are included as additions to the mass densities 
calculated as a function of the exospheric temperature \ased in the model. These 
variations are small and because they occur solely below 170-km altitude they will 
have little effect on orbital lifetime prediction calculations. 


A. 1.6 Seasonal-Latitudinal Variations of Helium 


Experimental results have shown a strong increase of helium above the winter 
pole. The mechanism for this migration is unclear; however, empirical equations 
which describe the phenomenon are included in the models. These variations influence 
the computed densities only at heights above approximately 600 km (Fig. A- 10). 


A, 1.7 Density Waves 

Ambient density waves have been detected throughout the upper atmosphere in 
the height range from 140 to at least 510 km [A- 4, A- 10]. In addition, traveling 
ionospheric disturbances (TIDs) have long been thought of as manifestations of 
internal gravity waves. The waves appear to be more prevalent at the higher lati- 
tudes near the auroral zone; however, TIDs have been observed at aU latitudes. 
Their vertical wavelengths apparently increase with altitude. Density increases on 
the order of 100 percent have been observed to occur over short distances [A- 10]. 
The waves apparently propagate from either south to north or north to south with 
maximum horizontal wavelengths on the order of 500 to 700 km. Although current 
models do not include variations associated with internally propagating waves, users 
should be cautioned that they have been observed. 
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Figure A- 9. The geomagnetic effect as derived from the drag of six 
satellites in May and Jxme 1967. [The plotted ATs are temperature 
residuals from the models when all variations except the geomagnetic 
effect have been taken into account. M.J.D. is the Modified 
Jxilian Day (J.D. minus 2 400 000.5).] 
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Figure A- 10. Observed and computed density variations caused by the helium 
migration, as derived from the drag on Satellite 1963 53A (Explorer 19). 
[The data points in A i log p are density residuals from the model when 
all variations except the helium variations are suppressed. A 2 log p 
is the difference 0-C in Ai log p. M.J.D. is the Modified 
Julian Day (J.D. minus 2 400 000.5).] 


A. 2 APPLICATIONS 

Requirements exist for models of the Earth’s neutral upper atmosphere which 

win : 


a) When used in lifetime prediction calculations , make predictions of lifetimes 
for satellites with orbital altitudes between approximately 180 and 2500 km within 
reasonable limits. 

b) Contain the short period fluctuations observed in upper atmospheric density 
which may influence the design and operation of reaction control systems (RCS), 
thrust augmented control systems (TACS) , control moment gyroscopes (CMC) , rendez- 
vous and docking systems, etc. 

c) Provide continuous values of atmospheric parameters from orbital altitudes 
down to about 25-km altitude for use in evaluating thermal protection systems (TPS) , 
control and guidance systems, etc., during reentry of space vehicle (e.g. , the 
Shuttle), and in accomplishing orbital debris analyses. 

d) Provide data on the composition and thermal structure of the upper atmos- 
phere and its temporal and spatial variability for use in interpreting and analyzing 
results of scientific and technical investigations; e.g., contamination surrounding ^e 
Spacelab. 

e) Provide data on the variability of the neutral atmosphere at mesopause 
heights for use in aerobraking maneuver studies. 
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A. 3 CRITERIA 

a) The Global Reference Atmosphere (A-13, A-14) model will be used in conjunc- 
tion with the MSEC Lifetime Prediction Computer Routine for calculating orbital life- 
time. This tirodel is a combination of the 4-D Model Atmosphere, the Groves Model, 
and the Jacchia 1970 (A-11) Atmosphere model. It provides continuous density, tem- 
perature, and composition data from the surface to 2500 km for analyses requiring 
these data. Copies of the computer program for this model and a computerized simpli- 
fied version of the Jacchia 1970 atmosphere model are available upon request from the 
Atmospheric Sciences Division, Marshall Space Flight Center, Alabama 35812, 

b) Personnel performing analyses requiring knowledge of any small-scale, short- 
time period fluctuations in atmospheric parameters such as those perturbations asso- 
ciated with internally propagating waves should contact the Atmospheric Sciences Divi 
Sion, Marshall Space Flight Center, Alabama 35812. 
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